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EISENHOWER AND GRASS RIVER LOCK MODELS* 


Martin E. Nelson,! M. ASCE, and 
Harvey J. Johnson,2 A.M. ASCE 
(Proc. Paper 1582) 


SYNOPSIS 


Models of the St. Lawrence Seaway Dwight D. Eisenhower and Grass River 
locks were constructed to a scale of 1 to 24.24, model to prototype. A type of 
hydraulic system consisting of intake ports in the upper gate sill, culverts 
and ports in both walls of the lock and lateral culvert diffusers in the lower 
approach channel was selected prior to model studies. Vortical activity over 
the intake ports and turbulence in the lock chamber and lower approach were 
reduced to acceptable proportions. Culvert and gate pressures, emergency 
gate and tainter valve operating loads, and hawser forces were measured. 


Prototype operation procedures were selected from results of the model 
tests. 


INTRODUCTION 


Model studies of the Dwight D. Eisenhower and Grass River locks, St. 
Lawrence Seaway, have been conducted by the U. S. Army Engineer District, 
St. Paul, in facilities provided by the University of Minnesota at the St. 
Anthony Falls Hydraulic Laboratory. The locks and other features of the 
American portion of the seaway were designed for the St. Lawrence Seaway 
Development Corporation by the U. S. Army Engineer District, Buffalo. Dur- 
ing the course of the model tests several conferences were held at the 
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laboratory with representatives of the Corporation; Office, Chief of Engi- 
neers; U. S. Army Engineer Division, North Central; and U. S. Army Engi- 
neer District, Buffalo. Mutually helpful conferences were also held by 
Corps laboratory personnel and engineers of the St. Lawrence Seaway Au- 
thority of Canada. 


Criteria for Design of Lock Hydraulic Systems 


The obvious function of a navigation lock is to facilitate the transit of 
water-borne traffic past abrupt changes in elevation of navigable rivers or ‘ 
canals. The design engineer holds as his primary objective a structure that 
will perform this function rapidly, safely and economically, insofar as these 
characteristics are compatible and mutually attainable. The key is found in 
the design of a sound and efficient filling and emptying system for the lock. 
The first criterion by which the performance of the hydraulic system is 
gaged is the time required to elevate a vessel from the lower pool level to 
the surface of the upper pool or the reverse of this process in the case of 
downbound traffic. Increasingly greater emphasis is being placed-on the im- 
portance of expediting river traffic along navigation channels and the avoid- 
ance of unnecessary delays in lockage. A monetary value of time gained or 
sacrificed in lockage, based on the operating costs of towing equipment and 
the value of cargo, is sometimes used to determine the amount of capital in- 
vestment that would be justified to achieve a more efficient lock operation. 
The second criterion requires that the forces acting on a vessel during a 
lockage operation shall not produce excessive stresses in the mooring lines 
nor endanger the safety of the craft or any part of the structure. This cri- 
terion is to some extent in conflict with the first. Obviously, as the capacity 
of the hydraulic system increases the amount of energy to be dissipated in- 
creases and more elaborate devices are required to reduce to satisfactory 
proportions the attendant disturbances in the lock chamber. The evidences 
of unsatisfactory hydrodynamic performance of a lock filling and emptying 
system are manifested in excessive strains in the mooring lines due to 
strong currents or surging of the water level in the lock chamber and exces- 
sive turbulence. Hawser forces are most critical with respect to large ves- 
sels while unrestrained turbulence can cause hazardous conditions for smal- 
ler craft. The third criterion to be held in proper perspective together with 
speed and safety is that of economy. In his efforts to achieve a high degree 
of efficiency and safety in the operation of the hydraulic system the designing 
engineer must not lose sight of the fact that construction costs generally in- 
crease in proportion to the refinement in hydraulic design. Consequently, 
the importance of the respective criteria must be weighed by the designer e 
whose judgment and experience must to a large extent enter into the selection 

of the optimum type of hydraulic system for a specific lock project. a 


Selection of Hydraulic System Type 


Lock filling and emptying systems in use today may be classified into 
three general types: (1) End filling systems admit the lockage water near 
the upstream end of the lock chamber by means of loop culverts around the 
service gates, sluices through the sill, or by operation of the gates. A simi- 
lar arrangement at the downstream end provides the means for emptying the 
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lock. (2) The most common type of hydraulic system used in locks of nomi- 
nal height consists of culverts in the lock walls with intake ports in the upper 
pool and short ports to distribute flow into the lock chamber. In the empty- 
ing operation the flow in the chamber ports is reversed and discharges into 
the lower pool through port manifolds or discharge basins. (3) Several locks 
have been constructed in recent years with filling systems that admit the 
lockage water into the lock chamber through openings in the floor. The sup- 
ply conduits with intake manifolds in the upper pool are located either in the 
lock walls or in the floor. From the culverts the flow enters the chamber 
through a system of lateral and port manifolds or directly through vertical 
ports as in the case of longitudinal floor culverts. The emptying operation 
for this type of system is performed with any of the various types used with 
wall port hydraulic systems. 

The end filling type of hydraulic system was not considered practical for 
the seaway locks because prior experience had indicated that satisfactory 
lockage conditions could not be attained in locks of the size in question with- 
out using excessively long periods of time to fill and empty the chambers. 
The floor lateral system undoubtedly would have met all the design criteria 
and provided satisfactory lockage conditions for the Eisenhower and Grass 
River locks. However, the wall port system generally can be constructed at 
a substantial saving in cost and the fortunate circumstance of a deep water 
cushion made available by the necessity for removal of overburden to a con- 
siderable depth in order to obtain adequate bearing for the lock walls made it 
appear possible to satisfy the design criteria with this type of system. In 
figure 1 are shown data obtained from several tests in prototype locks with 
wall port hydraulic systems which provide a guide as to permissible rates of 
filling under various depths of water cushion in the lock chamber. These 
data, which have been used in preliminary design of several new locks, indi- 
cated that the seaway locks with wall-port hydraulic systems could be ex- 
pected to operate satisfactorily. 

Model tests on locks have repeatedly indicated that hydrodynamic forces 
exerted upon craft in the chamber during filling or emptying operations are 
due to swelling of the water surface along the lock axis. Swells or surges are 
set up at the very beginning of the valve opening period when the initial pres- 
sure wave is transmitted along the culvert and into the lock chamber through 
the connecting ports or laterals. The degree of surging and consequently the 
magnitude of force imposed upon a craft in the lock chamber reaches a mini- 
mum value when the manifold system of ports or laterals occupies the ap- 
proximate middle third of the lock chamber walls. This feature was approxi- 
mated in the design of the seaway locks. 


General Description and Purpose 


The Eisenhower and Grass River locks located in the Long Sault Canal 
provide a means for navigation to pass the Long Sault Dam and power plant. 
The locks are 80 feet wide and 860 feet long with a 30-foot depth over the gate 
sills. The Eisenhower lock has a normal lift of 43 feet and the Grass River 
lock located downstream is designed for a normal lift of 49 feet. The inter- 
mediate channel or pool between the two locks is about 3 miles long and has a 
navigable width of 442 feet. 


A vertical lift gate is located in the upper end of the Eisenhower lock as an 
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emergency closure structure. The service gates of both locks are of the 
two-leaf mitering type. The hydraulic systems for the locks are almost 
identical consisting of a culvert in each wall drawing water from the upper 
approach channel through ports in the gate sill, filling or emptying the cham- 
ber through 20 ports in each wall and discharging water from the locks 
through lateral culvert diffusers in the floor of the lower approach channel. 

The model studies described in this paper were conducted for the purpose 
of developing filling and emptying systems for the Eisenhower and Grass 
River locks that will perform the lockage operations as rapidly as possible 
without violating the requirements for safety of vessels and structure. The 
criteria adopted for the design of the hydraulic systems specified a lockage 
time of about 8 minutes with a maximum stress not to exceed 10 tons in the 
mooring lines of a vessel of any size and a degree of turbulence that would 
not jeopardize small craft in the chambers or in the entrances to the locks. 
The scope of the model studies was broadened during the course of the pro- 
ject design stage to obtain data to aid in the design of the vertical lift emer- 
gency gate for the Eisenhower lock, the purpose of which is to prevent loss 
of the power pool in the event of damage to the lock service gates. Other 
incidental observations were made in the model to determine whether operat- 
ing or maintenance problems in the nature of cavitation, vibration, or erosion 
might be indicated and to measure the flood discharge capacity of the hy- 
draulic systems. 

The similarity of the two locks made it possible to conduct the tests in one 
composite lock model making minor modifications to convert from physical 
properties of one lock to those of the other. The lock models were construct- 
ed to a scale 1/24.24 of the proposed prototypes. However, dimensions and 
quantities used in this paper refer to the prototype. 


Intakes 


The intakes in the sill for the wall culverts, figure 2, are two manifolds, 
each consisting of three ports, 9 feet 6 inches by 22 feet at the top converging 
to a constricted throat area. Uniformity of flow among ports was accom- 
plished by making the constricted area progressively smaller as the flow 
distances diminished. The maximum variation from uniformity of flow in the 
intake ports was 2 percent. 

Vortices over the intakes presented a problem in the first model when the 
port dimensions were only 8 feet by 12 feet. The proximity of the miter gate 
and gate recesses in the walls was considered a contributing factor but satis- 
factory conditions were obtained by doubling the port area and enlarging the 
constricted sections. This reduced the maximum port velocity from 21 feet 
per second in the first model to 12 feet per second in the final model. All 
pressures measured in the water passages were found to be positive. 


Culvert Bends 


Flow from the intake manifolds passes through a 90-degree bend in each 
culvert. Two piezometers are located at the midpoint of each bend in oppo- 
site walls of the culvert to measure pressures and serve incidentally as flow 
metering devices for the purpose of comparing calibrations of the model and 

prototype. Pressures measured for the adopted filling operation were always 
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in the positive range at both piezometers. The formula derived for the flow 
meters was: 


Q = 1230x _/(pressure difference between opposite piezometers, c.{.s. 
in feet of water) 


Valves 


The valves which control the flow in the filling and emptying culverts are 
of the reversed tainter type, that is, with the trunnions upstream of the skin 
plate. The linkages between constant speed motors and the valves produce an 

. opening pattern that is not a straight line function of time; the valve opening 
rate increases slightly with respect to time. The culverts at the valves are 
12 feet wide and 14 feet high. Stop log recesses located about 20 feet up- 
stream and 30 feet downstream from the valves are provided to facilitate 
dewatering for inspection and maintenance. 

Culvert pressures downstream from the valves momentarily reached a 
slightly negative pressure at the halfway point of the valve opening cycle. 
The stop log recess in that region will be sealed to prevent the possibility of 
entraining air in the flow which may accumulate in the culvert and then be re- 
leased suddenly into the chamber with such force as to endanger small craft. 


Chamber Manifold 


Between the valves and the lock ports the culvert is enlarged from 12 feet 
by 14 feet to 15 feet by 16 feet. 

Water enters the lock through 20 venturi-shaped ports in each wall which 
are 3 feet 1 inch wide at the throat, 3 feet 7 inches wide at the wall face, and 
4 feet high and 15 feet long. All the ports are sloped downward 1 on 15 and 
the upstream 10 ports in each wall are turned 3.5 degrees in the upstream 
direction, the remainder being perpendicular to the lock wall. They are sub- 
merged about 34 feet below normal tail water and are about 8 feet above 
natural rock floor in the Grass River lock and about 26 feet above natural 
rock floor in the Eisenhower lock. The average jet velocity in the ports 
reaches a maximum of about 20 feet per second. 

The limitation of surface turbulence was difficult in the Grass River lock 
because of the relatively shallow water depth in the chamber when the inflow 
during filling is the greatest. This condition is shown graphically on figure 1 
where the solid line is an approximate envelop curve of acceptable turbulence 
for the chamber water depth and rate of filling in a number of existing locks. 
Note that the rate of filling for the Grass River lock exceeds that given by the 
envelop curve. This occurs for about 3 minutes of the 8 1/2-minute filling 
time. Experimentation in the model resulted in a design that reduced the 
’ turbulence to acceptable proportions. 

The ports were first tested in a single port model and then a refinement of 
design was accomplished in the composite lock model of the Grass River lock. 
The basic designs tested are shown in figure 3. The 1 on 15 port slope, flare 
of the walls, number and size of ports were selected as giving the best water 
surface conditions. The 3.5-degree horizontal angle was adopted because this 
feature indicated a reduction in hawser pulls. For stability reasons the lock 
wall was designed with a toe projecting into the chamber near the base of the 
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wall. A sloping toe was sufficient to satisfy the structural design but the 
model tests indicated a slight reduction in turbulence during filling opera- 
tions by building the projection in steps. 

A great number of baffles, deflectors and port shapes were tested but the 
only modification which performed as well as the design developed was an 
interim test port 5 feet 2 inches square at the large end and 23 feet long, 
projecting horizontally 8 feet beyond the lock wall. The jet issuing from this 
port entered the lock on a slope of 1 on 15. The same port length to width 
ratio and jet slope were approximated in the final port design. 


Discharge Manifold 


The emptying system is a continuation of each wall culvert of the filling 
system which terminates in a manifold of three smaller lateral culverts 11 
feet 6 inches high tapered from 6 feet 4 inches wide to 2 feet in width which 
diffuse flow in the lower approach channel through a total of 96 ports 4 feet 
3 inches high and 1 foot 8 inches wide. The lateral culvert walls were 2 feet 
thick. Thus, the ports were 2 feet long but a remarkable improvement in 
water surface conditions was obtained by doubling the length of the ports by 
means of projections on the culvert walls. Velocities along the bottom of the 
lower approach channel reached a maximum of about 5 1/2 feet per second. 


Lockage Operations 


Model hulls of two ships were used in this study. The George S. Humphrey 
is considered one of the largest carriers on the Great Lakes with a total 
length of 710 feet and a 75-foot beam, representing a gross displacement of 
37,000 tons when loaded to a maximum draft of 27 feet. The George Eaton is 
a smaller ship, length 257 feet, beam 43 feet, draft 16.5 feet, and displace- 
ment 5,000 tons. The forces required to hold these ships at various locations 
in the model lock chamber were recorded for a wide range of lifts and valve 
operations. 

A recorder built especially for this type of testing made a continuous 
record of water surface levels, valve movement, surging, and hawser pulls. 
These data are shown on figures 4 to 7 for the selected operation of the two 
locks. 

To avoid exceeding a maximum longitudinal hawser pull of 10 tons by the 
largest ship a valve opening time of about 4 minutes was selected for filling 
and 2 minutes for emptying. A lifting force of about 20 tons more than fric- 
tion and weight will be required to raise the valves. Vibration of the model 
valve occurred during the first foot of opening and reached a peak frequency 
of 1/5 cycle per second producing variations in loading of between 2 tons and 
12 tons, in prototype values. 

The Eisenhower lock will fill in 8.1 minutes with a valve opening time of 
4 minutes and will empty in 8.0 minutes with a valve opening of 2 minutes, 
while operating at a maximum lift of 43 feet. 

The Grass River lock will fill in 8.5 minutes when the valves are opened 
in 4 minutes and will empty in 8.3 minutes when the valves are opened in 2 
minutes, while operating at the maximum lift of 49 feet. 
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Emergency Gate 


The emergency gate in the upper approach of the Eisenhower lock, shown 
in figures 2 and 8, consists of a skin plate supported on the downstream side 
by horizontal plate girders spanning the 80-foot channel. The bearings of the 
gate roll on the downstream side of a recess in each wall and the gate when 
open will be submerged in a pit in the bottom of the channel. 

The gate model simulated 50 feet of the 80-foot channel width to a scale of 
3 to 100. 

The gate as designed with a sharp crest functioned satisfactorily. The 
force required to raise it varied from 364 tons, the submerged weight while 
resting in the recess of the channel bottom, to 880 tons near the top of the 
closure operation. The three conditions of flow and total lifting forces are 
illustrated in figure 8. The nappe becomes aerated as the crest reaches 
elevation 226 when air is drawn in through the wall recesses. The nappe 
falls on the top plate girder when the crest is near elevation 239. The loss 
of buoyancy when the nappe becomes aerated and the weight of water trapped 
on the webs of the girders are also shown on figure 8. No attempt was made 
to simulate prototype friction factors in the model and consequently the 
friction loads indicated by the model are excessive relative to those antici- 
pated in the prototype. 
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SYNOPSIS 


During the design of the Tennessee Valley Authority’s Kingston Steam 
Plant it became desirable to develop facilities which would divert cold waters 
which pass the site as thermally stratified density underflows and cause them 
to be drawn approximately 2 miles up a side channel at the same depth as the 
original channel, after which they would cross an area one mile long and 20 
to 25 feet higher than the original channel. The unique solution to this prob- 
lem consisted of a canal across the high ground area at the end of which was 
placed a skimmer wall structure which allowed only the cold waters in the 
bottom of the side channel to flow into the canal. A diversion dam was also 
placed underwater immediately downstream from the side channel. This pa- 
per gives the methods used in determining the required design. Much of the 
study concerned density current flows where analytical methods were not 
available and had to be developed. The efficacy of these methods and the 
adopted design is proven by the presentation of field operational data obtained 
after the structures were placed in use. The financial benefits attributable 

to the design were found to be high. 


INTRODUCTION 


During the construction of the Kingston Steam Plant, TVA engineers were 
faced with a challenging problem in developing a design which would make 
available for condensing purposes cold waters which passed the plant site as 
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thermal density underflows. The general site plan for the Kingston Steam 
Plant is shown in Figure 1. The plant is located on the neck of a peninsula 
formed by the Emory and Clinch River embayments of the Watts Bar Reser- 
voir. During the summer months the only appreciable flows in the reservoir 
occur as thermal density underflows in the Clinch River arm.(1) The source 
of these underflowing waters is TVA’s Norris Reservoir some 78 miles up- 
stream. Because of stratification in this large storage reservoir, relatively 
cold water is discharged through the powerhouse turbines from spring until 
late fall. Over-all steam plant considerations required that the condensing 
waters be drawn from the north or Swan Pond side of the peninsula and dis- 
charged on the south or Clinch River side. Under normal conditions the 
Emory River flows are negligable during the summer months. As a conse- 
quence the cold Clinch underflows push their way up the original Emory River 
channel. This fact is shown by Figure 2, which presents the maximum tem- 
peratures observed at Clinch River Mile 4.0 and Emory River Mile 2 over a 
5 1/2-year period. Except for the surface the two sets of curves are quite 
Similar, indicating that the Clinch River waters penetrate and thus control 
the Emory River temperatures in the immediate area of the Kingston Plant. 
Figure 2 also shows that in the Swan Pond Embayment the surface tempera- 
tures become very hot, reaching 95° F. on occasions. 

Figure 3 shows a section taken along the southern edge of Swan Pond at 
Section A-A of Figure 1. This section is typical for Swan Pond and shows 
that the water depth in this area is comparatively shallow and the water is 
thus susceptible to stagnation and marked temperature increases due to solar 
heating. In its natural state, the high bottom elevation of this area in essence 
forms a dam which would block out the cold underflowing Clinch River water. 
Under this condition, only the hot surface waters would be drawn into the 
pumps. 

For the best over-all plant design it was desirable to set the maximum 
allowable condenser water temperature at 75° F. The plant as originally 
conceived was to have a maximum condenser water demand of 2470 cubic 
feet per second. This value was used for all underflow diversion design pur- 
poses. Later changes reduced the maximum demand to 2310 cubic feet per 
second. There was no question about the normal Clinch flows being sufficient 
and of the proper temperature to meet the plant demand. The problem was 
thus one of developing a design which would economically make available to 
the pumps the cold underflowing Clinch River waters by bringing these up the 
Emory River channel and across Swan Pond to the condenser water pump in- 
takes. 


Proposed Design 


One obvious solution to the problem would have been to excavate a channel 
across Swan Pond so that the cold Clinch waters could be reached. This solu 
tion was not economical since over 1 1/2 million cubic yards of material, a 
high percentage of which would have been rock, would have had to be exca- 
vated under water. It was proposed, therefore, to provide an intake canal 
across the high Swan Pond area and into the original Emory River channel by 
means of dikes such as shown in Figure 1. A wall structure, as shown on 
Figures 1 and 3, would close the east end of the canal except for a series of 
openings along its bottom. These openings would allow the cold bottom water 
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to flow into the canal to replace water pumped from the canal for condensing 
purposes. In effect the wall would skim off the hot top waters and it was 
therefore termed a skimmer wall. Figure 4 shows the details of the final 
wall design. Figure 5 is a photograph of the finished structures. 

In addition, it was desired to investigate the feasibility of an underwater 
dam located on the Clinch River just downstream from the Emory River 
mouth. Since the releases from Norris Dam are intermittent and subject to 
extreme regulation,(1) especially during the summer months, it appeared 
that such a dam might be required to divert the Clinch flow into the Emory 
channel at a high enough level to make full use of the available cold water. 


Design Considerations 


The design required three separate but interrelated studies: (1) The study 
of flow conditions in the Emory River channel from its mouth to the skimmer 
wall, (2) the study of the diversion of underflow waters into the Emory from 
the Clinch, both with and without an underwater dam, and (3) the determina- 
tion of the basic dimensions for the skimmer wall openings. The first study 
determined the capability of the Emory River channel, and when combined 
with the second study it determined the anticipated temperatures. The 
second study also established the effectiveness of an underwater dam. 


Flow in Emory River Channel 


At the time of these studies, TVA had a large number of field measure - 
ments of the flow of thermal density currents obtained in Watts Bar and Fort 
Loudoun Reservoirs. An analysis of these data indicated that, in reservoir 
channels such as those at the project site and for thermal stratifications of 
the magnitude expected at the project site, a relationship exists between the 
maximum velocity at which density underflow can occur and the depth of the 
underflow. Figure 6 shows the data and the envelope curve established for 
use in the Kingston study. This curve indicates that, if at any given depth the 
velocity is equal to or less than the value given by the curve, thermal under- 
flow can be sustained. 

A second item of data available from the TVA field measurements was 
that, at water intakes, the velocity in the intake could be double the underflow 
velocity without breaking the stratification and pulling in warmer overlying 
waters. Velocities greater than double appeared to break the stratification. 

By using Figure 6 the capability of the Emory River channel was deter- 
mined in a preliminary manner. The section of minimum cross-sectional 
area in the channel was determined and a curve of depth versus velocity was 
plotted for the maximum condenser demand. This curve, cross-plotted on 
Figure 6, showed that a depth of 23 or more feet would be required to sustain 
a density underflow and the velocity at this depth would be 0.35 foot per 
second. The river channel bottom at the point of minimum cross-sectional 
area was at elevation 704. Therefore the maximum requirement of water 
would be drawn from below elevation 727. Reference to the Clinch River 
Mile 4 temperature curves of Figure 2 makes it apparent that an average 
temperature of 75° F. or less would always exist throughout the depth be- 
tween elevations 704 and 727. 

A more comprehensive survey of the capability of the Emory River chan- 
nel was made after preliminary dimensions for the skimmer wall openings 
were established. In this study the density underflow interface profiles were 
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computed by a modified standard step backwater computation method. A 
number of authorities have shown that stratified flow calculations can be 
made by using normal hydraulic equations if the gravitational term for the 


underflow is multiplied by the factor A where Y is the density of the under 


flow fluid and AY is the difference in density between the underflow and the 
overlying layer. 

For the Emory River channel backwater computations the Chezy equation 
was used for computing head loss for each reach. The equation was of the 


form: 
V=CyVRS 


Chezy’s C can be written as 


therefore C,, = Cy= 8g =C 


f 


In Chezy’s equation the hydraulic radius R is computed by using the cross- 
sectional area of flow and the wetted perimeter. For these calculations the 
total length around the flow section, including the width at the interface of the 
density underflow and the overlying water, was taken as the wetted perimeter. 
S is the slope of the energy gradient. The coefficient f is the Darcy-Weisbach 
friction coefficient for the channel reach under consideration. 

The field data from Watts Bar and Fort Loudoun Reservoirs indicated that 
an f value of about 0.016 should be used. This gives a Chezy C of 130 and for 
a a of 0.0016, which is the value expected in midsummer, a C,, value of 5.05. 
Subsequent operations have proven these values to be reasonable. 

These backwater calculations determined the total depth of the underflow 
at the confluence of the Emory and Clinch Rivers which was required before 
the diversion at the Emory River mouth could be studied. 


Division of Flow at Mouth of Emory River 


To determine the anticipated temperature of the underflow the amount of 
flow originating from the Clinch River underflow and the amount originating 
from the warmer overlying waters had to be determined. These amounts 
were found by considering the flow at the junction of the Emory River and the 
Clinch River to be similar to the branched flow problems in pipes. As a first 
approximation of the division of the density underflow, the discharge up the 
Emory River was assumed. Thus, the interface elevation at the confluence of 
the two rivers was obtained from the modified backwater curve computations. 
It was then assumed that the depth of underflow thus determined at the 
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confluence was the depth at which the underflow continued down the Clinch 
River. It was further assumed that the flow in the lower Clinch River was 
normal flow and the calculated interface depth was the normal depth. With 
the normal depth known, the discharge in the lower Clinch River could be de- 
termined. The density underflow in the upper Clinch River was then the sum 
of the two parts. The normal depth of the density underflow was computed 
from the equation 


3/ 2Q2 
SCy* we 


where W is the width of the density underflow interface. Since the width-to- 
depth ratio of the underflow channel is relatively large, the wetted perimeter 
was taken as two times the width of flow at the interface. This gave the num- 
ber 2 in the above equation. 

The effect of a submerged dam located just downstream from the mouth of 
the Emory River upon the division of flow was calculated by assuming that 
critical depth of flow‘would occur over the weir. The critical depth was de- 
termined from the equation 


2 
= 


where q is the flow per unit width. 

Figures 7 and 8 show schematized flow and temperature conditions and the 
results obtained for two computations. Figure 7 is for conditions with only 
the skimmer wall constructed, while Figure 8 shows conditions with the un- 
derwater dam added. The dam for these calculations had its top at elevation 
722, the value selected for construction. The percentage of the maximum 
condenser requirement that any particular Clinch River underflow could sup- 
ply for the two design conditions is shown in Figure 9. These curves show 
clearly the need for the underwater dam, especially when it is recognized 
that Norris turbines normally discharge an average of about 5000 cubic feet 
per second during the summer months. 

With the division of flow computed, temperature estimates were made 
based upon the curves of Figure 2. Figures 10 (a) and 10 (b) show the antici- 
pated inflow temperatures for Clinch River underflows of 1500 and 5800 cubic 
feet per second, respectively, with the skimmer alone and with both skimmer 
and underwater dam. These plottings show that at a discharge of 1500 cubic 
feet per second the submerged dam would reduce the temperatures obtainable 
from only the skimmer structure from 1 to 2 1/2 degrees Fahrenheit. Ata 
discharge of 5800 cubic feet per second the dam similarly would effect a re- 
duction of from 1 to 3 1/2 degrees Fahrenheit. Figure 10 (c) further shows 
the benefits from the submerged dam in making available at the intake mini- 
mum water temperatures at Clinch River discharges much lower than would 
be possible without the dam. With no dam it would require 9500 cubic feet 
per second of flow in Clinch River to give the same minimum temperatures 
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which, with the dam, are produced with Clinch River flows of only 4400 cubic 
feet per second—less than half of the quantity required without the dam. The 
higher discharges are above the normal operating range for Norris Dam, 
whereas the required 4400 cubic feet per second of flow with the dam is less 
than the usual operating discharges. Thus the dam makes it possible to 
achieve minimum intake temperatures that would not otherwise be obtainable. 
The underwater dam also provides a material benefit with respect to in- 
termittent releases from Norris. Daily fluctuations of the Norris discharges 
are not noticeable at the general location of the Kingston Steam Plant, but 
shutdowns of two days which occur normally over weekends result in the 
minimum water temperature rising for several hours in midsummer to ap- 
proximately 78° F. at the elevation of the skimmer wall intake. Longer 
shutdowns result in higher temperatures for longer periods. The underwater 
dam, however, will provide a storage pool of cold water which will require 
approximately 2.5 days to deplete. Thus, for a weekend shutdown condition 
the temperature rise would occur gradually throughout the period of draft and 
reach the 78° F. maximum only after 2.5 days. As a consequence the storage 


pool will materially reduce the effect of fluctuation in discharge from the 
Norris powerhouse. 


Skimmer Wall Design 


The basic dimensions for the skimmer wall openings were determined by 
first computing the area of opening required in order not to exceed the allow- 
able velocity and then determining the permissible maximum height of open- 
ing. The Watts Bar-Fort Loudoun underflow data indicated that the maximum 
average velocity that could be sustained in the Emory River channel immedi- 
ately upstream from the skimmer would be 0.35 foot per second without 
breaking the stratification. These data also indicated that the velocity 
through the intake openings could be 0.70 foot per second without disturbing 
the stratification. The necessary area of opening was determined, therefore, 
to be 3530 square feet. 

Physical characteristics at the wall site indicated that the bottom of the 
opening could be at elevation 700 and that a structure about 350 feet long 
could be used. It was considered desirable to keep the height of opening as 
small as possible to minimize the opportunity for the overlying waters to be 
drawn through the skimmer. Conversely, of course, it was desirable to keep 
the length of the structure as short as possible to minimize the cost. In co- 
ordination with the design engineers, a structure consisting of five bays each 
48 feet wide and 15 feet high with circular sheet piling piers 17 feet in diame- 
ter was adopted for final study. This is the design shown in Figure 4 and is 
the design adopted for construction. The backwater curves indicated that the 
15-foot high opening would be satisfactory but that any increase in opening 
height would result in some increase in the expected temperatures. 

Final check of the design was obtained by model studies which were spon- 
sored by TVA at the Massachusetts Institute of Technology, Hydrodynamics 
Laboratory. The results of these tests are being published concurrent with 
this paper in a paper by Harleman, Gooch, and Ippen, entitled “Submerged 
Sluice Control of Stratified Flow.” Upon the basis of the results from the 
M.I.T. tests, the maximum allowable velocity through the proposed opening 


would be 0.59 foot per second in the early spring and late fall when ar values 
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of 0.0005 are experienced and 1.06 feet per second in midsummer when the 


AY 


y value may reach 0.0016. These values generally confirm the satisfactori- 


ness of the adopted skimmer wall design. 

This design provides a flow area of 3600 square feet and the maximum 
plant demand is 2310 cubic feet per second. With the allowable velocity of 
0.59 foot per second from the M.I.T. tests for spring and fall, the flow 
through the skimmer wall opening without drawing in water from the overly- 
ing warmer water would be 2120 cubic feet per second. Since the maximum 
plant demand is 190 cubic feet per second more than this, this deficit over 
and above that available from the underflow would have to be supplied from 
the overlying water. In the early spring and late fall seasons this is not im- 
portant because the general over-all water temperatures during those periods 
are comparatively low and no significant loss in plant efficiency would result 
from the small amount of slightly warmer water being drawn through the 
skimmer into the intake canal. 

The velocity of 0.70 foot per second used in the design is far below the 
value of the allowable velocity of 1.06 feet per second found at M.L.T. for 
midsummer periods. Hence, during the major portion of the stratification 
season when cold water is at a premium, the structure is capable of supply- 
ing all of the demand from the coldest bottom waters. 


Performance and Benefits 


The skimmer wall and Swan Pond dikes were completed in October 1954 
and the underwater dam was constructed during the winter of 1955-56. Per- 
formance data are available, therefore, on the operation of the system with 
and without the underwater dam. The efficacy of the skimmer for widely 
varying degrees of stratification is proven by the two vertical temperature 
profiles shown in Figure 11. These profiles were measured immediately up- 
stream from the skimmer. The data in Figure 11 (a) show operation when 


a = 0.0011, while Figure 11 (b) shows a profile for a oY value of 0.0004. The 


arrows indicate the temperature measured in the pump forebay at the re- 
corder location marked with a triangle on Figure 1. In June the incoming 
solar radiation heats the water flowing in the Swan Pond canal and can raise 
the temperature by as much as one degree Fahrenheit. If this is taken into 
account it is apparent from the Figure 11 (a) data that only the cold bottom 
waters were passing through the skimmer wall opening. It is equally appar- 
ent from Figure 11 (b) that at very small density differences only the coldest 
bottom waters are drawn through the wall openings. 

On June 20, 1955, velocities and flow direction were also measured at a 
section across the Emory River channel at the temperature record site, 
shown in Figure 1 just down the river from the skimmer wall. The measure- 
ments are plotted in Figure 12. The data clearly show the cold Clinch River 
waters being drawn up the Emory River. The data also show the sensitivity 
of a density underflow interface. It will be noted that the interface occurs at 
a much higher elevation on the left or south bank than on the opposite bank. 
Measurements made at other stations down river showed that the channel 
curvature caused these interface variations. When it is remembered that the 
effective gravitational factor is only 0.001 of that normally encountered in 
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engineering problems it is not too difficult to appreciate these interface 
anomalies. 

The real proof of the worth of these structures is contained in Figures 13 
and 14 which show the Emory River surface, Clinch River underflow, and 
pump forebay temperatures for the stratification seasons 1955 and 1956, 
respectively. In 1955 the underwater dam was not constructed and the pump 
forebay temperatures were, as prophesied, several degrees warmer than the 
Clinch River underflow temperatures. From Figures 10 (a) and 10 (b) it was 
predicted that during August with a steady flow of 3000 cubic feet per second, 
the temperature at the intake would be about 73.5° F. Since about 1° F. of 
heating would occur in the Swan Pond canal the pump forebay temperature 
was thus predicted to be 74.59 F. Figure 13 shows that the actual measured 
temperatures ran from 73 to 75 degrees Fahrenheit. 

Figure 14 clearly demonstrates the worth of the underwater dam. Where- 
as, the pump forebay temperatures were 3 to 5 degrees Fahrenheit higher 
than the Clinch River underflow temperatures without the dam, with the dam 
the two temperatures were nearly the same. The predicted temperatures 
shown on Figure 10 prove to be slightly higher than the actual. In late June 
1956 a steady flow of about 2500 cubic feet per second existed and resulted 
in the temperature rising to 70.59 F. The predicted temperature for this 
condition, allowing for Swan Pond heating, was about 72° F. During August, 
when the Norris flows averaged about 5500 cubic feet per second, the forebay 
temperatures ran about 65.5° F. The predicted value was 69° F. 

The temperature benefits derived from the two major structures during 
1956 are graphically shown in Figure 15. 

Prior to operation, and as a part of the economic justification of the works, 
the curves of Figure 10 were used in conjunction with an anticipated Norris 
release schedule and plant operating cost characteristics to estimate the 
financial benefits that would result from the diversion structures. The bene- 
fits were based wholly upon the purchase cost of the additional coal (not in- 
cluding costs for handling and processing through the plant) that would have 
to be burned if the plant produced nine-tenths rated capability while using 
condensing water at the temperature available without the structures. The 
Skimmer wall and the dike system, without the dam, were estimated to save 
$75,000 per year. It was estimated that the underwater dam would add at 
least $24,000 per year to these savings. In actual operation, the measure- 
ments of temperature reduction show that the savings realized from the 
skimmer wall alone in 1955 amounted to $105,000. In 1956, with both skim- 
mer wall and underwater dam effective, the savings amounted to $155,000. 

If the costs for handling and processing coal through the plant were included— 
which it is reasonable to do—then the savings realized would be increased by 

about 20 percent. Furthermore, if the costs were evaluated upon the basis of 
additional energy sales, the benefits would be still greater. 

The construction costs of the skimmer wall were $251,000, of the dike 
system $191,000, and of the underwater dam, $31,000. These total $473,000 
for all the works. Assuming that the year 1956 is representative of the 
annual benefits from the system, then in the short course of only three years 
or less, the total benefits will equal the entire cost. Thus, the skimmer and 
underwater dam are very profitable additions to the steam plant facilities. 
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SYNOPSIS 


Two layers of liquid differing slightly in specific gravity are contained in a 
channel or reservoir upstream of a submerged sluice gate. The liquid down- 
stream of the gate is homogeneous and has a specific gravity the same as the 
lower, heavier layer in the upstream portion. The denser liquid is caused to 
flow through the gate opening under the action of a slight head differential. 
The problem is to determine the maximum discharge of the lower layer 
through the gate without steady state withdrawal from the upper layer. Experi- 


mental and analytical results are presented to define this condition of incipient 
drawdown. 


INTRODUCTION 


Engineering interest in density currents has increased rapidly since the 
early days of Lake Mead when observations were made of the intermittent 
passage of silt-laden underflows through the reservoir. At the present time— 
over twenty years layer—analytical and experimental information on the dy- 

namics and behavior of density currents is still relatively meager. Density 

currents or gravity underflows occur as a result of slight density differences 
which may be due to temperature differentials, dissolved salts or suspended 
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sediment. The resulting stratified flow may contain one or more interfaces 
or zones in which large density gradients are found. 

The control of density currents has long been considered economically de- 
sirable; as an example, various barriers have been proposed for the control 
of salt water intrusion in rivers, canals and tidal estuaries. At the present 
time, the profession is beginning to realize the possibility of large scale con- 
trol of density current phenomena of even more diverse nature. The purpose 
of this paper is to report on some experimental and analytical work dealing 
with gravity underflow control by means of a submerged sluice. 

The problem which resulted in this study arose during the design of the 
condenser water intake for the Kingston steam power plant of the Tennessee 
Valley Authority. During the summer months, a thermal stratification occurs 
in the reservoir from which the condenser water is drawn. Since the thermal 
efficiency of a steam plant increases as the condenser water temperature de- 
creases, there was economic justification for a control structure which would 
result in the supply of water having the lowest possible mean temperature. A 
submerged sluice or skimmer wall was proposed as a means of accomplishing 
this objective. -1 study of the effects of density currents in the T.V.A. system 
and a discussion of the Kingston site in relation to the thermal stratification 
has been given by Fry, Churchill and Elder.(1) A companion paper being 
published concurrently by Elder and Dougherty entitled “Thermal Density 
Underflow Diversion Works for Kingston Steam Plant” discusses the Kingston 
skimmer wall design in detail. 

Inasmuch as no attempt was made in this study to reproduce the geometric 
or geographical conditions at the Kingston intake, the reader is referred to 
the above papers for further details. 

The experimental and analytical studies described in this paper may be 
considered as a fundamental investigation of the submerged sluice gate con- 
trol of a simple two-layer system having a definable interface or density dis- 
continuity. In particular, it is desired to establish the conditions for maximum 
discharge from the lower-heavier layer without withdrawal from the lighter 


layer above. Fig. 1 is a schematic diagram showing the notation used in this 
investigation. 
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Experimental Equipment and Procedure 


The experimental equipment consisted of a parallel glass wall flume 4.5 
inches wide, 18 inches deep and 8 feet long connected to an entrance tank or 
reservoir 18 inches square and 36 inches deep. Between the flume and reser- 
voir plastic cylindrical quadrants in both the horizontal and vertical planes 
form a smooth transition section. A schematic diagram of the flume is shown 
in Fig. 2. The adjustable sluice gate can be placed in any longitudinal position 
along the flume. Throughout the test program, fresh water was used as the 
upper layer liquid and sodium chloride solutions of varying specific gravity 
were used to simulate the colder, lower layer. The addition of vegetable dyes 
to the two liquids made the interface between the two layers distinct. 

In preparing the equipment for an experimental run the flume and reser- 
voir were initially filled with fresh water to a depth of about six inches above 
the flume floor. The salt solution was then transferred by means of a small 
filling pump from the mixing tank to the bottom of the reservoir. By introduc- 
ing the salt solution at a very slow rate, it could be made to lift the fresh 
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water without appreciable mixing. The salt solution entered the reservoir 
through a small diffuser and baffle which served to reduce the entrance ve- 
locities and to maintain a well defined interface between the two layers. This 
step of the process (which required from three to four hours) was continued 
until the horizontal salt-fresh water interface was a distance h, above the 
bottom of the flume. At this point, the original free surface of the fresh water 
had been raised to a distance y above the flume bottom as shown in Fig. 1. 
The sluice gate remained closed during the filling stage and salt water was 
poured into the flume downstream of the gate until it was at a level only slight- 
ly below the free surface of the fresh water on the upstream side of the gate. 
The small difference in surface levels was necessary to avoid a reverse flow 
through the gate upon opening. After the required amount of salt water was in 
the flume, the filling pump was turned off and the larger circulating pump was 
started. Appropriate changes in the valves of the system were made to isolate 
the mixing tank and to open the circuit connecting the suction side of the main 
pump to the downstream end of the flume and the discharge side to the bottom 
of the reservoir. The sluice gate was slowly raised a distance b which was 
always less than the elevation of the horizontal interface h,. The control 
valve was opened gradually and the denser liquid began to circulate smoothly 
in a clockwise direction. As long as no water from the upper layer was drawn 
through the gate the density of the circulating lower liquid remained constant. 
The total discharge through the gate was obtained by means of a small Venturi 
meter in the pump discharge line. The discharge was increased slowly by 
small incremental valve adjustments. After each change several minutes 
were allowed to elapse while the interface adjusted to the new steady state 
condition. This procedure was continued until the point of incipient drawdown 
was reached. Incipient drawdown is defined as the condition at which a slight 
increase in flow results in a depression of the interface below the gate lip 
with a corresponding steady state withdrawal of liquid from the upper layer. 
For a given geometry, density difference and interface elevation, incipient 
drawdown therefore represents the maximum rate at which liquid can be dis- 
charged from the lower layer alone. 

At the time of incipient drawdown the discharge was recorded and samples 
and temperatures of both layers were taken so that the density difference 
could be obtained for this condition. Specific gravities were determined on a 
Westphal balance reading to four decimal places. In addition, the instantane- 
ous reservoir interface elevation, h,, and the profile of the interface at draw- 
down was obtained by point gage measurement at stations along the flume. 
During the course of a test, it was noted that the reservoir interface elevation 
increased slightly and the density of the lower fluid decreased by a small 
amount; hence, the necessity for noting these readings at the time of incipient 
drawdown. This density change was due to the increased volume of the denser 
fluid caused by transient discharges from the upper layer as the drawdown 
condition was approached. It was found after some experience with the flume 
that the gradual increase in the reservoir interface elevation proved to be a 
useful means of increasing the amount of data which could be obtained from a 
single filling of the equipment. In other words, for a constant gate opening, 
the condition of incipient drawdown could be obtained for several reservoir 
interface elevations during a run. Fig. 3 shows a typical plot of corresponding 
values of discharge and specific gravity differences as a function of the time 
required for a single test. At the conclusion of this run, the difference in 


specific gravities of the two layers was 78 per cent of the value at the begin- 
ning of the test. 
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Experimental Results 


A summary of the experimental runs and data for incipient drawdown con- 
ditions is given in Table I. The fourteen runs in series A were made with the 
sluice gate near the entrance to the flume (3" downstream from the point 
where the parallel wails begin) and the remaining fifteen runs in series B 
were made with the sluice gate 82 inches downstream from the entrance. Gate 
openings ranged from 1.0 to 2.5 inches and density differences between the 
two layers varied from a maximum of six-tenths of one per cent (.0060) to a 
minimum of one-quarter of one per cent (.0025). These values correspond to 
density variations caused by temperature differences between 20 and 40 de- 
grees (F) for fresh water. In addition, Table I gives the instantaneous value 
of the reservoir interface elevation and the discharge at the time of incipient 
drawdown. A dimensionless grouping of the pertinent variables (which ap- 
pears in the theoretical studies) is given in the last column. 

The following is a brief description of the observed flow conditions in the 
flume for a typical run prior to the time at which incipient drawdown took 
place. 

At the start of a run, the density interface was horizontal and coincident 
with the interface elevation in the reservoir end of the flume. After the es- 
tablishment of a small steady state flow through the gate a longitudinal line 
defining a zone of high velocity gradients immediately became evident. This 
line, which is called the separation streamline, is shown schematically in Fig. 
4. It was noted that the line of high density gradients (interface) remained 
essentially horizontal. For convenience, the wedge shaped region delineated 
by these two lines is called the separation wedge. The outline of the wedge 
was clearly observed by introducing dye below the density interface upstream 
of the gate. Any coloring injected below the separation streamline was rapid- 
ly carried away by the flow while that above remained and served to make the 
wedge visible. Due to the shear stresses at the separation streamline a low 
velocity circulation pattern was developed in the wedge region. As the flow 
rate was successively increased the curvature of the separation streamline in 
the vicinity of the gate increased and the length of the separation wedge gradu- 
ally decreased due to entrainment with the moving lower layer. The disap- 
pearance of the wedge was also hastened by transients accompanying the suc- 
cessive changes in discharge. At the time of incipient drawdown, the density 
interface coincided with the separation streamline. Fig. 5 shows a series of 
interface profiles at the drawdown condition for runs 63, 65-68. In this series, 
the density difference and gate opening were constant while the reservoir 
interface elevation increased. At flow rates beyond the point of incipient 
drawdown, the density interface passes through the gate below the lip and is 
maintained for some distance downstream. The two-phase flows were not 
investigated quantitatively due to the difficulty of maintaining steady state 
conditions with the existing equipment. 


Theoretical Considerations 


A completely general analytical approach to the submerged sluice gate dis- 
charge of a stratified fluid is complicated by the presence of the separation 
wedge previously noted. Hence, the theoretical considerations will be con- 
cerned with flow conditions at incipient drawdown when the separation wedge 
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TABLE I 


SUMMARY OF RUNS AT _INCIPIENT DRAWDOWN 
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has disappeared and the interface has the form shown in Fig. 1. In the analo- 
gous single-fluid submerged-sluice problem, the normal analytical approach 
is to write the total head equation between a section of parallel flow upstream 
of the gate and the section through the vena-contracta below the gate. There- 
fore, the surface elevation difference across the sluice is introduced as a 
primary variable in the discharge equation. This procedure is not advan- 
tageous in the two layer system under consideration, since (as will be shown 
later) the presence of the upper liquid greatly reduces the gravitational driv- 
ing force required to produce flow and results in surface elevation differences 
across the gate which are too small to measure with any degree of accuracy. 
Thus, it is desirable to avoid introducing this term which, aside from the 
measurement difficulties, is not particularly useful in practice. The down- 
stream section for evaluation of the total head is therefore chosen to coincide 
with the plane of the sluice gate. It is recognized that in choosing this section 
the curvilinear nature of the flow will have to be considered in the total head 
equation. 

Referring to the definition diagram, Fig. 1, the total head equation is written 
between the reservoir and the vertical section at the gate using the flume 
bottom as a datum. It is assumed that the stratification originates in an area 
of such extent that the flow velocities in the reservoir region may be neglect- 
ed by comparison with those in the vicinity of the gate. If AY represents the 
difference in specific weights of the two fluids, then Y and ¥ - AY are the spe- 
cific weights of the lower and upper layers respectively. The total depth of 
the two layers upstream of the gate is y and h, equals the elevation of the 
density interface in the reservoir. The mean velocity Vj at the gate is given 
by the discharge per unit width dividided by the gate opening b. Using the 
above notation, the total head equation becomes: 


(= Ab) (¥ — he) (b-aty—d) 
+ hy 


2 
+ Ahm) + SM (1) 


The left side of the equation represents the piezometric head and also the total 
head in the reservoir in accordance with the assumption of negligible reser- 
voir velocities. The first two quantities on the right side represent the mean 
piezometric head in the plane of the gate. Due to the curvature and inclination 
of the streamlines in this region the piezometric head varies between the gate 
lip and the flume bottom. The length Ahp is therefore introduced to express 
the amount by which the mean piezometric head differs from the value b which 
would be obtained for hydrostatic pressure distribution. The kinetic energy 
head at the gate is represented by the third quantity on the right side of Eq. (1) 
where a@ is the kinetic energy correction factor required by the non-uniform 
velocity at the gate section and the subscript i refers to the critical condition 
of incipient drawdown. Finally, hg is the total resistance loss in the approach 
channel between the reservoir and the sluice gate. The total loss is a function 
of the boundary shear on the bottom and side walls of the channel and the 
frictional resistance along the interface, 

Expanding the terms containing (7 - AY), Eq. (1) reduces to: 


2 
ab ab &(Vp) 
= th, + = + he (1a) 


which can be put in the following dimensionless form after dividing each term 


by 
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It should be noted that Eq. (2) represents the interrelationship of the basic 
variables h;, b, A¥/Y and Vp only at the critical point of incipient drawdown. 
Prior to the drawdown condition, the total head equation should contain a term 
accounting for the additional pressure at the gate due to the separation wedge. 
Beyond the incipient point the equation should again be changed to account for 
the two-phase flow through the sluice gate. Eq. (2) is therefore useful in 
correlating the experimental data which was obtained at precisely the flow 
condition for which the equation is applicable. The term containing Vp is 
readily seen to be a generalized Froude parameter in which the acceleration 
of gravity is reduced by AY/?Y due to the bouyant effect of the upper layer of 
fluid. Defining the reduced gravitational constant as gt = g(A¥/Y), then 


Vp): 
(3) 
g'b 
and Eq. (2) becomes: 
hy Ah 
b ob, 2 Ab, (4) 


Eq. (4) can be shown to reduce to the elementary concept of free surface flow 
on a broad crested weir in which the critical depth (b) is assumed to exist 
over the entire length. For the case in which the upper fluid is air AY = 7, 
neglecting curvature Ah,, = 0, for uniform velocity distribution @ = 1 and neg- 
lecting friction hg = 0; then hy/b = 3/2 since Fj = 1. 

In order to simplify the correlation of the experimental data with the above 
analysis, emphasis will be placed primarily upon the situation in which the 
sluice gate is located near the reservoir entrance. This is equivalent to the 
assumption that the head loss term in Eq. (4) is negligible. Approximately 
one half of the experimental runs (section A of Table I) were conducted with 
the gate three inches from the flume entrance and may be assumed to satisfy 
this condition. 

Various aspects of the problem of pressure distribution in curvilinear free 
surface flow have been treated by Rouse, (2) Jaeger(3) and more recently by 
Khafagi and Hammad.(4) The latter have developed expressions for pressure 
distributions appropriate to the type of curvilinear flow obtained near the 
critical sections of free surface Venturi flumes. The influence surface incli- 
nation (streamline convergence) and surface curvature were analyzed inde- 


pendently with the combined mean departure from hydrostatic conditions being 
given by an equation of the form: 


2 
Ahm = + #,(0)) (5) 


where R is the radius of curvature and @ the angle of inclination of the inter- 


face (or surface) as shown in Fig. 6. The Khafagi-Hammad functions for R 
and @ are given by: 


3 2 b 


where k 2 


1+ 2e" 2.12R (6) 


#(R) = ( 
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f,(@) where 5 tan( + 4 (7) 
Substituting Eq. (5) into the total head Eq. (4) and noting that hy = 0, the follow- 


ing is obtained if R and 6 are evaluated at the drawdown condition: 


2 2 
hy F 
=) (4,(R) + + (8) 
Since it is known that for given boundary conditions the surface configuration 
of free surface flow is a function of the Froude parameter it may be assumed 
that the geometric functions given by Eqs. (6) and (7) are therefore dependent 
upon Fj. For convenience, a new parameter 7 is defined such that, 


2 
F 
= (f,(R) + t,(@)) (9) 


When R =~ and @ = 0, the pressure distribution must be hydrostatic and it 
can readily be shown that for this condition 7 reaches a maximum value of 
unity. According to the foregoing reasoning, the dependence of 7 on Fj can be 
investigated if R and 6 can be obtained experimentally. In two of the experi- 
mental runs at the drawdown condition photographs of the interface at the gate 
were obtained which were sufficiently clear for evaluation of R and @. The 
experimental data and the plotted values of 7) and Fj are shown in F“g. 7. 
Some indication of generality and an extension of the relation of higher values 
of F, can be obtained by making use of some experimental results of Khafagi 
and Hammad for free surface flows having combined inclination and curvature 
in a Venturi flume. These values are also shown in Fig. 7 and compare favor- 
ably with the stratified flow results. 

Combining Eqs. (8) and (9) the basic total head equation is given by: 


2 
+ 

b 2 (10) 
where the dependence of 7 on Fj is given by Fig. 7. 

Eq. (10) can be plotted using hy/b and F; as independent variables and a 
family of curves depending on the kinetic energy correction @ will be obtained 
as shown in Fig. 8. The corresponding experimental data from part A of 
Table I are also shown for comparison with the analytical results. In general, 
the experimental data follow the trend indicated by Eq. (10). For the low 
values of Fj the flows are laminar and values of a between 1.75 and 2.0 are 
not unusual for a rectangular channel having a width to depth ratio of approxi- 
mately 1 to 2. For increasing flow rates a comparison of the experimental 
data with the curves of Eq. (10) indicates a progressive decrease in a toa 
minimum value of approximately 1.10. The decrease in the kinetic energy 
correction factor is attributed to increasing curvature and acceleration rates 
and the onset of turbulent flow. No attempt was made to measure o@ due to the 
three-dimensional character of the flow and the low velocities involved. 

The experimental data for the sluice gate 82 inches from the flume 
entrance (part B of Table I) are also plotted in Fig. 8 for comparison with the 
above. In accordance with Eq. (4) the theoretical curves for this condition 
should be displaced upward due to the friction loss in the flume. For identical 
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values of Fj the vertical displacement between points obtained at the flume 


entrance and at station 82 should equal the dimensionless term i 2 . No pre- 
cise methods for determining the energy loss in non-uniform subsurface flows 


have as yet been developed. Only the resistance problem for laminar uniform 
underflows has been treated in any detail. 5 


Illustrative Example 


As an illustration of the use of Fig. 8 in determining the drawdown con- 
dition, the computations for the Kingston skimmer wall will be presented. 
TVA had determined that in early spring and late fall A’/Y values as low as 
0.0005 exist while in mid-summer the density differences may produce values 
of 0.0016. The proposed TVA design was for a skimmer wall opening b = 15 
feet. Their measurements and calculations indicated the elevation of the inter- 
face upstream of the skimmer wall (hy) to be about 24 to 25 feet (zone of maxi- 
mum rate of change of density with depth). For the relatively high Reynolds 
numbers of the prototype flow the kinetic energy factor might reasonably be 
assumed at a = 1.10. From Fig. 8, for h,/b = 1.63, the corresponding value 
of F; = 1.20. Hence, in the spring and fall with g' = 0.0005 (32.2) = 0.016 ft./ 
sec.’, the critical drawdown velocity (Vp); = 0.59 ft./sec. is obtained directly 
from Eq. (3). For the mid-summer condition when g' = 0.0016 (32.2) = 0.052 
ft./sec.* the critical drawdown condition will occur at a velocity of 1.06 ft. / 
sec. The maximum discharge per foot of width from the colder underlying 
layers would be therefore qj = (Vp)i (b) = 8.8 cfs/ft. and 15.9 cfs/ft., re- 
spectively. If a discharge in excess of the critical values given above is ob- 
tained, steady state withdrawal of the warmer overlying fluid should begin. 


Summary and Conclusions 


The experimental and analytical work described herein indicates that the 
control of stratified flows involving very small density differences is both 
possible and practical. The primary aim in the presentation of the experi- 
mental data for the small scale laboratory work was to achieve a correlation 
which would permit application of the information to the design of engineering 
structures. The family of curves in Fig. 8 delineate a zone of incipient draw- 
down which depends upon the kinetic energy correction factor. The region to 
the right of the curves represents a combination of the parameters h;/b and 
Fj which result in simultaneous draft of fluid from both the upper and lower 
layers. In the region to the left of the curves, flow from the lower layer only 
will occur at less than the maximum rate. 

In engineering applications it can reasonably be assumed that high values 
of @ will never be obtained. The uncertainty in the choice of the proper kinetic 
energy factor, however, is by no means as important as that involved in the 
assumption of a discrete density discontinuity for the determination of h,. In 
practice the assumed interface might be logically chosen to coincide with the 
elevation at which the rate of change of density with depth is a maximum, 
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MECHANICAL ANALOGS AID GRAPHICAL FLOOD: ROUTING* 
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The dampening effect of valley storage on a wave moving through a river 
reach is a familiar phenomenon, and hydrologic literature abounds with 
descriptions of routing methods—analytical, graphical, electric analogs, and 
mechanical devices.(1,2,3,4,5) A similar dampening is apparent if one com- 
pares the rear-wheel track of a moving vehicle with that made by a front 
wheel, although closer examination will reveal that the two effects are not 
alike in all respects for a vehicle with fixed wheel base. Nonetheless, study 
of the two phenomena reveals some rather interesting results, pointing the 
way for development of mechanical analogs and simplified graphical solutions 
to flood routing problems. 


Linear Storage-Discharge Function 


Let us first consider the case of a linear storage-discharge relation, i.e. 


S =KQ (1) 


where S is the reach storage, Q is the reach outflow, and K is a constant hav- 
ing the units of time. Figure 1 shows schematically the requirements for an 
analog which will satisfy Equation 1. In this diagram, the time scale for out- 
flow is shifted K time units to the left of that for the inflow, and the tongue 
connecting the inflow and outflow pens is adjustable so as to maintain the time 
displacement between the two at K units. This is the underlying principle of 
the patented Harkness flood router.(4) It will be seen that the slope of the out- 
flow hydrograph dQ/dt is equal to the difference between inflow and outflow 
divided by the storage factor K. Since I - Q is the time rate of change in 
storage 
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(2) 


which is equivalent to Equation 1. 


A simple, graphical solution for Equation 1 is also apparent from Figure 1. 
If the outflow up to any initial time is plotted K units to the left of the inflow 
hydrograph to be routed, a straightedge connecting simultaneous values of I 
and Q has the slope dQ/at, and the outflow can be extended for a reasonable 
time with this slope.(3) Moving to successive later times, the inflow hydro- 
graph is routed segment by segment. The routing period (i. e., the length 
with which the outflow curve is extended in one step) can be varied at random, . 
depending upon the rate of change of dQ/dt. Moreover, the outflow and inflow 
time scales can be made to coincide by placing the straightedge so as to con- 
nect the outflow with simultaneous inflow lagged K time units (see Figure 2). 


In this case, it is preferable to have a straightedge calibrated in time units 
to assist in lagging the inflow. 


Muskingum Storage Function 


The Muskingum method of routing developed by McCarthy(1) is based on 


the assumption that storage is proportional to a weighted value of inflow and 
outflow, i.e. 


S=K[xI +(I-x)Q] (3) 


where x is visualized as a weighting factor which takes into account wedge 
storage. Figure 3 illustrates schematically the requirements for a mechani- 
cal analog satisfying the Muskingum storage function. In this case the tongue 
is jointed at some point, the location of which is determined by the relative 
weights assigned to I and Q in the storage function. In addition to tracing the 
inflow hydrograph, the inflow portion of the tongue must always be pivoted so 
as to have the slope dI/dt. Then 


1-Q = K(1-x)$@ (4) 


which is equivalent to the Muskingum storage function. 
Graphical Muskingum routing can also be accomplished with the aid of a " 

jointed straightedge of the general character shown in Figure 4. Actually, 

experience has shown that the Muskingum equation provides no more flexi- 

bility than does the linear storage-outflow relation employed in connection 

with a lagged inflow. Since graphical solution of the simpler equation can 3 


readily incorporate fixed lag L, there is little to recommend use of the 
Muskingum function. 


Variable Lag and Storage Factor 


There are numerous reaches for which neither the Muskingum nor the 
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lag-and-route (constant K and L) techniques provide satisfactory results. In 
fact, it is quite generally found that the results of these approaches can be 
improved upon if K and L are permitted to vary with flow in the reach. It 

can be shown that inflow and outlfow volumes do not check when K is assumed 
to be a function of inflow and, therefore, K must be related to outflow. Lag 
can be made a function of either inflow or outflow. That is, the inflow hydro- 
graph can be lagged as a function of the inflow discharge and the result routed 
with a variable K, or the inflow can be routed with K as a function of outflow 
and the resulting hydrograph lagged. While reasonably comparable results 
can be obtained in either manner, practical considerations, both in develop- 
ment and application, favor the assumption that lag is a function of inflow. 

It is not the purpose of this paper to describe in detail the actual develop- 
ment of routing procedures, but a brief outline for the variable K and L tech- 
nique will aid in following the subsequent discussion. If storage is a function 
of outflow only, then it follows that the peaks and nteeoe of the routed outflow 
fall on the inflow hydrograph (i.e., when e. 2065 ° = OandI=Q). Therefore, 
a good approximation to the lag an. can be derived by analysis of a series 
of inflow and outflow hydrographs (adjusted for local inflow) as shown in 
Figure 5. The time difference between the outflow peak and the correspond- 
ing value of inflow is the lag to be used for that inflow. The mean line 
through a plotting of such points constitutes a first approximation of the re- 
quired lag function. 

Replotting the inflow hydrographs lagged in accordance with the relation 
derived as above, the values of K can be obtained as shown in Figure 6. Plac- 


ing a straightedge tangent to the outflow curve (at any point except where 
dQ 


—» O),K is read as the time difference between the lagged inflow (at 


time of selected outflow) and the straightedge (Equation 2). Having determined 
K at each of a number of points for selected rises, the relation K = f{(Q) can be 
derived by plotting K vs Q. 

Relations for K and L “developed in the manner described are usually high- 
ly satisfactory. Nevertheless, comparison of routed and observed hydro- 
graphs may indicate improved results through minor adjustments. The analy- 
sis for determination of K and L may show there is little to be gained by per- 
mitting these factors to vary with flow, and in this case the simpler solution 
is to be preferred. Limited experience has shown that K = f(Q) and L = f(I) 
are quite similar for a particular reach—if L or K varies appreciably with 
flow, then both vary ina similar fashion. - 

Graphical routing with variable K and L, although more complex than when 
using constant factors, is more rapid than | many less flexible methods. It is 
greatly facilitated by the aid of a transparent template of the type shown in 
Figure 7. The left-hand edge of the cutout section represents L = f(I), and 
the right-hand edge represents K = f(Q), both as measured from the vertical 
base line. The template is moved horizontally along a straightedge to deter- 
mine a segment of the outflow hydrograph centered at the base line. Use of 
the template is illustrated in Figure 8. Intersection of the inflow hydrograph 
with the left-hand edge of the template represents lagged inflow at time indi- 
cated by the base line. Eye extrapolation of the outflow curve to the base line 
determines the value of K on the right-hand edge of the template. The point 
“A” thus obtained is connected by straightedge to the previously-determined 
segment of the outflow hydrograph and a reasonable extension to the curve is 
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constructed. This process is repeated until the routing operation is com- 
pleted. 

The mechanical analog for accomplishing this form of routing must 
naturally be more complex than that for the simple, linear function. The 
time difference between the inflow tracer and the outflow pen must be varied 
in accordance with L = f(I), while the time difference between the outflow pen 
and the anchor point for the outflow tongue must follow K = f{(Q). Modifica- 
tion for variable lag is reasonably simple, but variations in K must be geared 
to the outflow pen, and this presents quite a different problem. The required 
mechanism can be readily designed, but vertical movement of the outflow pen 
would be forced to “work” the mechanism, which would tend to induce side 
thrust on the wheels and thus possibly distort the outflow hydrograph. 

Although a more difficult task, it was believed an operator could soon ac- 
quire skill to follow the inflow hydrograph and simultaneously make mechani- 
cal adjustment for variations in K. Figure 9 is a photograph of an analog 
which has been designed and constructed accordingly. The left-hand knob is 
used for maintaining the crane assembly at the intersection of the inflow and 
lag curves while the carriage is moved to the right along the time axis. The 
right-hand know is used to vary K in accordance with the routed outflow. 
Preliminary operational tests indicate this pilot model to be quite satisfac- 
tory, but it does have certain limitations which will require minor modifica- 
tions for improved efficiency and convenience. 


SUMMARY AND CONCLUSIONS 


Graphical techniques for routing directly on the plotted hydrograph charts 
have been discussed for (1) the case of simple, linear storage-discharge 
function, (2) the case of the Muskingum equation, (3) the case of linear stor- 
age function in conjunction with constant lag (translation), and (4) the case of 
both lag and the storage factor varying as complex functions of flow in the 
reach. The fourth category represents what is believed to be one of the most 
flexible methods of routing yet devised. Amazingly reliable results have been 
obtained for reaches which are incompatible with storage functions of the 
Muskingum and other recognized methods of routing. The time required for 
procedure development compares favorably with any other known technique. 
Moreover, the graphical solutions illustrated are operationally rapid and re- 
liable. The basic features of mechanical analogs comparable to the graphical 
solutions described have also been discussed. An analog which will facilitate 
variable K and lag has now been constructed and is undergoing tests. 
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TIDAL MOVEMENT IN THE CAPE COD CANAL, MASSACHUSETTS?# 


B. W. Wilcox! 
(Proc. Paper 1586) 


SYNOPSIS 


Variations in the form of the tide wave at selected points, analyses of ob- 
servations to obtain harmonic constants, and a practical method for predict- 
ing the tidal currents in the canal are discussed. Diagrams are included to 
show the shapes of the tide curves and also the agreement between predic- 
tions and observations of the tidal current at the railroad bridge. 


INTRODUCTION 


For many years the U. S. Corps of Engineers has been engaged in the im- 
provement of the Cape Cod Canal. The land cut has been widened and deep- 
ened. Dikes have been built and new channels dredged in Buzzards Bay. 
Fig. 1 shows the canal and the greater part of Buzzards Bay. The improve- 
ments have changed the tide and the tidal current in the canal necessitating 
new observational data to establish the pattern for predictions which are of 
considerable importance to the navigator. 

This discussion of tidal movement divides naturally into two parts: the 
vertical rise and fall of the water which is ordinarily called the tide; and the 
horizontal flow or tidal current. 

Automatic tide gages furnished by the Coast and Geodetic Survey have 
been operated by the canal engineers at a number of selected places during 
the years 1955-1957. During this same period tidal current observations at 
both ends of the land cut were made by the Coast and Geodetic Survey. 


Note: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1586 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 


a. Presented at meeting of the ASCE, New York, N. Y., October, 1957. 
1. Chief, Predictions Branch, Tides and Currents Div., Coast & Geodetic 
Survey, U. S. Dept. of Commerce, Washington, D. C. 
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Tides 


The difference in range and phase of the tide in Buzzards Bay and Cape 
Cod Bay produces a periodic change in the slope of the water in the canal 
which in turn causes the periodic movement of water in the canal. Slope 
lines were drawn for 12 lunar hours between stations 625,379,200 and 35. 
These slope lines determined the high and low water profiles of figure 2. 
Observations at station 452 and 320 fit well into the profiles thus furnishing 
a check on their validity. The high and low water profiles in figure 2 show 
the slight decrease in mean range from 3.9 feet at Abiels Ledge, Buzzards 
Bay to 3.5 feet at the entrance to the land cut (RR. Bridge) and then the sharp 
rather uniform increase to 8.7 feet at the mooring basin (Sandwich) near the 
Cape Cod Bay entrance to the canal. 

The great difference between the tides at the two ends of the canal is 
brought out clearly in figure 3. The tide at Sandwich (station 35) is of the 
regular semidaily oceanic type with the duration of rise equal to the duration 
of fall. The distortion appearing in the Buzzards Bay wave is due to what is 
commonly called “shallow water effects.” The duration of rise is not the 
same as the duration of fall. Note the change in the low water part of the 
tide wave as it progresses from station 625 to station 379. At Abiels Ledge 
(station 625) the rise of the tide takes longer than the fall whereas at the 
railroad bridge (station 379) the duration of rise is less than the duration of 
fall. Before the change in the channel and the building of the dikes the shape 
of the curve for station 379 was similar to that for station 625. 

The tide changes continually as the tide producing forces vary with the 
motions of the moon, sun and earth. Inasmuch as the equations of motion of 
moon, sun and earth are known it is possible to compute the tide producing 
forces for any desired time and hence the resulting tide. Present theory is 
not sufficient for mathematical prediction of the tide without observational 
data. Observations are necessary because each body of water reacts to the 
tide producing forces in its own way. 

The well known harmonic equation(1) of the periodic tide is 
y = Zo +L fH cos (at + Vo + u - K) 


in which 
y = height of tide at time t 
Zo = height of: mean sea level above datum of predictions 
H = mean amplitude of any constituent A 
f = factor for reducing mean amplitude H to year of prediction 
a = speed in degrees per solar hour of constituent A 
t = time in solar hours reckoned from some initial epoch 
Vo+u_ = theoretical phase (equilibrium argument) of constituent A 
when t =0 
K = difference between theoretical and actual phase of constituent A 


as determined from observations 


The number of constituents involved will vary with the complexity of the 
tide. 
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The mean amplitude (H) and the epoch (x) for each constituent are the 
harmonic constants to be obtained from observed hourly heights of the tide. 
The problem is simplified by the fact that the speeds of the constituents 
which may be involved are known in advance. There are a number of ways 
to make the analysis. The one used by the Coast Survey involves the assign- 
ment of each observed hourly height to a constituent hour. There are 24 con- 
stituent hours in a constituent day. Twenty four sums (one for each constitu- 
ent hour) are obtained from the entire series. Modern electronic computing 
machines are now used in the summation of long series. The 24 hourly sums 
for constituent A contain the augmented amplitude of A but the constituents 
of other periods are largely cancelled. Short period meteorological effects 
which are contained in the observations are filtered out in the process of 
analysis. From the 24 hourly sums the amplitude H and the epoch x of A 
are readily obtained. Some refinements are required to eliminate inter- 
ference of constituents which the process of summing can not complete. 

The harmonic constants thus obtained may be used in making predictions. 
This is an arduous task unless special machines are used. Many of the 
principal maritime nations have had analog computers built especially for 
tidal predictions. By the process briefly outlined above the Coast and 
Geodetic Survey has analyzed the tidal observations made in 1955 at stations 
379 and 35. The principal tidal harmonic constants are listed in table 1. 


Table 1 
Tidal Harmonic Constants(1) 


M2 S2 N2Q Kj Oj Pj M4 Mg 


Sta. 35 H ft. 4.08 0.62 091 044 0.35 0.14 0.12 0.14 
K deg. 326 002 286 135 111 135 106 243 


Sta. 379 H ft. 162 0.28 0.46 0.27 420.18 0.09 0.29 0.07 
K deg. 268 276 244 113 128 113 139 236 


The tidal harmonic constants for these two stations adjusted by current ob- 
servations are used in predicting the tidal current at station 379 - the method 
will be discussed later under currents. 


Currents 


Strong tidal currents flow back and forth periodically through the land cut, 
the westbound current being the stronger. A velocity of 5.3 knots has been 
observed in midchannel at the railroad bridge. In entering Cape Cod Canal 
land cut from Buzzards Bay it is natural to call the current from Buzzards 
Bay into the land cut the flood current; but the canal may also be entered 
from Cape Cod Bay and it is therefore just as natural to regard the west- 
ward current as the flood current. In the canal, clearly, upstream and down- 
stream cannot be used in designating the flood and ebb current, and so the 
time relations between current and tide must be used. The observations 
show that the strength of current setting eastward at the RR. bridge occurs 
about three hours before local high water; that is, on a rising tide. Hence 
the flood current at the RR. bridge is the one setting eastward from 
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Buzzards Bay into the land cut. In fact the eastward current is the flood cur- 
rent throughout the canal for the tide and current data at the Cape Cod Bay 
entrance also show that the eastward current is on the rising tide. 

The following table gives the average time intervals referred to the 
Greenwich transit of the moon for various phases of the tidal current and 
the average maximum velocities for the entrances to the land cut as deter- 
mined from observations in 1955 and 1956. 


Table 2 


Slack Maximum Slack Maximum Velocities 
Flood Ebb Flood Ebb 
hours hours hours knots knots 


8.25 11.06 2.22 4.68 4.0 4.5 


hours 


RR. Bridge 
Sta. 379 


Mooring Basin 
Sta. 35 


8.03 10.86 


2.03 4.45 


2.4 2.7 


Observations at station 379 were made with Price meter at depths 10 feet 
and 20 feet below the surface. The currents at the two depths were so nearly 
equal that the values were meaned for table 2. At Station 35 the currents 
were measured at 15 foot depth by a Price meter. The time of the current 

is practically the same throughout the canal which verifies the theory that 

the flow is almost entirely governed by the changing surface head in the 

canal. The velocity at the RR. bridge is much greater than at the Mooring 
Basin. This is due to the tide being about two hours later at the Mooring 
Basin than at the bridge. 

Advance information about the currents in the canal is extremely important 
to the navigator. Daily predictions of the slacks and strengths of current con- 
tained in the annual tidal current tables published by the Coast Survey are 
based on observations of both tides and currents. The differences in tidal 
range and the times of high and low water at the two ends of the land cut 
cause the water surface at one end to alternately rise above and fall below 
the level at the other end, thus creating a periodic reversing current in the 
canal. Theoretically, disregarding friction or inertia, the velocity of the 
current would vary as the square root of the surface head, being zero when 
the surface is at the same level at both ends and reaching a maximum when 
the head is greatest. Actually there will generally be a lag of some minutes 
in the response of the current movement to the surface head which must be 
determined by observations. The problem of predictions therefore, reduces 
to one of predicting the surface head between the two ends. This may be done 
by means of harmonic constants derived from the tidal constants for the 
entrances. 

For any one constituent let T represent time as expressed in degrees 
reckoned from the phase zero of its Greenwich equilibrium argument. Also 
let Yq and Yb represent the height of the constituent tide for any time T as 
referred to the mean level at stations 379 and 35 respectively and let Yq = 
surface head = Yg-Yp. Tidal harmonic constants for stations 379 and 35 are 
obtained from analyses of observations as outlined under tides. H, is the 
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mean amplitude and G, is the Greenwich epoch of a constituent for station 
379; Hj} and Gp are the corresponding harmonic constants for the same 


constituent at station 35. Formulas for the heights and surface head may 
now be written 


Yq = Hg cos (T-Gg) for station 379 
Yb = Hp cos (T-Gp) for station 35 
Yq = Ha cos (T-Ga) -Hp cos (T-Gp) 


Harmonic constants (Hg and Gg) for the surface head may be obtained by tak- 
ing the vector difference between H,,G, and Hp,Gp. Hence the surface head 
for any one constituent may be written Yg = Hg cos (T-Gg). 

Separate computations are made for each of the tidal constituents. The 
values obtained for Hg and Gq are the corresponding amplitudes and 
Greenwich epochs in a harmonic expression for the continually changing dif- 
ference in elevation of the water surface at the two ends of the land cut. 

For predictions of the current the amplitude Hyg must be modified 
through the formula (Velocity)2 = constant (C) x Hg. Furthermore the epochs 
Gg need to be changed to allow for the lag in response of the current to the 
changing difference in water level at the two entrances. 

Following the above outlined method the tidal harmonic constants from 
table 1 were combined vectorially and adjusted for the observed current at 
the RR. bridge. The constants are listed in table 3. 


Table 3 


Tidal Current Harmonic Constants 
M2 No Kj Py M4 Mg 


RR. Bridge Hft. 8.73 1.65 1.62 0.52 0.47 0.17 O.51 0.18 
Sta. 379 K deg. 185 223 150 353 281 352 189 117 


With the modified harmonic constants direct readings from the tide 
predicting machine would give velocity squared. However a special square 
root dial on the machine makes it possible to read the velocities directly. 

Figure 4 shows a comparison of predicted tidal current with observed. 
The period was picked at random. The differences are small with one ex- 
ception. Part of the deviation of predictions from observations is due to non 
periodic variations in the observations and part is due to the omission of 
factors other than surface head involved in the current. But it is evident 
from Figure 4 that the other factors have a very small effect and that the 
predictions through surface head are dependable. 


RE FERENCE 
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DISTRIBUTION OF SEDIMENT IN LARGE RESERVOIRS4 


Whitney M. Borland,* A.M. ASCE 
and Carl R. Miller,** A.M. ASCE 
(Proc. Paper 1587) 


SYNOPSIS 


In planning and developing multi-purpose reservoirs it is necessary to 
make an allocation of storage space for sediment accumulation. This sedi- 
ment accumulation, in turn, influences the design of the dam because of the 
effect on active storage capacity requirements, outlet sill elevations, recrea- 
tional facilities, and backwater conditions. It is, therefore, necessary to 
predict the probable distribution of the sediment in the reservoir with time. 

Two methods of predicting the probable sediment distribution are given 
and an illustrated example shown. One method is basically mathematical 
(Area-Increment Method), whereas the other is a mathematical procedure 
developed on the basis of actual experiences (Empirical Area-Reduction 
Method). The area in which additional study is needed to evolve procedures 
that will give consistently reliable predictions of sediment distribution under 
varying conditions is pointed out. 


INTRODUCTION 


Unfortunately, when a dam is constructed to store water for useful pur- 
poses, sediment being transported by the water is also stored. Unlike water, 
the sediment cannot be effectively released through the reservoir or put to 
useful purposes. It is inevitable, in most cases, that the useful storage ca- 
pacity behind a dam will be reduced in time by sediment accumulation until 


Note: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1587 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 
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replacement storage is necessary. In many countries reservoirs construct- 
ed years ago have become completely filled with sediment deposits rendering 
them practically useless. As most large reservoirs in this country are rela- 
tively new, only a small number have yet been substantially reduced in ca- 
pacity. However, future generations will be faced with replacing storage 
capacity lost to sediment in our major reservoirs. Replacement of many 
reservoirs may not be easy since the most economical reservoir sites are 
currently being used. Removal of sediment from a reservoir by mechanical } 
means is not economically feasible at present but increases in the value of 

an acre-foot of water may change this. Various methods have been suggested 
and tried for removing, reducing, or by-passing sediments. Although none 
has yet supplied us with a solution, particularly in large reservoirs, perhaps 
some day when the need becomes acute a method will evolve. Of course, the 


economic value of a reservoir may become less with new methods of power 
production and water supply. 


The Problem 


In present day analysis we seldom look beyond one or two generations in 
planning the water supply to be provided by a structure and often the struc- 
ture size is dictated by available water supply and economical considerations. 
Indeed the pay out period on all projects built by agencies such as the Bureau 
of Reclamation is usually less than 100 years. However, in planning the 
storage allocations it is the Bureau’s policy to provide for 100 years of sedi- 
mentation unless it can be demonstrated by economic analysis that a lesser 
allotment is justified. 

There was a time when it was believed that sediment always deposited in 
the bottom elevations of a reservoir rather than depositing throughout the 
full range of reservoir depth. It is now fully realized that sediment deposits 
will spread throughout the reservoir reducing the incremental capacity at all 
elevations (Figure 1). Factors that influence the mode of deposition include: 


1. Reservoir shape 

2. Sediment characteristics 

3. Reservoir operation 

4. Sediment-reservoir volume ratio 
5. Inflow-capacity relationship 


There are other factors, of course, such as salt cedar infestation, amount 
and location of tributary inflow Within the reservoir, etc., but these, when 
existing, require special considerations. 

Of immediate concern in planning and designing a dam and reservoir, is 


the effect of sediment accumulation on design considerations. Major design 
considerations include: 


1. Elevation to which sediment will accumulate at the dam in a given 2 
period of time. This information will affect the design elevation of the river 
outlet and power penstock sills. Sometimes an adjustable stoplog type trash- 
rack arrangement is used which provides flexibility in the outlet elevation to 


allow complete use of the water storage at alltimes. Stoplogs are added as 
the sediment accumulates. 


2. Reduction of active storage capacity. If this loss can be predicted the 
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structure height can be set to provide the required water supply storage. 
Often, the effect of sediment accumulation on the water supply is evaluated 
through the operation studies. 


3. Deposits of sediment in the upper reaches of the reservoir.(1) These 
deposits can result in a significant increase in the water surface elevations 
upstream. This will, in turn, endanger upstream installations and develop- 
ments, particularly during periods of high water inflow. The dangerous con- 
dition results because the reservoir water surface is usually high when the 


inflowing stream is at flood stage. This condition must be considered in the 
designs. 


4. Effect on recreational developments around the reservoir. The sedi- 
ment deposits could influence the location of such facilities, especially when 


the reservoir is one that will have considerable drawdown exposing sediment 
deposits. 


5. Reservoir planned solely for sediment accumulation. In this case the 
probable sediment accumulation is of primary importance. A considerable 
portion of the total sediment stored may be above crest elevation of the 
reservoir. 

It is therefore, evident that estimating only the amount of sediment inflow 


to a reservoir is not sufficient. We must try to predict how and where this 
sediment will deposit. 


Methods for Predicting Sediment Distribution 


Having discussed the need for predicting the manner in which sediment 
will be distributed, two methods that can be used in making this prediction 
will now be presented and illustrated. The first method which is strictly 
mathematical is known as the Area-Increment Method, and was developed by 
E. A. Cristofano.* The second method, which is a mathematical procedure 
developed on the basis of actual occurrences in large reservoirs, is identified 
herein as the Empirical Area-Reduction Method.(2) Both procedures involve 
the adjustment of the original surface areas to reflect the decrease in area 
with sedimentation. (Figure 2). 


Area Increment Method 


The area increment procedure is based on the assumption that the sedi- 
ment deposition in a reservoir can be approximated by reducing the reservoir 
area at each reservoir elevation by a fixed amount. This involves a series of 
approximations. By using the average end area or prismoidal formulas the 
reservoir capacities are computed on the basis of reduced surface areas un- 
til the total reservoir capacity below maximum normal water surface is 
essentially the same as the predetermined capacity obtained by subtracting 
the sediment accumulation with time from the original capacity. Using 
Alamogordo Reservoir on the Pecos River in New Mexico(3) as an example, 
the step procedure for the Area-Increment Method follows. The procedure 


* Former Hydraulic Engineer, Project Office, Bureau of Reclamation, 


Albuquerque, New Mexico. 
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| 
| 


1587-4 HY 2 April, 1958 


illustrated is an adaptation developed by W. T. Moody.* The basic informa- 
tion for Alamogordo Reservoir as obtained in the 1944 resurvey is: 


Original capacity at El. 4275—156,750 acre-feet 
Capacity after sedimentation—132,170 acre-feet 
Sediment accumulation—24,580 acre-feet (1936-1944) 
Original Reservoir depth at dam—125 feet (El. 4150 to 4275) 
Elevation of sediment at dam in 1944--4190 (approximately) 


Step Procedure—Area-Increment Method (See Table 1) 


The area increment procedure can be expressed by the following basic 
equation: 


Vs 


Ao (H - ho) + Vo 


where: Aj» = area correction factor in acres which is the original reser- 
voir area at the new zero elevation at the dam, 


Vo = sediment volume below new zero elevation (acre-feet), 


Vg, = sediment volume to be distributed in the reservoir (acre- 
feet), 


H = reservoir depth at the dam—streambed to maximum normal 
water surface (in feet), 


ho = depth in feet to which reservoir is completely filled with 
sediment—new zero elevation 


This equation states mathematically that the total volume of sediment Vg, 
consists of the portion which is uniformly distributed vertically over the 
height H - ho, plus the portion below the new zero elevation of the reservoir. 
Using this equation the step procedure for the Area-Increment Method is: 


Step 1. From the equation, hp can be obtained by trial and error as fol- 
lows: 


Vg = 24,580 acre-feet 
H = 125 feet 


Vo, hg, and Ap are obtained by assuming h, and reading Vo and Ap from 
the original area-capacity curves (Figure 6) at that indicated elevation. 
First trial—assume ho = 25 feet 
then Ao = 150 acres Vo = 1,600 acre-feet 
so 24,580 = 150 (125-25) + 1,600 + 16,600 
Second trial—assume h, = 34 feet 
then Apo = 235 acres Vj, = 3,300 acre-feet 
so 24,580 = [ 235 (125-34) + 3,300] = 24,680 
The new zero elevation at dam = 4150 + 34 = 4184 
Area correction factor = 235 acres 


Step 2. Compute accumulative sediment volume by applying the area cor- 


Engineer, Technical Engineering Analysis Section, Division of Design, 
Bureau of Reclamation, Denver, Colorado. 
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correction factor at each depth increment (Col. 4) and computing volumes by 
average end-area formula (Col. 5). Computed volume should check the pre- 
determined sediment volume within 1 percent. 


Step 3. Reduce original area at each increment by the area correction 
factor of 235 acres to give revised area table, (Col. 6 = Col. 2 - Col. 4). 


Step 4. Determine the revised capacity table by reducing the original 


capacity at each increment by the sediment accumulation, (Col. 7 = Col. 3 - 
Col. 5). 


The results obtained are compared with the actual resurvey curves in 
(Figure 6). 

The basic equation for this method is very useful in making a quick esti- 
mate of the depth of sediment at a dam over a given time period. By refer- 
ring to the original area-capacity curves and knowing the amount of sediment 
that will accumulate in x years, we can, with a few approximations, balance 
the equation and thereby derive the height to which sediment will accumulate 
at the dam. Once this is established the area reduction can readily be applied 
at each interval of depth up to maximum normal water surface and the proba- 
ble area-capacity curves after sedimentation can be determined. 

It should be noted that there is no smooth reduction in area from the 
original to the revised condition between maximum normal water surface and 
the next increment above in the flood control pool. This can be corrected, if 
desired, by allowing some sediment to deposit above maximum normal water 
surface and by making the area correction factor progressively smaller 
starting a few intervals below the maximum normal water surface elevation 
and extending a few intervals above. The surface area for the before and 
after condition should be the same at the uppermost elevation. Some addi- 
tional trial and error will be necessary here to get the computed sediment 
volume to equal the predetermined volume. 

The Area-Increment Method is only an approximation of the probable loca- 
tion of the sediment deposits. It is known that sediment does not generally 
reduce the surface areas by the same amount at every elevation. However, 
comparisons have been made between the results obtained by this method and 
the actual resurvey results and good agreement has been obtained in most 
cases. It should be noted that the applicability of the method decreases with 
an increase in the sediment/capacity ratio. Also the accuracy decreases 
with any deviation of the reservoir in question from a more or less standard 
type. Usually the arbitrary stipulation is that if the 100-year sediment 
accumulation exceeds 15 percent of the original capacity a more exacting 
method of distributing the sediment should be applied. In any case, the area 


increment method is a good first check procedure for determining the prob- 
able sediment distribution. 


Empirical Area-Reduction Method 


The Empirical Area-Reduction Method for determining the probable sedi- 
ment distribution is accomplished through two main steps: 


1, Classify the reservoir using 4 basic standard type curves which were 
developed from actual resurvey data. 


2. Make a trial and error type computation using the average end area or 
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prismoidal formulas until the capacity computed equals the predetermined 
capacity. 


The resurvey data for 30 reservoirs have been used to develop 4 standard 
type curves of percent sediment deposit versus percent reservoir depth 
(Table 2 and Figure 3). The reservoirs included in establishing the 4 stand- 
ard types vary widely in location, type operation, and sediment characteris- 
tics. Capacities range from 40,000 to 30,000,000 acre-feet (with a few ex- 
ceptions). The depth is based on the difference between streambed elevation 
at the dam (zero capacity) and elevation at maximum normal water surface. 
The general classification resulting from the analysis of the data from the 
30 reservoirs is as follows: 


Reservoir Standard 
“M” type classification 
1.0-1.5 Gorge IV 
1.5-2.5 Hill Hl 
2.5-3.5 Flood plain-foothill 
3.5-4.5 Lake I 


(“M” is the reciprocal of the slope of the line obtained by 
plotting reservoir depth as ordinate against reservoir capacity 
as abscissa on log-log paper (See Figure 4). 

It should be noted that a lake type does not necessarily have 
to be out on the plains nor does a gorge type have to be in the 


mountains. The gorge type (IV) is the poorest defined of the four 
classifications.) 


In applying the standard type curves to a proposed reservoir, it is first 
necessary to determine the type that fits the reservoir involved. First 
selection of type is obtained by plotting depth versus initial reservoir ca- 
pacity on log-log paper (Figure 4). This usually gives a straight line rela- 
tionship or in some cases 2 straight lines.(4) The slope of the line (or lines) 
gives an indication of the type. Generally a curve showing rapid increase in 
capacity with depth indicates a large basin type reservoir whereas a small 
increase in capacity with depth indicates a gorge type reservoir. This in- 
formation classifies the reservoir according to shape. The standard type 
curve can then be selected and adjusted as necessary to reflect other condi- 
tions, such as anticipated reservoir operation, type sediment inflow, or 
capacity/inflow ratio. The type curve of percent depth versus percent sedi- 
ment deposit finally adopted may be a combination of 2 types or could fall in 
between 2 of the 4 basic types, but generally 1 of the 4 basic types will fit the 
conditions involved. It can be seen that considerable judgment is sometimes 
involved in this selection. 

Having selected a type sediment deposit curve the rest of the procedure is 
mathematical. The 4 standard type sediment versus depth curves have been 
converted to area design curves for use in the computations (Figure 5). The 
area design curves have relative sediment storage as the ordinate and rela- 
tive depth as the abscissa. The area under the curve in each case equals 
unity. Other curves can be produced as long as the area under the curve is 
unity. The conversion from the standard type sediment deposit curves to the 
area design curves has been made by Moody by applying the equation: 


Ap = C p™(1-p)” 
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where: Ap represents a dimensionless relative area at relative distance 
p above streambed 
C, m, and n, are dimensionless constants which 
are determined by the type of reservoir 


Numerical values of m, and n, were computed by trial and error using a 
least square procedure so as to make the analytical curve the best possible 
fit to the field survey data. With m, and n, determined C is fixed by the con- 
sideration that the total area under the *e curve must be unity. The charac- 


teristic constants C, m, and n determined for the 4 types of reservoirs are 
as follows: 


Type Sediment Storage 
near 

I 3.4170 1.5 0.2 Top 

0 2.3240 0.5 0.4 Upper middle 

I 15.882 it 2.3 Lower middle 

IV 4.2324 0.1 2.5 Bottom 


There is no logical progression of these constants as the silt factor moves 
from top to bottom. 
The procedure from this point can be handled in tabular form. Again 


Alamogordo Reservoir data are used in the illustration. The steps are as 
follows. 


Step Procedure—Empirical Area Reduction Method (See Table 3) 


Step 1. Determine the relative depth for each increment at the dam (Col. 
4). This is merely the ratio of the incremental depth to the total depth (maxi- 
mum normal water surface to streambed—El. 4275 to 4150). This is deter- 
mined from values in Col. 1. 


Step 2. Determine values of relative sediment area (A,) from the standard 
type curve selected (Figures 3 and 5). For Alamogordo a type II has been 
determined. Enter the curve with relative depth for each increment and read 
corresponding Ap values. (Col. 5). 


Step 3. Select a first approximation of the probable sediment elevation at 
the dam after sedimentation. Areas at and below this elevation will equal 
areas in Col. 2. Sediment areas for each depth increment above the estimat- 
ed new zero elevation are obtained by dividing the original area at zero ele- 
vation (Col. 2) by the corresponding Ap value (Col. 5) and multiplying this 


ratio (K) by the A, values at each succeeding increment. For example—re- 
ferring to Table 


Assume zero elevation of 4190 
Surface area = A, = 300 acres 
Ap at El. 4190 = 1.125 


300 
K = 1.125 266.7 


The new area at each succeeding elevation is the Ap at that elevation 
times 266.7. 
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Step 4. With the sediment areas established (Col. 6) the incremental sedi- 
ment volumes can be computed by the average end-area formula and entered 
in Column 7. Columns 6 and 7 represent the first approximation. If the sedi- 
ment volume summation in Column 7 exceeds or is less than that desired a 
second approximation is made. Second and third approximations are shown 
in Columns 8 through 11. Note that the new zero elevation may not occur at 
an even contour interval. In the example shown the new zero elevation is 
4183.5. When making the evaluation at odd increments the original area 
capacity curves (Figure 6) and Figure 5 are referred to. 


Step 5. Accumulate the sediment volumes shown in Col. 11 and enter in 
Col. 12. 


Step 6. Determine the revised areas for after sedimentation (Col. 13 = 
Col. 2 - Col. 10). 


Step 7. Determine the revised capacities (Col. 14 = Col. 3 - Col. 12). 


Usually 2 or 3 trials are all that are necessary to obtain the desired 
balance. If the summed volume is within 1 percent of the predetermined 
volume, sufficient accuracy has been attained. 

There are several advantages to the Empirical Area-Reduction Method. 
It permits the classification of a proposed reservoir into a standard type 
that immediately gives an indication of the manner in which the sediment will 
distribute. A portion of the sediment can readily be deposited above maxi- 
mum normal water surface (flood control space) in the computations, by 
selecting a top elevation above maximum normal water surface in determin- 
ing relative depths. (In a large reservoir the sediment deposits above maxi- 
mum normal water surface will seldom exceed 10 percent of the total ina 
50- to 100-year period, except for those dams which are for sediment reten- 


tion only.) It gives a very good approximation of the height of sediment at the 
dam at various time intervals. 


CONCLUSION 


A comparison of the Area-Increment, Empirical Area-Reduction Method 
and actual reservoir resurvey area and capacity curves demonstrates the 
adequacy of the methods, (Figure 6). However, in the two methods demon- 
strated it is evident that a procedure in which all the major factors affecting 
sediment disposition is still needed. Also, several major factors affecting 
reservoir sediment distribution, such as density currents, delta development, 
and density of deposits with time, but not discussed herein, need additional 
study and evaluation to develop procedures for predicting the effect of these 
items on conditions in and upstream from the reservoir. 

Current work of the ASCE Task Force on Reservoir Sediment Distribution 
includes an evaluation of the effect of the various factors on sediment deposi- 
tion and distribution, and the nature of delta development. Investigations into 
these problems are also being carried out by agencies such as the Corps of 
Engineers, U. S. Geological Survey, and Bureau of Reclamation. Factors in- 
fluencing the sediment deposition in numerous reservoirs are being collected, 
consolidated, and analyzed, and mathematical relationships and criteria will 
be developed therefrom to give procedures that will yield consistently reliable 
results in predicting sediment disposition in proposed and existing reservoirs. 
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Table 1 


Carlsbad 
Alamogordo 


Stream-bed elevation at Das: 4150 Annual sediment inflow: 


Max. water 
Spillway c 


igina iginal edisent Revised 
Elevation Pec Capacity Ao Volume Area Capacity 
(feet) (Acre) (Acre feet) (Acre) (Acre feet) (Acre) (sere feet) 
i 2 3 6 7 
215 5 156, 150 235 15 100 


4270 
4250 
4230 


4220 


surface elevation: 4275 Period of sedimentation: 
rest elevation: Total sedigent for riod: 


AREA INCREMENT METHOD 


133,500 


235 
235 
68,750 235 


9T,000 


47,750 
32,750 


22 ,000 


4210 650 1b, 250 235 


&200 450 8,750 235 


*Determined new tero elevation 


Basic data 


4190 300 5,000 235 
#4164 235 3,300 235 
4180 200 2,500 200 
4170 
4160 
4150 


AREA INCREMENT PROCEDURE 


Ve = 24,580 ac.ft. 
H = 125 ft. 

Vo, ho, Ao = obtained by approximating ho 
and reading Vo and Ao from 

Original area - capacity 
curves (Fig. 12) 


(sediment volume) 


Basic equation 


assume ho = 25 feet 


then 24,580 = 235(125-34) + 3,300 


235 
The 
inc 


Ve = Ao(H-ho) + Vo 


Ao = 150 acres 
Vo - 1600 acre feet 


2h ,580 = (150)(125-25) + 1600 
24,580 16,600 


ho = 34 feet 
Ao = 235 acres 
Vo = 3,300 acre feet 


24,580 26680 


The area correction factor is therefore 
acres and the new zero elevation is 4184.0. 


resulting area and capacity data are shown 
Olugne 6 and 7. 


SEDIMENT DEPOSITION COMPUTATIONS 


3,600 ac.ft. 
6.8 years 
24 


23,510 


580 ac.ft. Checked b 


COLUMN DESCRIPTIONS 


April, 1958 


Computed by: CRM Date 
6-15-54 


JML Date 6-21-54 


3865 


2965 


2215 49,900 


18,600 


10,200 


415 


215 1,690 


65 290 


ie) 


ie) 0 


Basic data. 


Area correction factor as 
determined above. 


Sedigent volume as 
determined by average- 
end area formula, as 
based on area correction 

factor in Coluan 4. 


Revised areas obtained by 
Subtracting Column 4& 
Column 2. 


Revised capacity obtained 
by subtracting Column § a 
Coluan 3. 


| 
4100 110,000 
2450 18,810 
1750 p 235 16,460 1515 31,300 
1250 235 1&,110 1015 
yoo 235 11,760 665 
9,*10 
7,060 
4,710 
3,300 
2,500 
1,000 
250 
5. 
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Teble 
ALAMOGORDO RESERVOIR 


Sediment Deposition 
Empirical Area Reduction Method 


) 


volume 
(acre-feet) 
(acre-feer 


‘ 
v 
a 


4183: : : 
4163.5: : 


Sediment to be deposited « <4,600 acre-feet 
Standard Type II selected 


Column Descriptions 
From original area-capacity data. 
Relative depth--Ratio of contour interval depths to total depth (spillway crest t 
elevations in Column l. 

Ap--Values obtained from Figure 5 after proper stendard-type curve (Figure 3) has been selected. 
Enter curve with relative depth values tc obtain Ap values. 
Approximations of zero elevation at the dam after sedimentation. 
zero elevation et each increment will be equal to those 
elevation above zero elevation are obtain 
(Column 2) by the corresponding Ap value 

values at each contour elevation ebove. 


> stream bed) from 


Areas at and below approximated 

in Column 2. New areas for each contour 
z the original area at zero elevation 

nd multiplying this ratio (K) by the Ar 


Exarple: Assumec zero El. « 4190 
urface area = 300 acres 
st 


= 1.125 
= 266.7 


ote that zero elevation may not occur at an even contour interval (see 
Delta sediment volume determined from average-end ares formule: 


n/2 AL) = V 


supplement data) 


Summation of tnis column should equal (within 1 percent) the total sec 
Accumulated sediment volumes from Column 11. 
Revised areas obteined by subtracting values in Column 


2» 10 from those in Column 
Revised capecity obtuined by subtracting values in Column le from those in Column 


iment being deposited. 


* 
“aw e@* ° 4 © #2 a go. & 
(3) (3) (b) (5) 3(€) (7) : (8) : (9) 3; (20) (21) : : (23): (4) 
4270 «34,100 :133,500 3 2429 : 134 724,150 <:3,97 2109, 300 
$ : : 71,615 : 21,675 : 
4260 :3,200 : 97,000 : .& : 93 : 194 : 201 :2¢,505 :3,000 : 74,50 
: 68,750 : 31.090 :290 : : <27 : 237 : 720,315 :2,c10 
bebo 31,750 47,750 : :1.155 : : 247 : : 25¢ : 217,850 :1,490 : 29,900 
be30 «21,250 =: 3¢,750 : 21.235 3325 : omy : 267 215,235 : 980 : 17,52 
: : <2,000 .5€ 31-255 3334 : : 2A : : 271 712,545 : 630 : 
b210 ; 650 : 14,250 ; :1.240 : : 265 : 9,850 : 380 : 4,400 
«68,750 621-195 : : 256 : : 7,220 : 190 : 1,530 
$ 3 23,090 : 2,415 2,505 : : 
: 300 : 5,000 .32 3300 : : : : 243 : : 4,715 : fo : 285 
3 3 22,500 : 22,170 12€,¢15 : : 
4180 200 : ;: 31.0c 4 : : <00 : : 2,500 
: 21,500 : 71,500 : 21,500 : : : 
4170 ; 100 ; 1,000 : : 100 : : 100 : : 1,000 ; 
: 50 : 250 : ; .63 750 : : : : 50 : : 250 : >: 
2 : 3 3 : 250 : : 250 : 
H H : £9,509 : #4040 : $ 
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1, 2, 3- 
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€, 8, 10. 
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McILROY NETWORK ANALYZER* 
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SYNOPSIS 


Illustrates time-distribution of demand rates, design rates and equalizing 
storage requirements for certain Philadelphia districts. Efficacy of ; 
perfectly-balanced system characteristics and utilization of the McIlroy 7 
Analyzer for design in Philadelphia are discussed. Merits of different types 


of computers are compared. 


INTRODUCTION 


Inasmuch as the conditions and requirements of a given water system are 
generally somewhat unique, samples of data for some of the high service 
districts in the City of Philadelphia will be presented as a background for 
the descriptions of design practices and experiences which follow. After 
discussing design experiences with the McIlroy Network Analyzer, compari- 
sons will be made with digital and analog computers, from the standpoint of 
use by both utilities and consultants. 


Rates Used in Distribution Design and Their Justification 


In most of the discussion of design criteria in this paper, data from a 
recently completed study of the City of Philadelphia’s Roxborough High 


Note: Discussion open until September 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1588 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 2, April, 1958. 

*Presented at the Annual Convention of the ASCE, New York, N. Y., 
October 14, 1957. 
**Research Engr., Philadelphia Water Dept., Philadelphia, Pa. 


***Water Distribution Design Engr., Philadelphia Water Dept., Philadelphia, 
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Service District will be used as an example. The Roxborough High Service 
District is predominantly residential, with a moderate mercantile develop- 
ment and practically no industry. 

The average Roxborough consumption rates for each day in 1955 and 1956 
have been arranged in order of magnitude and appear in Figure 1. The 
“normal,” or demand equalled or exceeded half of a year, is very close to 
the yearly mean, or average, for both years. The maximum-day demand of 
record occurred in 1955; the 1955 data are generally more diverse. 

Hourly demand variations over a 24-hour period for both a typical day 
with an average demand equal to the 1955 average annual demand, and the 
maximum day of 1955, are shown in Figure 2. Here, the normal is definitely 
greater than the mean or average. The curves are similar in shape. Note 
that demands included between 0.8 and 1.2 of the average demand in both 
cases are representative of about 12-hours (50%) of the 24-hour period, and 
rates above 1.2 of the average demand for each case are representative of 
only about 6-hours (25%). 

Figure 3 is to some degree a composite of the previous two figures, for 
1955. In addition to the daily average demand for the 365 days, the maximum 
and minimum hour demands have been graphed. To be especially noted is 
the fact that variations in daily average demands within 90% to 110% of the 
average annual demand cover 83%, or about 300 days, of the year. 

Because of the moderate variations about the average annual demand for 
the major share of the year, design has been based primarily upon the maxi- 
mum and minimum hour demands of the maximum and average day. Since 
Philadelphia is supplied entirely from river sources by pumpage (annual 
average consumption, 370 MGD), efficient pumping at average day rates is 
particularly important in all districts. The maximum hour of the maximum 
day is normally the critical demand in terms of service pressures at points 
of high elevation. 

The variations, daily and hourly in Figures 1 to 3, give a good indication 
of intermediate pump needs for adequate coverage of demand fluctuations. 


Uniform Pumping Storage Requirements 


In an effort to provide more uniform service pressures and to improve 
supply reliability, equalizing storage is planned for all high service districts. 
One of the major objectives in providing equalizing storage is to permit 
pumping rates as near as possible to constant or uniform 24-hour rates. 

Figure 4 shows the storage requirements for five districts recently 
studied. Average day and maximum day storage requirements are similar 
in trend. In Figure 4(b) the greatest equalizing storage requirement shown 
is 14.8% of the maximum day demand. The zero storage pumping ratios 
represent maximum hour demand rates. In contrast, Schmid(1) reports a 
25% uniform pumping equalizing storage need for Royal Oak, Michigan, for 
a maximum day demand of 26.4 MGD (where the maximum hour demand is 
1.99 times the maximum day demand). Percentagewise, much less storage 
is needed for uniform pumping in the Philadelphia districts shown. 


Utilization of the McIlroy Analyzer in Distribution Design 


There are approximately 2900 miles of pipe in the system, ranging in size 
from 4'' to 93''. However, the McIlroy Analyzer of the Philadelphia Water 
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Department, installed in January, 1956, has been used mostly in connection 
with arterial mains. Field measurements of Hazen-Williams C-values for 
arterial mains are made by the Pitometer Company. To date, the rate of ob- 
taining field measurements has set the pace in rehabilitation design. 

Distribution rehabilitation design includes piping, equalizing storage, 
pumping stations, routing and districting. Many design variables and sys- 
tem complexities become involved. For example, in the Roxborough High 
Service District previously mentioned, the existing district is presently fed 
by one pumping station; revisions include two standpipes recently constructed 
and a third future standpipe at another location, in addition to a second, en- 
tirely new, pumping station and a modification of the district boundary. 
(Many of these revisions are necessary for abandonment of an existing filter 
plant). 

The following conditions were investigated using the analyzer for present 
and future demand rates in the Roxborough High Service arterial system: 


1) Minimum and Maximum Hour of Maximum Day - latter usually deter- 
mines reinforcement needs (new mains or cleaning and lining). 


2) Minimum and Maximum Hour of Average Day. 
3) Arterial pressures for fire flows. 
4) Pumping station outages. 

For present demand rates, also investigated were: 


5) Effects of premature change in district boundary (prior to new, second 
pumping station), and 


6) Schedule for satisfactory service under which mains can be laid or 


cleaned and cement lined, connections can be made and valves installed 
or rehabilitated. 


Uniform pumping was one of the major design objectives, from the stand- 
point of improved service, efficiency and amenability to automatic control. 
Storage volumes and rates were computed for perfect 24-hour uniform pump- 
ing. Using constant pumping rates and the computed storage water levels, 
pumping station discharge heads were obtained with the Analyzer (after piping 
system had been designed). The resulting discharge heads for one of the two 
Roxborough High Service Pumping Stations is shown in Figure 5, plotted in 
terms of pump T.D.H. Hypothetical pumps are shown for the purpose of illus- 
tration. On the maximum hour of the average day (lower point) a greater 
quantity than the average rate will be pumped, with a lesser quantity than 
computed coming out of storage. On the minimum hour of the average day 
(upper point) a lesser quantity than the average rate will be pumped, with a 
lesser quantity than computed going to storage. Adjusting the storage curve 
in this instance would not noticeably change the storage volume used, and 
would modify only slightly the water levels computed for uniform pumping. 

If the curve for pump No. 1 was raised, then less than the storage volume 
computed for uniform pumping would be used; for a lowered curve, more 
storage volume would be used. For the maximum day under practically all 
conditions shown, the pumping rate would be less than the average of 9.5 MGD, 
and the difference would have to be made up by the second pumping station or 
by a storage depletion. 
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In Figure 6 is plotted part of the maximum day storage output-versus- 
water-level data for three Roxborough High Service standpipes, computed 
for uniform pumping and assuming a specific filling schedule. The Maximum 
Hour rate out of storage with uniform pumping is 8.0 MGD. In the table below 
the graph, are the T.D.H. for pumps of each of the two stations, as designed, 
for uniform rates on the Maximum Hour of the Maximum Day (the second and 
third cases are outages of each of the two stations). In the “initial balance” 
columns are shown T.D.H. and rates for the sources indicated, which were 
obtained by extrapolating system test points to the pump curves; the storage 
rates shown would, by default, be required to balance the rates. A revised 
storage rate requires a new storage curve. Since uniform pumping can be 
approached but seldom perfectly attained, the only way in which a perfectly 
balanced system analysis could be achieved would be by making a series of 
runs starting at an assumed full storage time, following the demand curve 
hour by hour and balancing heads with pump characteristics. This is much 
too involved for the minor modifications which would be obtained in this 
instance. The balancing of pumped systems with and without equalizing stor- 
age will be discussed later in terms of computer operation and analysis. 


Other Computers for Analysing Distribution Systems 


With the advent of the McIlroy Analyzer and the successful programming 
of digital computers for network analyses, any discussion of approximation 
methods should be confined to the simplest of distribution systems. 

Features of the McIlroy Analyzer have been presented previously.(2,3,4,5,6) 
Suffice it to say that the principal feature of the McIlroy Analyzer is a system 
of resistances through which voltage drop (analogous to head loss) is propor- 
tional to the 1.85-power of the current (current flow is analogous to flow rate). 
The McIlroy Analyzer thus provides a “model” or direct, working analog of 
the system under study. In operating the Analyzer, “drawoffs” are set at 
prescribed full load values and system inflow rates (from pumps or storage) 
are brought up from zero to system capacity. Often, in design analyses, 
needed system revisions are evident at inflow rates well below the design 
rates; obviously needed piping improvements are then made prior to any 
formal or complete test runs. Trial changes can be made in a matter of 
minutes. For example, changing the distribution of inflow between two pump- 
ing stations can be accomplished by adjusting two or three dials. The most 
valuable attribute of the McIlroy Analyzer in this respect is the elimination 
of unrealistic or superfluous or discordant arrangements with quick, simple, 
direct dispatch. 

Electronic digital computers have been employed in gas and water distri- 
bution network analyses. The punched card program for a low-pressure gas 
system using a digital computer, with a wired program, has been presented 
by Sickafoose.(7) (However, a 2.00-power rather than 1.85-power is involved). 
An outline of the programming for a water network analysis using a stored 
program type of digital computer has recently been presented by Hoag and 
Weinberg.(8) In both of these references, a form of the Hardy Cross relaxa- 
tion technique was included in the analysis program, the computer perform- 
ing the trial-and-error iterations. The digital computer thus adjusts flow 
rates initially assumed by the designer for each pipeline into a balanced 
distribution to almost any degree of desired accuracy in loop headloss 
balance. Other tests of the same system can often be made by replacing only 
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a few of the punched cards from the initial test data instructions. The ma- 
chine itself is tied up only a very few minutes in any one run. However, 


tabulation and manipulation of data instructions is tedious and time-consuming. 


Further, the results are normally presented as tabulated head loss and flow 
rates (with direction) for each pipeline, and these must be summed 
algebraically before the analysis can be evaluated. An error in programming 
or data instructions is not evident until the run is made and the analysis is 
completely summated and checked, and can be very difficult to locate. 
A punched data card is needed for each loop (identifies mains and flow 
directions in mains comprising the loop), and for each main (“C”, diameter, 
length, initially assumed flow rate) for each loop (mains often common to two 
or three loops). About 50 to 60 punched cards are needed for the basic pro- 
gram instructions, almost regardless of the size of the network being studied; 
the basic program instructions can be used over and over with different net- 
works. The possibility of errors in data instructions is evident, particularly 
with medium to large networks. Further, the “model” visualization feature 
of the McIlroy Analyzer is completely absent from the digital computer. 

Analog computers are being used in water-hammer analyses,(9,10) but 
devices for the simulation of the non-linear characteristic of pipelines in an 
analog circuit are cumbersome, and network analyses for either a steady or 
a transient state must of necessity be restricted to simple systems. As with 
the McIlroy Analyzer, meter adjustment and wiring of components must often 
be charged against “machine time,” since in both cases the machines are 
usually not then available for other work. (The McIlroy Analyzer is a 
specialized type of steady flow, or passive, analog computer). 

Network analyses can be run, by a series of trial runs with intermediate 
computations and adjustments, on electric service power-system calculating 
boards with linear resistors. (Refer to first four references listed in 


Reference (3) ). Only a very few of these boards have been used for water 
network analysis. 


Analyses with Balanced Pumping Station Characteristics - 
Without Equalizing Storage 


Shown in Figure 7-A are two hypothetical pumping station pump charac- 
teristics. For direct pumpage the demand would be supplied by the combined 
output of the two stations. For any given pumpage rate the discharge hydrau- 
lic gradient, h, of a given pump is fixed by its characteristic curve. For an 
output QR there is a singular hp for a given Roxborough pump. Since Qw 
must equal the demand less QR, by default there is a specific hy for the 
pump delivering Qw. Hence, for a given demand rate, either Qa or Qw 
is a function of hy-hp, or hp-hy. As an illustration, for an 18.9 MGD 
demand in Roxborough High Service (1955 Maximum Date Rate), the com- 
bined pump characteristics for the two stations have been plotted in 
Figure 7-B. (A maximum difference of 20-feet is the practical limit in 
this district. “Pumps b and d” represent the only realistic combination 
in this illustration, if only from the standpoint of efficiency). 

Information of this type (Figure 7-B) could be fed into a digital or analog 
computer, and provided that the system including pump characteristics has 
the capacity of being balanced, a perfect balance would be achieved. With the 
current model of the McIlroy Analyzer, balance can be achieved only by suc- 
cessive meter adjustments, converging on pump characteristic balance 
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curves similar to Figure 7-B. From observations of direct-feed low-pressure 
gas network analyses using a McIlroy Analyzer, achievement of a balanced 
system is not particularly difficult. Experience with direct-pumped system 
analyses at the Philadelphia Water Department is limited inasmuch as all 

high service areas have or will have equalizing storage, and the only district 
which will be served by direct pumpage in the future has not yet been 

analysed. 

Where the design of new pumps is at issue, balance curves can be derived 
from network analyses and pumps selected to match or approximate the 
balance curves. 

Analysis of a multiple gravity-supply system without equalizing storage 
would be similar, but more complicated if the water levels varied over a 
wide range. 


Analyses with Balanced Pump and Equalizing 
Storage Characteristics 


In Figure 5 it was noted that a program of uniform pumping with equalizing 
storage would not match perfectly a set of typical pump characteristic curves. 
If the pumps and the equalizing storage are existing system components (not 
a subject of design) and the system is inherently capable of balance (a very 
restricted case) a step-by-step integration of the demand curve andthe _ 
presentation of an automatically derived storage curve could be achieved 
with both a digital and an analog computer. The normal analog computer 
would be restricted to a simple system because of the complex assemblies 
necessary for each pipeline. The digital computer would be limited, in ex- 
tent of complexity of system to be analysed, to the patience and exasperation 
limits of the design engineer and to the ingenuity of the programmer. Use of 
a digital computer for a balanced system with pumped supply and equalizing 
storage has not been mentioned in the literature. 

Step-by-step integration of a 24-hour demand curve using the McIlroy 
Analyzer would be tedious and would require considerable meter adjustment 
to get each hour of demand in balance. 

In general, when the design of either the characteristics for a pumping 
station or for equalizing storage are under investigation, the necessity or 
desirability of simultaneously achieving a perfectly balanced system is 
highly questionable. This is particularly true when new pumping stations 
and new equalizing storage facilities are both to be designed. 

In attempting a perfect system balance, using any type of machine, some 
very pertinent questions must be raised: 


1) Is the field data on pipe lengths, diameters, locations and C-values 
practically perfect? 


2) Are drawoffs as evaluated by population indices and field flow-data 
typical and near-perfect? 


3) Will or can the pumps, tanks and other improvements defined by the 
designer be exactly reproduced in procurement and construction? 


4) Will the system always be operated on identically the same schedule as 
presumed by the designer? 


Inasmuch as righteous affirmative answers to all these questions will lack 
complete conviction, attainment of a perfectly balanced system analysis is 


ASCE WATER DISTRIBUTION DESIGN 1588-7 


often not justified. The objective here is to promote a quality or level of 
analyses reasonably consistent with the quality and accuracy of the data. 
Therefore, attempting to attain a perfectly balanced system with a digital 
computer would seldom be warranted despite the capacity of the machine to 
handle the problem. Another factor worthy of consideration is the large seg- 
ment of additional design time involved in balance refinements, using any 
type of machine. 


Merits of Various Machines for Network Analysis 


The adaptation of electric service power-system calculating boards with 
linear resistors, by means of computation and manual adjustment iteration, 
is not very satisfactory. The calculating board is similar to a digital com- 
puter to the extent that the final result is all that is obtained; it lacks most 
of the analog features of the McIlroy Analyzer. More important, the operating 
staff of only a very few of these calculating boards have run water network 
analyses. 

The current types of analog computers available for analyses are too 
restricted in pipeline capacity to handle typical complex network problems. 
Until a small and simple non-linear resistance component is developed for the 
analog computer, exploitation of this excellent tool will be quite limited in net- 
work analysis. 

There are at present only 15 McIlroy Analyzers in use, 7 of which are fully 
employed on gas network analyses, one on air and only 3 on water. Only 4 
machines are readily available to private consulting engineers. 

However, various models of advanced types of digital computers are in 
use in most metropolitan, business and university areas. Because of its 
more widespread availability, the digital computer is the machine which 
holds the most promising future for use by private consultants. As soon as 
a standard set of instructions is published (particularly for the stored- 
program type of computer), together with an outline for setting up data for 
problem mains, loops and special operating instructions, the digital com- 
puter will become an extremely valuable design tool for all consultants. 

With a set of standard program instructions available, the time required 

of the designer to set up and evaluate a test run should not differ appreciably 
from the time required for the same test run on a McIlroy Analyzer, pro- 
vided that special instructions (such as balanced pump characteristics) are 
not too involved. If few runs are made in a given study, or if the designer 
has little experience with the procedure, it will no doubt take more of his 
time to do a test run on the digital computer. A Texas computer firm(11) 
will perform digital computations via correspondence. This service is a 
step in the right direction. Waiting time between test runs will not be re- 
duced until either similar services are available across the country or until 
a standard program is published for general use. Analysis by means of the 
McIlroy Analyzer has outstanding advantages over analysis by digital com- 
puter, if the design engineer is present during the test run. There is, on the 
other hand, no particular advantage to the design engineer in being present 
during the performance of iterations by a digital computer. 

Opposed to the private consultant in general practice to whom a network 
analysis is only an occasional problem, is the distribution designer or planner 
of a public utility. It is the steadfast opinion of the authors that the McIlroy 
Analyzer is the superior device for distribution design of gas or water 
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networks when such work constitutes an important part of engineering 
functions. 

The analog features of the McIlroy Analyzer are outstanding where dis- 
tribution design is a vital or full-time occupation. The design engineer gets 
the “feel” of the whole system as the required demand is approached. Sys- 
tem revisions can be made in a matter of minutes. Decisions concerning 
revisions or alternate schemes can be made most effectively by immediate 
appraisal of the results as modified or evaluated in terms of observed, 
mentally recorded trends. Proportionate head-loss in mains is indicated 
by the comparative brightness of the resistors. Thus the resistors and flow 
and head meters provide a special auxiliary service as instantaneous, simple 
visual aids. 

With the digital computer, on the other hand, the designer gets only a 
tabulated summary which is meaningless until analysed in toto. The elapsed 
time between conversion of the problem into simple arithmetic and interpre- 
tation of the final results contributes nothing to the problem evaluation. Pos- 
sible revisions or alternate schemes are not readily suggested. Physical 
identity with the problem is lost and the solution is reduced to abstract num- 
bers in tabular form. 

The Mcllroy Analyzer of the Philadelphia Water Department is used an 
average of only one or two full days per week. The remaining work week is 
devoted to studying and evaluating data for analyses (about 50% of rehabilita- 
tion design work-load), preparing reports, advising on small main local grid 
changes, and additions consistent with overall distribution planning. The 
Water Distribution Design Unit is keeping up with the Pitometer Company 
C-value surveys and the capital budget program. The machine could be used 
more if the 2-1/5-man staff was increased, but to double machine use it would 
be necessary to double data analysis and the whole operation would become 
quite involved. 


CONCLUDING REMARKS 


To repeat, the accuracy of network analyses, regardless of type of ma- + 
chine employed, is of the order of the accuracy of the data fed into the ma- 
chine. The value of the results obtained depend entirely upon their proper 
utilization and interpretation by the design engineer. 

The McIlroy Analyzer was tailor-made for Water Department network 
distribution design and is heartily recommended for use by medium and 
large water departments. The more prevalent digital computers will 
probably become increasingly adapted to network design by private con- 
sultants and will probably be extensively used as soon as general instructions 
are made available and until such time as there are a much greater number 
of McIlroy Analyzers available on a rental basis. If a simple pipe-simulating 
component for analog computers is developed, water-hammer characteristics 
of networks can then be determined with little difficulty and vastly reduced 
approximation. 

As evidence of the progress made in water network analysis procedures, 
one need only consider that as recently as two years ago the question at issue 
was the comparison of “hand” trial-and-error calculations versus utilization 
of the McIlroy Analyzer. However, the gas industry remains well ahead of 
the water field. Six of the monthly issues of Gas, in 1955, gave pertinent 
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discussions related to a series of symposia on the use of various computers 


in gas distribution network piping design (particularly January, May and 
September). 
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SYNOPSIS 


A review of important work performed in the field of sediment densities is 
given in this paper. The factors affecting the initial density of newly deposited 
sediments and the rate of consolidation of deposited material are discussed. 
Results of extensive sampling of deposited sediments are summarized, and 
methods of estimating sediment densities for particular sediment conditions 
and methods of reservoir operation are presented. Consolidation theories 
developed in the field of soil mechanics are presented, with a discussion of 
their possible application to the sediment problem. 


Purpose 


This review is presented as a contribution to the work of the Task Force 
on Reservoir Sedimentation, Hydraulic Division, ASCE, which is making a 
review of work in the field of reservoir sedimentation. 

The purpose of this paper is to present a review and summary of the work 
that has been done in the determination of densities of deposited sediments. 
Some of this work has been presented as incidental to sediment studies of 
more general nature and, therefore, has not been given the attention it de- 
serves. Other work has been presented in Government publications that are 
not widely available to the non-Government engineer. Still other work is 
recorded in literature of other fields of engineering which is not widely read 
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by engineers working in the field of sedimentation. No original contributions 
are made in this paper. 


The Problem 


The determination of the volume of sediment expected to accumulate in a 
reservoir is a problem which occurs in the planning of all reservoirs. Where 
the expected sediment volume is large, it becomes necessary to set aside a 
reasonable amount of space for sediment storage, so that the purpose for 
which the reservoir is intended can be served without undue shortages during 
the useful economic life of the reservoir. Even when the sediment volume 
anticipated is so small that it constitutes no problem it is still necessary in 
most cases to demonstrate proof, through studies or verification by measure- 
ments. 


Estimates of sediment accumulation in reservoirs consist of three phases: 
(1) Estimate of sediment inflow to the reservoir 

(2) Estimate of sediment inflow that is retained in the reservoir 

(3) Conversion of retained sediment from a basis of weight to volume 


Much attention has been devoted to the first phase, the estimate of sedi- 
ment inflow. Intensive research has been performed on development of 
proper instrumentation for measuring sediment load in the field, and for 
analysing the sediment characteristics in the laboratory. Much effort has 
been expended in the routine observation of a network of sediment sampling 
stations throughout the country and in processing the data obtained for gen- 
eral use. Many studies have been made of the proper method of interpreting 
these data to arrive at an estimate of the sediment inflow. The frontal attack 
made on this phase has been very essential, very effective, and has justified 
its cost. 

The second phase involves an estimate of the trap efficiency of the reser- 
voir. This phase was covered in another paper sponsored by the Task Force. 
The third phase is the subject of this paper. 

The sediment inflow that is retained in the reservoir is determined in 
terms of weight, since all field observations are made on the basis of sedi- 
ment concentration by weight. To complete the estimate of sediment accumu- 
lation in a reservoir, there remains the very short, but highly important, step 
of converting the weight of the retained sediment inflow into volume. This re- 
quires an estimate of the space occupied by a unit weight of sediment, which 
is variously referred to as specific weight, unit weight, or density. In this 
paper, the term “density” is used to define this conversion factor. Density 
is expressed as the weight of solid material per cubic foot of mixture of 
solids and water. 

If the density of the deposited sediment is not correctly estimated, the es- 
timated volume of deposited sediment will be in error by the same percentage 
as the error in the estimated density. The density of deposited sediment in 
reservoirs has been observed to vary from 18 to 125 lbs. per cubic foot. 
Therefore, it can be seen that the determination of a proper density value 
is quite important. Unfortunately, while many observations have been made 
of density under various conditions, relatively little concentrated study has 
been devoted to relating density to the various casual factors. 
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Equipment for Measuring Sediment Density 


This discussion of equipment is not intended to cover all the types of 
samplers that have been developed for collection of samples of deposited 
sediment. A description of a number of mechanical samplers is given by 
a U. S. Government interagency publication.(1) Two pipe samplers developed 
by the U. S. Bureau of Reclamation are typical of the mechanical samplers. 

The first pipe sampler was developed about 1947 and used in a number 
of reservoir surveys through 1949. It consisted of a one and one-half-inch 
diameter pipe about eight feet long with a check valve near the upper end and 
threaded at the bottom end for attaching sample containers, which were about 
eight inches in length. A rope was attached to the top of the sampler for 
lowering it through the water into the sediment at the reservoir bottom. 

The second sampler was designed and fabricated in 1950. It consisted of 
a standard two-inch diameter pipe with a removable cutter head at the bottom 
for soil penetration, and a check valve and removable lead collar weights at 
the top. Either a five-foot or ten-foot pipe section may be used. A plastic 
tube inserted in the pipe sampler is used to secure the sample. This sampler 
is dropped over the side of the boat and allowed to fall to the reservoir bot- 
tom. 

The latest instrument developed by the Bureau of Reclamation is the 
radioisotope densitometer, which was developed and tested under the sponsor- 
ship of the Sedimentation Subcommittee of the Federal Interagency Committee 
on Water Resources.(2) It consists of a probe containing cobalt 60 as a source 
of gamma rays and a self-ionization chamber for the detector. A lead cylin- 
der is provided for radiation shielding and weight to drive the probe into the 
sediment. To utilize the phenomenon of Compton scattering of gamma rays 
for sediment density measurement, the radiation scattered from the source 
of gamma rays back to a detector is measured. When the probe is placed in 
the sediment, the intensity of the radiation scattered by the surrounding 
media to the detector is a function of the electron density of the saturated 
sediment. Preliminary field tests of this instrument have shown that the 
basic method of measuring sediment density with gamma ray scattering is 
sound and practical. Further tests of the instrument are contemplated. 


Factors Influencing Sediment Density 


It can be assumed, with reasonable accuracy, that all sediment particles 
have the same specific gravity, about 2.65. The density of deposited sedi- 
ment thus depends upon the amount of void space occupied by water. This 
is influenced by the size of material, degree of sorting and gradation, and by 
the amount of consolidation that has oceurred. 

The rate of consolidation in itself varies with a number of factors. Engi- 
neers in the field of soil mechanics have developed well proven theories for 
consolidation of continuously submerged materials. However, engineers in 
the field of sedimentation have not yet solved the problems of compaction of 
either submerged or exposed sediments. 


Definition of Sizes 


Numerous classifications have been used to describe the various size ranges 
of sediments. Except where otherwise specified, the classification of the 
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U. S. Bureau of Soils is used in this paper. These are as follows: 

Sand - particles with sizes between 0.05 and 1.0 mm 

Silt - particles with sizes between 0.005 and 0.05 mm 


Clay - particles with sizes less than 0.005 mm 


Variation in Observed Densities 


Many samples of deposited sediments have been collected for purposes of 
determining the density of the deposits. These samples have covered a wide 
range of sediment characteristics and conditions of deposition. Comprehen- 
sive summaries of these data are given by Lane and Koelzer(3) and by 
Hembree, Colby, Swenson, and Davis.(4) The first named summary presents 
data on over 600 samples collected at nearly 100 locations by various investi- 
gators throughout the world. Hembree, Colby, Swenson, and Davis gave data 
on about 225 samples from 38 locations in the United States. The variation 
that exists in sediment densities is demonstrated by the data in Table 1, 
selected from these summaries. These data indicate a general relationship 
between particle size and density. However, it can be seen from the com- 


parative data on similar size classifications that factors other than size in- 
fluence the density. 


Consolidation Theories 


In the field of soil mechanics, consolidation theories for submerged 
materials have been well developed, starting with the pioneering work of 
Terzaghi. Consolidation of submerged material takes place primarily be- 
cause of imposed load. In the consolidation process, water must be squeezed 
from the deposits; therefore, the permeability of the deposits is a significant 
factor. Thus permeable sands have a higher rate of initial consolidation than 
impermeable clays. 

Terzaghi first advanced the theories under which consolidation under 
constant load varies with time. Figure 1 shows the results of typical 
laboratory tests, carried out at Swarthmore College, which demonstrates 
Terzaghi’s theories. Many similar curves are given in standard references 
in soil mechanics literature. The percent consolidation is the amount of 
compression at any time divided by the total primary compression. Primary 
compression is that which results from squeezing out the excess water from 
the voids as the load is being transferred to the soil particles. The time of 
primary settlement varies from a few seconds for coarse material to years 
for very fine material. Secondary compression starts when all the load is 
taken by the soil particles. The rate of secondary consolidation is extremely 
low, with settlement being due primarily to readjustment of the particle 
structure. 

The relationship of imposed load to consolidation has been developed 
through years by many engineers working in the field of soil mechanics. 
Figure 2 demonstrates the relationship between pressure and voids ratio, 
for clay and sand sizes. The voids ratio is the volume of voids divided by 
the volume of soil particles and is thus inversely proportional to density, 
assuming a given specific gravity. A greater degree of consolidation occurs 
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with clay than with sand, but clay sizes maintain a higher percentage of voids 
than the coarser sand sizes. These standard curves developed by soils me- 
chanics engineers appear to be significant to the sediment engineers. Con- 
clusions similar to those reached from these curves have been reached by 
Parker Trask, in the field of geologic sedimentation, and by Russel and 
Dickey, in investigating petroleum shales. 

Trask(5) states that compaction of sediments, with resulting decreases 
in porosity, is clearly related to the imposed load. For sediments to com- 
pact, the water must be squeezed from the deposits; therefore, the permea- 
bility is also a big factor. Thus, permeable sands have a high initial rate of 
compaction compared with impermeable clays. All sizes compact quickly at 
first, and then the rate of compaction declines. Trask concludes that the es- 
sential factors in compaction are load, grain size, permeability, and time. 

Russel and Dickey(6) discuss the compaction of fine silts and clays, stating 
that these sizes, when settled out of still water, have porosities greater than 
75 percent, the porosity increasing with decreasing particle size. The fine 
character of the particles causes them to settle into a cellular structure. 

As the pressure of the overlying sediments increases with continued deposi- 
tion, compaction occurs through collapsing of the cells and expulsion of the 
water. Soft clays are stated to have high porosities, but very low permeabili- 
ties, so that the water is expelled slowly. They show curves which indicate 
an increase of density with depth of deposits, in petroleum shales from 600 
feet to 5,000 feet of depth. 

Happ(7) states that, in alluvial sediments, there is no doubt that air drying 
is the major factor in determining density of unsubmerged sediments, at least 
of the finer sediments. He cited cases of two samples of similar fine 
material, one collected a day or so after exposure and the other after at 
least four months of exposure. Density of the first sample was 48.3 pounds 
per cubic foot, and that of the latter sample was 94.2 pounds per cubic foot. 

From the above, it is apparent that the primary factors influencing the 
rate of compaction of particles of a given size are, (1) the weight of overlying 
sediment, (2) the degree of exposure to drying, (3) particle size, (4) permea- 
bility (probably related to particle size), and (5) time. 

Some of the foregoing theories are verified by the results of the compre- 
hensive sedimentation survey of Lake Mead, carried out in 1948-1949.(8) In 
this survey, the deposited sediments were intensively sampled. The Lake 
Mead deltas are geographically graded, the coarsest particles accumulating 
at the head of the deltas, with the finest particles being carried further out 
into the lake. As the delta has built downstream, new beds of coarse material 
have deposited over old beds of fine materials in some areas. The sand de- 
posits, although almost continuously submerged, are firmly packed. On the 
other hand, the silt and clay deposits have an extremely high water content 
which gives them the consistency of a soft ooze. Only the samples of the 
most deeply buried layers are solid enough that they do not flow of their own 
weight. There are also numerous gas cavities in some zones of the finer 
deposits. The water content of the new finer deposits is commonly more 
than 80 percent. As additional material is deposited, part of the water is 
driven out. Water content in a vertical section has been observed to vary 
from 16 percent in the sand at the top to 75 percent in the silt and clay of 
the bottom. The water content in all sections decreases with increasing 


depth. The density has been observed to vary from 18 to 108 pounds per 
cubic foot. 
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Initial Density 


The sedimentation engineer is interested in the variation in density 
throughout the economic life of the reservoir. He, therefore, wants to know 
the initial density, the density at the end of a given period, say 100 years, 
and the average density throughout the period. From his viewpoint the 
density at the end of any reasonable short time, say up to a year, can be 

' considered initial density. Since primary consolidation of large sized 
material has been demonstrated to occur very rapidly, with only minor 
secondary consolidation occurring thereafter, the data from practically all 
samples of coarse material can be assumed to represent initial density, 

a regardless of the time the material has been deposited. This will not be 

true of the finer materials. Initial densities of sediments, as determined by 

various investigators, are shown in Table 2. 

The values of density developed by the California State Division of Water 
Resources (9) have gained considerable acceptance for coarse sediments. 
These values were developed after a very thorough study of the voids and 
water retaining capacity of sediments, in which more than 200 samples were 
analyzed. It should be noted that the sizes are identified by the maximum 10 
percent size range, whereas all other data in the table are in terms of the 
complete size range or the median size. 

Trask(10) separated sediments into various size ranges and determined 
their density after what he termed a “very short settling time.” The data 
of Hembree, Colby, Swenson, and Davis are from samples collected near the 
surface of submerged sediment deposits, which were believed to be formed 
within a few years of sampling time. 

Data for 65 samples taken of alluvial deposits in the Middle Rio Grande 
Valley are recorded by Happ. Although he found great variation for similar 
sizes, he found a general relationship in which density increased with median 
particle size. He attributed most of the variation to differences in compac- 
tion, and drew envelope curves which he believed represented the initial and 
compacted densities. His values in Table 2 are taken from an envelope curve 
of minimum densities which represent what he believed to be initial densities. 

A comparison of the laboratory values determined by Trask with the field 
samples determined by the other investigators indicates higher densities for 
the field samples. This could be due to some consolidation occurring prior 
to the field sampling, even though the deposits were assumed to be relatively 
new and unconsolidated. 


Variation of Density with Depth of Deposit 


The variation of density of marl sediments with depth of deposit was 
determined in Lake Niedersonthofen, Bavaria, By Reisinger.(11) These ob- 
servations reflect the combined effect of consolidation due to time and the 

weight of overlying sediment. The results are shown in Table 3. 


Variation of Density with Drying 


Happ states that his envelope curves are representatives of the approxi- 
mate magnitude of the factor of compaction by drying in the Middle Rio 
Grande. His results are shown in Table 4. The higher compaction shown 
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- VARIATION OF DENSITY WITH DEPTH OF SEDIMENT 


IN LAKE NIEDERSONTHOFEN, BAVARIA 


Depth 
(meters) 


Upper Layer 


WON OUP Wh 


10 


Initial Density 


(lbs. per cu, ft. ) 


85 
78 
72 
68 
59 


48 


Density 
(lbs. per cu, ft.) 


21.6 
39,6 
41.5 
34.8 
42.1 
45.1 
52.8 
53.9 
56.0 
73.6 
74,0 
80,4 
85,2 
97.7 
92.7 
84.0 
87.2 
85.3 
92.8 
88.1 
89.6 


(BY HAPP) 


Terminal Density 
(lbs. per cu. ft.) 


104 
102 
98 
97 
95 


94 


Ratio of Compacted 


to Initial Density 


/ 
| 
11 
12 
13 
14 
15 
16 
17 
18 
20 
TABLE 4- ESTIMATED COMPACTION OF MIDDLE RIO GRANDE SAMPLES 
1, 30 
1. 33 
1, 36 
1.43 
1.61 
|_| |_| 1,98 
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for smaller sizes bears out the consolidation curves shown on Figure 2. 


Consolidation of Sediment with Time 


Lane and Koelzer studied the sediment density data available and developed 
a formula relating density to particle size, time, and method of reservoir 
operation. Their classification of reservoir operations was as follows: 


(a) Sediment always or nearly always submerged - includes reservoirs 
normally kept full and emptied only for a short time, at infrequent intervals. 


(b) Normally a moderate reservoir drawdown - includes cases where 
reservoirs store over long periods, with small drawdown in normal years. 


It also includes run-of-river hydroelectric plants, where the basin is used 
for pondage. 


(c) Normally considerable drawdown - includes reservoirs with relatively 
large drawdown in normal years. 


(d) Reservoir normally empty - includes retarding basins and reservoirs 
that are empty in normal years. 


The Lane and Koelzer formula is shown in Table 5. Values of Wy were 
determined by examination of over 600 samples, giving consideration to 
degree of drying considered compatible with the method of reservoir opera- 
tion. Values of K for the first three methods of operation were chosen to 
give the same density over a 1,000 year period as would be reached under 
condition (d) in one year, where the reservoir is normally empty and the 
sediment would have ample opportunity to dry out. The logarithmic correc- 
tion for compaction was used to give a rapid consolidation in the first few 
years, with rapidly decreasing rates later. This follows the general trend 
found by soils mechanics engineers, although it seems probable that equa- 
tions could now be developed which more closely approximate the primary 
and secondary settlement rates shown on Figure 1. For material covering 
a range of particle sizes, it was recommended that the equation be used with 
appropriate weights applied in proportion to the percent by weight in each 
size classification. An illustration of this computation is shown on Table 5. 

Miller(12) of the U. S. Bureau of Reclamation has made studies which he 
suggests indicate that the initial densities in the above formula are too high. 
He believes that the values determined by Trask, as shown in Table 2, are 
more applicable, particular for the finer sizes. Trask’s values are for the 
always submerged condition. It should be noted that Trask’s values are for 
density after “a very short settling time,” in contrast to density values of 
Lane and Koelzer which are represented as being after one year. 

The Bureau of Reclamation uses Miller’s modification, considering that 
the “always submerged” condition is more representative of initial densities 
and that Trask’s values for clay sizes are more applicable. 


Methods of Computing Average Sediment Density 


The data presented lead to computation of sediment density by applying 
specific densities to different size grades. Miller presented one method and 
Hembree, Colby, Swenson, and Davis presented two methods of applying these 


ASCE DEPOSITED SEDIMENT RATES 1603-13 


TABLE 5 
LANE AND KOELZER FORMULA 


1_g10T 


W «density after T years in lbs. per cu. ft. 


W) = density after 1 year in lbs. per cu. ft. 


K =a constant for each sediment class and operation 
condition to reflect consolidation 


T = number of years of consolidation 


Values of W and K 


Sediment always submerged 
or nearly submerged 


Normally a moderate 
rescrvoir drawdown 


Normally considerable 
reservoir drawdown 


Reservoir normally empty 


Example of computation of density of deposit "always submerged" 


Deposit consists of sand 
silt 
clay 


20 percent 
40 percent 
40 percent 


Density of deposit after 100 years consolidation ~ 0. 20 (93) 


+ 0.40 (05+5.7 log)g 100) + 0.40 (39 + lo log}g 190) = 74.0 lbs. /cu. ft. 


MILLER EQUATION 


Average density of all deposits after T years of compaction = 


T 
Wave = W, + 0.434 K (loge T) 
T-1 


/ 
where 
f 
Reservoir Operation [wie | 
a. 
93 0 65 5.7 30 16.0 
93 0 74 2.7 46 10.7 
93 0 79 1,9 | 60 6.9 
93 | 82 | 0.0] 78 | 0.9 
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data to suspended sediment data to arrive at an overall value representing 
the average for the suspended sediment load. 

The Lane-Koelzer equation gives the density of the first years deposition 
after T years of consolidation. With sediment continuing to enter the reser- 
voir, the average density of all sediment deposited would be the average of 
all the densities from one to T years of consolidation. Miller integrated the 
equation from one to T years and developed the equation shown at the bottom 
of Table 5. Waye is the approximate density for the sediment deposits after 
T years of compaction. 

In the first method by Hembree, Colby, Swenson, and Davis, median parti- 
cle sizes for each sample were plotted against the suspended sediment load. 
A curve was drawn through all points on the curve so that about 50 percent 
by weight of the particles in each size range would be finer than the median 
size indicated on the curve. Median particle sizes were taken for predeter - 
mined size ranges and combined with their values of density given in Table 2 
to obtain an average density for the entire sample. 

The second method of Hembree, Colby, Swenson, and Davis was to divide 
the sediment into size fractions and compute the density from the percentage 
in each size fraction and the density of each size. To determine the average 
percentage of sand, silt, and clay in the suspended sidement, correlation 
curves between percentages of each class and suspended sediment discharge 
were prepared. Frequency tables of suspended sediment discharge were ap- 
plied to the correlated values to determine the tonnage in each classification 
and the average percentage of each classification thus computed. 


CONCLUSIONS 


It is concluded that adequate data are available for determination of densi- 
ties of coarse sediments, in which consolidation is not a problem. However, 
it is the opinion of the writers that the problem of estimating sediment densi- 
ties for finer sediments has not yet been satisfactorily solved. The Lane- 
Koelzer method, with the modification by Miller, is a rational approach, but 
does not appear to utilize all the techniques of the soil mechanics engineer 
that are applicable. It appears that the “always submerged” condition could 
be readily solved by application of the consolidation theories of solid me- 
chanics. While this does not represent the only type of sediment deposition, 
it is certainly one of the more important. The effect of drying on exposed or 


partially submerged deposits needs further study before such techniques can 
be applied. 
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III) by J. N. Bradley and A. J. Peterka. (Proc. Paper 1403, Octo- 
ber, 1957. Prior discussion: none. Discussion closed.) 


The Hydraulic Design of Stilling Basins: Stilling Basin with Slop- 
ing Apron (Basin V) by J. N. Bradley and A. J. Peterka. (Proc. 
Paper 1405, October, 1957. Prior discussion: none. Discussion 
closed.) 


The Hydraulic Design of Stilling Basins, by J. N. Bradley and A. 

J. Peterka. (Proc. Papers 1401-1406, October, 1957. Prior dis- 

cussion: none. Discussion closed.) 


The Hydraulic Design of Stilling Basins: Hydraulic Jumps on a 
Horizontal Basin (Basin I) and Stilling Basin with Sloping Apron 

(Basin V) by J. N. Bradley and A. J. Peterka. (Proc. Papers 1401 

and 1405, October, 1957. Prior discussion: none. Discussion 

closed.) 
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The Hydraulic Design of Stilling Basins: Small Basins for Pipe 
or Open Channel Outlets—No Tail Water Required (Basin VI) by 
J. N. Bradley and A. J. Peterka. (Proc. Paper 1406, October, 
1957. Prior discussion: none. Discussion closed.) 

by John R. Argue..... 


Discharge Characteristics of Thin-Plate Weirs, by Carl E. 
Kindsvater and Rolland W. Carter. (Proc. Paper 1453, December, 


1957. Prior discussion: none. Discussion open until May 1, 1958. 
by Turgut Sarpkaya 
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OCEAN WAVE FORCES ON CIRCULAR CYLINDRICAL PILES* 


Closure by R. L. Wiegel, K. E. Beebe, and James Moon 


R. L. WIEGEL,! K. E. BEEBE,” and JAMES MOON. 3—Anyone who has 
dealt with large ocean waves would not be surprised at a large amount of 
scatter in measurements of phenomena associated with these waves. No two 
ocean waves are identical, they only appear to be so due to the rather crude 
methods used to describe wave characteristics, considering that they are 
non-uniform in three dimensions. This is especially true when both swell 
and locally generated wind waves are present. 

Mr. Stelson makes the statement that there was an area of deficient under- 
standing regarding inertia drag. There is much that is not known on inertia 
drag, especially when associated with large velocities and the resulting phe- 
nomena; however this is not what he refers to. Mr. Stelson has either mis- 
understood the definition of the term Cm, which he has called the added mass 
constant—which it isn’t—, or he has misunderstood the problem in general. 
We are dealing with a fluid which accelerates about a cylinder at rest, not 
with a cylinder accelerating through a fluid at rest—which is what he dis- 
cusses. The force necessary to accelerate a mass through a fluid is associ- 
ated with the mass itself and what is termed the added mass (for a right cir- 
cular cylinder the added mass constant is unity). The inertia force exerted 
by a fluid accelerating about a submerged body is associated with the mass 
of the displaced fluid and the added mass; the coefficient of mass as given by 
Lamb for this case is Cy = 1+ 1 = 2 for a right circular cylinder.* This 
development of the coefficient of mass was used because of the ease in under- 
standing the physical phenomenon. Of far greater importance is the fact that 
MacCamy and Fuchs(1) have solved the diffraction problem for a right circu- 
lar cylinder with a diameter small compared with the wave length, using the 


linear theory for waves in shallow water. When the resulting equation is put 
in the writers’ terms it is 


which shows a value of 2 for what has been called the coefficient of mass. 


. Proc. Paper 1199, April, 1957, by R. L. Wiegel, K. E. Beebe, and James 
Moon. 


1. Associate Research Engr., Inst. of Eng. Research, Univ. of California, 
Berkeley, Calif. 

2. Graduate Research Engr., Inst. of Eng. Research, Univ. of California, 
Berkeley, Calif. 


3. Cons. Engr., Signal Oil & Gas Co., Los Angeles, Calif. 
See Reference 14 of original paper. 
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The scatter cannot be accounted for any one of the reasons stated by the 
writers; however, each of these reasons accounts for a portion of the scatter. 
Undoubtedly there are also other reasons. 

Some advances are being made in understanding the basic phenomenon. 
For example, data in the paper presented by John McNown and G. H. 
Keulegan, “Vortex Shedding and Resistance in Unsteady Flow” at the October 
1957 meeting of ASCE showed a strong dependency of the coefficients of mass 
upon the wake for a right circular cylinder in standing water gravity waves. 
Further, tests made by A. D. K. Laird at the University of California using a 
right circular cylinder pair extending through an air-water free surface as 
the mass of a long pendulum oscillating in a direction approximately parallel 
to the free surface showed that variations in forces as great as 50 to 100 
percent could occur when large eddies in the wake of one cylinder hit the 
other cylinder in a certain manner. 


REFERENCE 


1. MacCamy, R. C. and R. A. Fuchs, Wave forces on piles: a diffraction 


theory, U. S. Army, Corps of Engineers, Beach Erosion Board, Tech. 
Memo, No. 69, 17 pp. December 1954. 


ASCE 1616-7 


THE ESTIMATION OF THE FREQUENCY OF RARE FLOODS* 


- Closure by Benjamin A. Whisler and Charles J. Smith 


; BENJAMIN A. WHISLER,! M. ASCE, and CHARLES J. SMITH,” A.M. 
ASCE.—The writers welcome the comments of Messrs. Beard, Benson, and 
Williams in their discussions of this paper. Several questions have been 
raised which the writers feel will be answered only after more data become 
available and after the suggested procedure and other procedures have been 
applied to the data and the results compared. Since there appears to be no 

general approach to the problem of estimating the frequency of rare floods 

Be that proves to be satisfactory under all conditions, the writers believe that 

some completely new thinking, at present unknown to them, is needed. In = 

time, with enough people working on the problem, a satisfactory approach 

may be found. 


a. Proc. Paper 1200, April, 1957, by Benjamin A. Whisler and Charles J. 
Smith. 

1. Head, Dept. of Civ. Eng., Pennsylvania State Univ., Univ. Park, Pa. 

2. Associate Prof., Dept. of Civ. Eng., Pennsylvania State Univ., Univ. Park, 
Pa. 
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SYSTEMATIC CHANGES IN BEDS OF ALLUVIAL RIVERS 


Discussion by Vito A. Vanoni and Hsin-Kuan Liu 


VITO A. VANONI,! M. ASCE.—The authors are to be commended for 
bringing their important observations of bed configuration of the Mississippi 
River to the attention of the profession. Such detailed observations, tedious 
and time-consuming as they may be, are necessary to the advance of our un- 
derstanding of the mechanics of alluvial streams. In this regard, the writer 
not only agrees with the authors conclusion(1) that the study of river bed be- 
havior might throw new light on some theoretical aspects of river hydraulics, 
but he feels certain that it will do so. 

The authors’ belief that sand waves in alluvial rivers form a major source 
of resistance to the flow is borne out by laboratory studies, as shown by the 
data(A) in the following table. These data were obtained from experiments in 
which the bed of the flume was covered with loose sand with a geometric 
mean grain size of 0.152 mm and geometric standard deviation of sizes of 
1.76. The flow rate was set for each experiment and the slope was varied to 
obtain uniform flow. It will be seen that the friction factors with dune- 
covered beds are from two to five times as much as those for the flat-bed 
condition when no dunes exist on the bed. It will also be noted that as the 
velocity and the sediment-discharge concentration increase the friction factor 
decreases. Detailed observations of all dunes were not made in these experi- 
ments; however pictures and visual observations did not show any great 
change in size or shape of dunes until the sand wave and smooth phases were 
approached, i.e., runs 6, 7 and 8. 

The behavior of dunes in the flume seems to be opposite to that observed 
by the authors in the Mississippi River where dunes increase in size with 
stage. However, in the flume the increase in velocity was achieved by in- 
creasing the pump speed and changing the slope, while in the natural river 
the slope remains essentially constant as the stage varies. Therefore, the 
two cases, i.e., river and flume, are not exactly the same. Any quantitative 
data that the authors have on the size of waves at high stages would be of 
considerable interest. 

It would be interesting to obtain simultaneous measurements of friction 
factor and dune configuration for the river, and it is hoped that the authors 
will have an opportunity to make such observations. It is also hoped that such 
observations can be extended to the highest flood stages. In the absence of 
such data one can only speculate on the effect of dune size and configuration 
on the friction factor. One is inclined to suppose that the friction factor 


a. Proc. Paper 1331, August, 1957, by Walter C. Carey and M. Dean Keller. 
1. Prof. of Hydraulics, California Inst. of Technology, Pasadena, Calif. 
* Letters refer to bibliography at end of discussion. 
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would increase as the dune size increases. Since the highest dunes were ob- 
served at the highest stages, this would mean that the friction factor in- 
creases as the stage increases. However, measurements of the Manning 
roughness coefficient n, by Eden(B) for the Mississippi River in the vicinity 
of Memphis, Tennessee, showed that as the stage increased beyond about 30 
feet, n stayed approximately constant or decreased slightly. As the stage 
dropped below about 30 feet, Eden found that n increased rather rapidly and 
attained the highest value at the lowest stages. Eden also observed that the 
largest dunes occurred at the highest stages. He observed dunes as high as 
34 feet, which agrees well with the maximum height reported by the authors 
for stages well below the flood stage. 

A careful observation of the profiles presented by the authors seems to 
show that more of the bottom was covered by small dunes at low stages than 
at high stages. It is possible that the small dunes actually offer more re- 
sistance to flow than the larger waves, which seem to have relatively smooth- 
er surfaces than the configurations of smaller height which occur at the lower 
stages. 

The authors explain the loop in the rating curve in part by postulating a 
lag in adjustment of the bed configuration. It seems unlikely to the writer 
that such a lag exists since the rise and fall of the flood occur slowly, re- 
quiring days and weeks to take place. During this time the stages are rela- 
tively high, as are the velocity and sediment transporting capacity of the 
stream, and it seems to the writer that the stream should be able to wash 
away the larger dunes and readjust the bed with very little lag. The fact that 
the authors find that crossings adjust very rapidly lends support to the writ- 
er’s belief that the lag is small. 

By following the idea advanced by Brooks,(C) that both flow rate and sedi- 
ment transport rate are independent variables, one can explain the loop in 
the rating curve without resort to the idea of lag advanced by the authors. 

In line with the idea proposed by Brooks, one visualizes the dependent vari- 
ables, such as velocity and friction factor, as adjusting themselves to the in- 
dependent variables of flow rate and sediment transport rate. For instance, 
let us compare conditions in the stream for the same flow rate during the 
rising and falling stages of the flood. On the rising stage the sediment load 
is usually high,(D) so the stream must adjust itself to transport this load. 
This is done by smoothing out the bed and achieving a relatively high velocity. 
On the falling stage, the sediment load is usually smaller than on the rising 
stage, so the stream can transport the load with a lower velocity and larger 
depth. The stream then adjusts to this condition by forming rougher dunes 
than it did on the rising stage. 


BIBLIOGRAPHY 


A. Vanoni, Vito A. and Brooks, Norman H., “Laboratory Studies of the Rough- 
ness and Suspended Load of Alluvial Streams,” California Institute of 
Technology, Sedimentation Laboratory Report No. E-68, December 1957. 


B. Eden, Edwin W., “A Study of Bed Movement and Hydraulic Roughness in 


the Lower Mississippi River,” Masters’ Degree Thesis, University of 
Iowa, June 1938. 


C. Brooks, Norman H., “Mechanics of Streams with Movable Beds of Fine 
Sand,” Proceedings, ASCE, April 1955, Vol. 81, Separate No. 668. 


_| 


ASCE DISCUSSION 1616-11 


D. Leopold, Luna B. and Maddock, Thomas, Jr., “Relation of Suspended- 
Sediment Concentration to Channel Scour and Fill,” Proceedings of the 
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Darcy-Weisbach Friction Factor f, for Sediment-Laden Flows 


0.241 ft Deep ina Flume 10.5 in. Wide 


U Cc f 
Run Average Sediment Friction Bed 
No. Velocity Discharge Factor Condition 
fps Conc. 


gr/t 


Dunes 
Dunes 


Dunes 


Dunes 
Dunes 
Dunes 
Dunes 
Sand wave 
Sand wave 
Flat 

Flat 

Flat 


0. 
0 
0 
2 
3 
3 
3 
5 


HSIN-KUAN LIU* A.M. ASCE.—The authors are to be congratulated for 
presenting information on the sand waves of a natural river. Such informa- 
tion is highly valuable to the profession. 

For several years, the writer has studied sand-wave motion under labora- 
tory conditions. The data collected from laboratory experiments are general- 
ly for steady, uniform flow in narrow straight channels. It is obvious that 
this type of flow is different from what has been observed by the authors 
from the Mississippi River. Ordinarily the flow phenomenon of natural 
rivers is so complex that it is not possible to understand the hydraulics of 
rivers based upon field investigations alone. On the other hand, conclusions 
based upon investigation under simplified laboratory flow conditions should 
not be applied to prototype flow conditions without due consideration of the 
complex problems involved in natural rivers. In the laboratory it is easier 
to study analytically the primary forces which affect the movement of 


* Asst. Prof. of Civ. Eng., Dept., Colorado State Univ., Fort Collins, Colo. 


A 
7 
2b 0. 80 30 124 
2a 0.85 59 115 a 
2 0.91 82 115 | 
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3 1.04 103 
4 1.20 075 
5 1.38 050 } 
6 1.55 038 
7 1.69 . O29 
8 1. 83 . . 024 
1 2.06 
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sediment. As a result, various concepts of mechanics of sediment transport 
have been formulated. It is known that these primary forces such as tractive 
force, lifting force, turbulent stress, vary with the flow conditions. The flow 
conditions in natural rivers are somewhat different from those in the labora- 
tories. For example, in natural rivers the flow is non-uniform and unsteady, 
the scale of turbulence is large, and the geometric scale of a river is large; 

these are not so in the labroatory. Despite the fact of these differences, the 

authors’ conclusion based upon field investigations are shared to a large ex- 

tent by the writer who has been acquainted only with laboratory conditions. 

The authors are correct in stating that sediment transport in the bend re- 
gion is different from that of a straight reach. In the past, most of the labora- 
tory study of sediment transport has been conducted in straight channels. 
Laboratory study of sediment transport in curved channels has not been done 
extensively. In view of the fact that river bends play an important role in the 
river hydraulics, it is obvious that studies of sediment transport around the 
bend should be emphasized both in the field and in the laboratory. 

The statement that sand waves occur generally on the bed of a river ina 
systematic manner agrees with the laboratory data. This fact does not sup- 
port the hypothesis that sand waves and ripples are caused by turbulence as 
suggested by some hydraulic engineers. Needless to say that turbulence 
inter-reacts with the sand waves which in themselves cause turbulence. 

That the size of sand waves becomes larger with increases in discharge 
and smaller with decreases in discharge may be owing partly to the effect of 
change of depth, partly to the effect of change of velocity.(a) Studies(b) made 
by W. H. Powless indicated that the sizes of the sand waves are also affected 
by the rate of rise and fall of the stage. 

The writer agrees with the authors that ‘the major sources of resistance 
arises from systematic river-bed changes, which the river makes of its own 
volition, and modifies to meet changing discharge requirements.” In other 
words, the bed roughness may have correlation with the bed configuration. 
However, the bed configuration is not an independent variable, but it is de- 
pendent upon the river itself. 

No information regarding the rate of the sand-wave motion was given by 
the authors in their paper. According to laboratory data, the writer found 
that 


N 


C = KV 


(a) 


in which C is the velocity of the sand waves in ft per sec. 
V is the mean velocity of the flow in ft per sec 
N is the pure number 


K is an empirical coefficient, probably is a function of the flow and 
the bed material. 


According to laboratory data, the N-value is 5 and K is in the order of M X 
10-4, in which M varies from 1 to 10. According to data collected by Wolman 
on the Du Noir Creek, Dubois, Wyoming, equation (a) still holds true as 
shown in Fig. A and N-value is also 5. For natural rivers there is very 
limited information regarding the rate of sand-wave motion. Table I was 


prepared by the writer by using data on the Mississippi at Fulton as reported 
by the Mississippi River Commission. (?) 
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DISCUSSION 
Table I. Rate of Travel of Sand Waves 


rate of travel 
Gage ft/day ft/sec 
ft x 10-4 


22 to 24 rising 
18 to 22 rising 
13 to 18 rising 
16 to 15 falling 
15 to 13 falling 
13 to 12 falling 


Mean velocity 
Date of Surveys 


1880 
July 7 to July 13 
July 3 to July 7 
June 23 to July 3 
June 11 to June 15 
June 15 to June 19 
June 19 to June 23 


Data in Table I are plotted on Fig. A. Those of falling stage are plotted 
separately from those of rising stage. Data are insufficient for drawing 
definite conclusions. However, there are several points that can be dis- 
cussed. 


a) The K-value for the sand waves found in the Mississippi River is much 
smaller than that of laboratory data. 


b) There are not enough data to determine the N-value. 


c) The mean velocity on a rising stage needs to be higher than that on a 
falling stage for a constant rate of sand wave motion. Regarding items c, W. 
H. Powless stated in his report(b) that the rate of sand wave motion is much 
greater on a rising than on a falling stage. Additional information of sand 
wave motion in Mississippi River is reported by E. W. Lane(C) and E. W. 
Eden. (qd) 

The Mississippi River data seem to indicate that the sand waves are in- 
fluenced by the rising and falling stages of the flows. In order to understand 
this phenomenon, a laboratory study of sediment transport by introducing a 
hydrograph into the flow system will be very helpful. 

The writer had the opportunity to discuss the authors’ paper with Mr. E. 
W. Lane. Mr. Lane thought that the authors’ information is highly valuable, 
however, he pointed out that most of the data presented by the authors were 
collected in a reach where the slope of the Mississippi River is appreciably 
affected by the sea level. Furthermore, the Mississippi River in this reach 
is deep and narrow. In order to understand fully the sand wave motion in the 
Mississippi River data should be collected also in other reaches where the 
flow conditions are different. 

In conclusion it can be stated that, alluvial rivers, such as the Mississippi, 
carry not only large amount of water but also large amount of sediment. In 
order to understand the nature of these rivers, investigation and study have 
to be made not only of the motion of the water, but also of the motion of the 
sediment carried by the river. To a certain extent the study of the surface 
water is easier because it can be seen from above the water surface, and ex- 
tensive studies have been made in an effort to understand its hydraulics. The 
study of sediment transported by an alluvial river is obscured by the fact that 
the movement is under water. For this reason, special effort has to be made 
in order to understand the behavior of alluvial streams. Such understanding 
is essential in order to control and develop alluvial rivers for the welfare of 
mankind. 

The writer wishes to express his appreciation toward Messrs. M. G. 
Wolmanand E. W. Lane for their help in the collection of data and of refer- 
ences, respectively. 
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22 2.54 5.2 

18 2.08 4.6 
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22 2.54 4.2 

16 1.85 4.1 
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TURBULENCE IN A DIFFUSOR BOUNDARY LAYER®# 
Discussion by Theodor S. Strelkoff 


THEODOR S. STRELKOFF.*—Before the problem of turbulent boundary 
layers can be solved completely, it seems apparent that a wealth of data will 
be required, data which relate turbulence characteristics to mean flow 
parameters, and which are obtained from flows confined by various geomet- 
rical configurations. This fact makes the present paper most welcome. 
With reference to the terms in the turbulent energy balance, the produc- 
tion term, when describing conditions in flow with a pressure gradient, es- 


pecially anywhere near separation, should include in addition to wer at 
least u' a variable which may be half the size of the shear stress term. 


Although another applicable term, v' x » partially nullifies the effect of the 
shear stress. (1) 

(The expression for W, as given in the paper apparently suffered from a 
misprint as it was being non-dimensionalized into its final form, because in- 
stead of u4 one reads u4.) 

It would be of some interest to know what method was used in placing “the 
plane of the hot wire at the wall of the diffusor.” As it is generally difficult 
to locate the probe accurately, any information on this point would be helpful 
to other researchers. In addition, one would wish to know what corrections 
if any were applied near the wall to compensate for heat losses to the wall 
itself. 

If the true three-dimensional momentum thickness had been plotted, com- 
parison with Uram’s data(2) would still have been convenient, and in addition 
the Karmd4n equation, modified to include normal stress terms, could have 
been employed for calculation of the shear stress at the wall to be compared 
with the value obtained at station 12 with the pitot tube method. 

It would also have been helpful if a description of the “history theory,” 
referred to in Fig. 6, had been supplied, together with indications of how the 
calculations were made and what the various symbols stood for. Schubauer 
and Klebanoff,(1) basing their reasoning on a simplified form of the Reynolds 


equation, which neglects viscous stresses, derived the relation op =9T 


Ox Oy’ 
where T is the turbulent shear stress. It is interesting that at station 12 a= 
can be computed as 2.2 psf/ft, in fair agreement with oP = 2.56 psf/ft. The 


a. Proc. Paper 1393, October, 1957, by J. M. Robertson and G. L. Calehuff. 
* Research Associate, Iowa Inst. of Hyd. Res., State Univ. of Iowa, Iowa 
City, lowa. 
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value of 2.2 psf/ft was obtained from the stated value of Rg = 25, 000, by as- 
suming an air temperature of about 80° and taking the slope of the shear 
stress coefficient function at y = 0 as 0.016 per inch. Actually, this computa- 
tion should be carried out with variables whose values are directly known, as 
an error of 20° in air temperature, for example, would cause a 10% error in 
OT 

If the graphs of u' had been normalized through division by the velocity U; 
instead of U, the peak, noted by other researchers studying adverse pres- 
sure gradient boundary layers, would have appeared at each station. The 
data of Schubauer and Klebanoff(1) peaked in this manner; that of Sandborn 
and Slogar(3) exhibited a somewhat similar peak at the station farthest down- 
stream. Both sets of data applied to two-dimensional boundary layers. 

It would have been informative to see a graph similar to figure 8 with data 
from all stations superimposed as an indication of the applicability of the 
inner-law to turbulent fluctuations. 

Further studies with the same apparatus seem to be definitely in order as 
time permits. In particular, careful investigation of U, u', and Uv very near 
the wall should prove profitable. 
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100 FREQUENCY CURVES OF NORTH AMERICAN RIVERS 


Discussion by G. N. Alexander, Max A. Kohler and Ralph W. Powell 


G. N. ALEXANDER* A.M. ASCE.—The first conclusion made by the author 
in his synopsis reads - 


“1. That frequency curves can generally be drawn as straight lines.” In 
the light of the paper this statement is interpreted as meaning that the most 
suitable distribution for flood flows is the log-normal. The author carried 
out a sampling (Monte Carlo) experiment where he assumed a log-normal 
population of 1000 floods and drew twenty random samples of size fifty. He 
noted that two of them were far from linear and rightly states that this “is no 
conclusive proof that frequency lines have to be straight.” In fact it is just 
the opposite. From this experiment it is clear that samples, even if drawn 
from the ‘true’ population may be far from linear. What is more, even if a 
sample is linear it may still be a poor representation of the population. Even 
in a record of 1000 years the order statistics will show a pattern of scatter- 
ing similar to a small sample. Presumably, one would obtain an even greater 
scattering than that found by the author if the samples of fifty were drawn 
from a ‘sample’ of 1000 from a very large population, rather than the ‘true’ 
population of 1000. 

Statisticians, and engineers familiar with quality control, have recognized 
these sampling effects, but many hydrologists are still not sufficiently aware 
of the effects of sampling and hence apt to ascribe them to other causes. The 
problem is to differentiate between sampling effects and characteristic ef- 
fects, and the writer believes that is best done if we first have some concept 
of the effects of sampling. 

These effects can be shown by using a ‘band’ instead of a line to represent 
the distribution. The width of the band would approximate to the interval be- 
tween the outer pair of samples obtained by the author from his 20 samples 
if we wish to cover about 95% of the samples from his population of 1000. 
Hence if we were to use a 5% significance test one would ascribe fluctua- 
tions within the band to sampling effects rather than seek natural causes. 

The fourth conclusion in the author’s synopsis reads - 


“4. That an increase in the size of the drainage area tends to flatten the 
frequency line.” 

The author has tried to prove this by an a priori argument comparing the 
slope of the frequency curve of annual floods obtained for stations a or b on 
tributaries identical both physiographically and climatologically with that for 
c just below their confluence. In the standardized log-normal distribution as 


a. Proc. Paper 1395, October, 1957, by E. Kuiper. 
* Hydr. Engr., U. S. Geological Survey, Washington, D. C. 
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plotted by the author, the slope may be measured by the ratio of the standard 
deviation o to the mean /, i.e. the coefficient of variation C.V. 

Now under the above assumptions the mean flood at c, mu (c) say, will be 
twice the mean flood at a, uw (a) or 


H (c) = 


2 uw (a). (1) 
If in the first instance we assume the floods at c are a random conbination of 
the annual floods at a and b then for uncorrelated distributions, the variances 
(square of the standard deviations) may be taken as additive (this ignores 

their non-normality) 


Var (c) = 2 Var (a) 


(2) 
or o (c) = oa (a) 


or from (1) & (2) 


Hence, in this case slope of c is flatter than slope of a. 
However, since we have tributary streams there will be some correlation 


P say between their flows and the relation between a, b, & c may be written 
as 


Var(c) = Var(a) + Var(b) + 2p.[Var(a).Var(b)]® 


For identical tributaries 


Var(a) = Var(b), 


Var(c) = 2(1+/P). Var(a). 


Hence, using (1) and (3), 


C.V.(c) 


= 4 /2UIFp) 
C.V. (a) 


giving as before s if P=0. If P = 1., i.e. perfect correlation the slopes are 


identical. 
However, the floods at c will not generally be a combination of the “annual” 
floods at a and b as assumed above. We could divide each annual flood at c 
into its a ‘and b « components” and hence obtain the distributions of floods in 
aandb causing the annual floods at c. The distributions of these components 
at a will usually have a wider range and hence greater variance than the an- 
nual floods at a, and will also have a lower mean. Hence the “component” 
distribution at a will have greater C.V. than the “annual” distribution at a. 


C.Ve(c) = 
C.V.(a) 
Hence 
(3) 
| 
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Thus although the combination of the components will result in annual floods 
at c with a lesser C.V. than the component floods in a or b, this C.V. may or 
may not be less than the C.V. of the annual floods at a or b. 

Actual analyses are therefore, required to justify the general thesis that 
the slope decreases with the drainage area: Such analyses now being under- 
taken by the U.S. Geological Survey do support this thesis. 

The author raises the question of suitable parameters—stating “it would 
be desirable if the inclination of the frequency line could be expressed in 
some numerical parameter.” He has selected the ratio of the 1% flood to the 
mean annual flood. In dealing with the log-normal distribution a decision has 
to be made concerning the stage, if any, in the analysis at which the original 
data will be transformed to log-units. The statistical analysis is consider- 
ably simplified and significance tests depending on normality may be used 
more readily if the transformation to log-units is made at the outset. In this 
case it is logical to use the independent and efficient parameters, namely, 
the mean Y and variance o2 of the common logarithms as the parameters to 
describe the distribution. It is readily shown (1, p. 161) that in a “true” log- 
normal distribution the mean of the logarithms is also the geometric mean 
of the original variate. This value is also the 50% or median value, with a 
mean return period of two years, irrespective of whether we use the loga- 
rithms or the original variate. It therefore seems more logical to measure 
“central tendency” by the median or geometric mean rather than by the 
arithmetic mean of the original observations used by the author. The rela- 
tion between a, Y & o is given by (Ref. 1, p. 161) 


As a measure of ‘dispersion’ the standard deviation, in some form or 
other, is almost universally used. For some curious reason, statisticians 
who derided as old fashioned the engineer’s use of the “probable error” with 
its 50% interval were apparently quite happy to measure in percentage sig- 
nificance levels, rather than using the natural unit, viz.; the standard devia- 
tion. As probability paper is not normally printed to facilitate measurement 
in standard deviations, and statistical tables are usually constructed with 
significant levels shown percentage wise, it is generally more convenient to 
use percentages as the author has done. 
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MAX A. KOHLER,* M. ASCE.—The study reported by Mr. Kuiper was di- 
rected toward the development of a technique for using discharge records at 
numerous points in an area when developing the frequency curve for a short- 
record station. This generalized approach to frequency analysis may be ex- 
pected to yield improved results (over single-station analysis), and is perhaps 
to be recommended when the nature of the problem does not justify detailed 
study of rainfall-runoff relations,(1) etc. The approach is in some respects 
similar to that widely employed by the U. S. Geological Survey,(2) but 
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Mr. Kuiper attempts to introduce “climate” as a factor, whereas the Survey 
restricts each regional study to an area of climatic homogeneity. 

Although much is said of climate in the paper, the discussion under “Pos- 
sible Application” would indicate that the frequency curve is derived from 
only two independent factors—mean-annual flood and drainage area. Dividing 
the 100-year flood by the mean-annual and the mean-annual flood by the 
drainage area only tends to obscure the high correlation which exists in the 
basic data. 

The upper plotting of the accompanying chart (derived from data in Table 
7) shows that the 100-year and the mean-annual floods are highly correlated 
and that the residual error is not appreciably correlated with drainage area. 
Although not immediately evident from the plotting, study of the data shows 
the relation is largely independent of physiography or climate. On the other 
hand, the plotting in the lower portion of the accompanying chart shows clear- 
ly that the mean-annual flood is highly correlated with both drainage area and 
physiography (or climate). 

Treating the entire 100 stations as constituting a single population may 
enhance the reliability of the relation between the 100-year and mean-annual 
floods, but the data must be classified by basin, or region, in deriving the re- 
lation between drainage area and mean-annual flood. At least this is the case 
until such time as we are able to explain the basin-by-basin variation in terms 
of measurable factors. Even so, the mean basin-curves in the lower part of 


the chart present an interesting array. Without belaboring the issue, it can 
be said that: 


1. As would be expected, the mean-annual flood increases with increasing 
annual precipitation and runoff (compare Kansas and Ohio River Basins). 


2. Proceeding downstream, the rate of change of mean-annual flood with 
drainage area changes abruptly when a major tributary of different 
climatic and physiographic characteristics enters (see points 27, 42 
and 74, which are for Keokuk, Kansas City and St. Louis, respectively). 
In other words, the relation between drainage area and mean-annual 
flood for mainstem stations may differ appreciably from that for sta- 
tions on adjoining tributaries. 


3. Climate and physiographic characteristics affect the exponent as well 
as the coefficient in the equation relating area and mean flood. 


REFERENCES 


1. Paulhus, Joseph L. H. and John F. Miller, “Flood Frequencies Derived 
from Rainfall Data,” Journal of the Hydraulics Division of ASCE, Vol. 83, 
No. HY 6, December 1957. 


2. Dalrymple, Tate, “Regional Flood Frequency,” Highway Research Board, 
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RALPH W. POWELL, |! M. ASCE.—Mr. Kuiper is to be thanked for as- 
sembling this large mass of data and giving the profession the benefit of his 
one hundred frequency diagrams. However, his treatment of the data raises 


1. Visiting Prof., Applied Mechanics, Univ. of Kansas, Lawrence, Kansas. 
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several questions. In the first place, would it not have been better to use the 
median annual flood, rather than the mean annual flood for the denominator 
of the ratio represented by R? It would seem that if the sample is unbal- 
anced by having too large or too small a number of large floods, it would 
produce more error in the mean than in the median. 

The second question regards using the eight-tenths power of the drainage 
area as denominator of the ratio represented by F. This evidently came 
from Fuller’s formula, but it is doubtful if it applies here. In fact, by 
Fuller’s formula all the R’s would be 1 + 0.8 log 100 = 2.6. The writer 
plotted F against R on semi-logarithmic paper and drew a regression line by 
eye with the equation log F = 2.5 - 0.5 R. There was a great deal of scatter 
and the larger drainage basins fell below the line and the smaller ones above. 
Stations with approximately the same R were then grouped together and for 
each group the mean annual flood plotted against drainage area. It was found 
that the mean annual flood tended to vary as about the two-thirds power of the 
drainage area. Values of K' = Q/A2/3 were then tabulated, where Q = mean 
annual flood in c.f.s. and A = drainage area in Sq. mi. This is easily done by 
using tables of two-thirds and three-halves powers from a hydraulics hand- 
book. (Moving the decimal point three or six places in A moves it two or four 
places in A2/3). F' was then plotted against R and the regression line drawn 
by eye has the equation log F' = 3.30 - 0.55 R. The scatter was reduced only 
a little, but there no longer seemed to be any systematic variation with drain- 
age area. As mentioned above, it is possible that using the median flood 
rather than the mean might reduce the scatter somewhat. However, it is 
quite evident that there are other factors determining R besides the relation 
of the mean or median annual flood to the drainage area. 

As the author points out, it is probable that one of these is the climate as 
represented by the amount of the mean annual rainfall. Perhaps some mea- 
sure of the variability of the rainfall would be a still better criterion. An- 
other factor mentioned by the author is the average permeability of the soil. 
But to try to develop a formula involving these and any other factors would 
require a vast amount of research. It raises the question as to whether it 
might not be just as feasible to try to develop a formula for estimating the 
one per cent flood directly. Attention is called to the fact that Pettis(1) pro- 
posed such a formula many years ago. It would be interesting if some one 
had the time to test this formula against the data the author has assembled. 
More recently Potter(2) proposed a formula for the ten year (ten per cent) 
flood for small drainage areas. Perhaps some day some one can investigate 


whether it would apply to larger drainage areas, or how it would have to be 
modified to do so. 
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THE HYDRAULIC DESIGN OF STILLING BASINS: 
HYDRAULIC JUMPS ON A HORIZONTAL APRON (BASIN 1)@ 


Discussion by E. A. Elevatorski and A. Rylands Thomas 


E. A. ELEVATORSKI,* J.M. ASCE.—The authors have presented useful 
information on the elements of the hydraulic jump formed in a rectangular 
channel on a horizontal apron. To the designer, the height and length of the 
hydraulic jump are of particular significance in determining the necessary 
dimensions of a stilling basin or channel transition. 

D 
In Figure 5, the authors’ relationship of > versus Fy verifies, without 
1 
doubt, the applicability of the pressure-momentum formula to the hydraulic 
jump. However, the extrapolation of the lower curve seems to be in error, 
since the undular-type jump could not be formed unless the Froude number 
exceeds 2, regardless of the tailwater depth. 

As an alternative to the employment of Figure 5 for finding the downstream 
depth, Dg; a dimensionless relationship between Do, Dj, and the critical 
depth, De can be found. This seems to have merit since the Froude number 
need not be computed for each discharge. The following analysis presents 

D, Do 
the relationship between D and D for a given discharge. 
c c 

For a rectangular channel, the difference in hydrostatic pressures at sec- 
tions (1) and (2), referring to the beginning and end of the jump, respectively, 
must be equal to the change of momentum per second of the mass flow. 


Equating and transforming, Dj and Dg in terms of q, the discharge per foot of 
channel width becomes 


D, (D, + D,) = (1) 


Since dD.” = q” for a rectangular channel, equation (1) becomes 


3 
D, D, (D,+D,) = 2D, (2) 


To obtain dimensionless ratios, divide eath term of equation (2) by D,°. 


- Proc. Paper 1401, October, 1957, by J. N. Bradley and A. J. Peterka. 
Asst. Prof. of Civ. Eng., Univ. of Arizona, Tucson, Ariz. 
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Therefore | Do 2 (3) 


D 


To assist the designer, Figure 1 presents the relationship between > 
D 
and D- for the range normally encountered by the designer. For a given dis- Y 
c 


charge q(c.f.s/ft. of basin width), the critical depth can be computed from 
equation (4) 


The longitudinal element of the jump is, without doubt, the most difficult 
element to measure. Certainly this is partially due to the differences in 
opinion of exactly where the terminus of the jump lies. The writer prefers 
the following definition: 


“The origin of the jump is the point where the roller turbulence begins 
and the water becomes white and foamy, due to the large amount of en- 
trained air. Downstream, the terminus of the jump is defined as the 
point where no return flow is observable.” 


This seems to be in harmony with the definition of the authors who judged 
the length of the jump as it would best represent the terminus of a conven- 
tional stilling basin. 

Contrary to the conclusions of the authors, the writer is of the opinion that 
the length of the jump is independent of the Froude number. Past studies 
clearly show that the length of the jump, Lj, bears a definite relation to the 
jump height, Dj- To support this, the writer has taken the liberty of using the 
authors’ published data to plot Lj versus Dj = Dg - Dy in Figure 2. Despite 
the fact that the Froude number varied from 2 to 20, a good correlation was 
obtained. From the plot, the correlation can be summarized by equation (5). 


L; = 6.6 (Dy - D,) (5) 


Equation (5) compares favorably with the results obtained by other investi- 
gators and listed below: 


Investigator 


Length of Jump Formula 


Smetana L, =6 (Dy - D,) 

wu? L; = 6.8 (Dy - D,) 
Safranez® L, = §.2 (D, D,) 
Page2 L, = 5.6 D 


j 2 
The writer wishes to emphasize that equations (3) and (5) are only 
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applicable for hydraulic jumps formed in rectangular channels on horizontal 


L 


aprons. For trapezoidal basins and channels, the ratio = is much greater 
j 


than the corresponding ratio for a rectangular channel. Experiments by 
L, 
Hsing(1) indicate that = increases as the side slope of the channel decreases. 
j 
It is strongly recommended that model studies be conducted to ascertain the 


elements of the hydraulic jump in a trapezoidal basin before the particular 
design is adopted. 
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A. RYLANDS THOMAS,* M. ASCE.—The authors deserve great credit for 
their comprehensive study, supported by extensive experimental data, pre- 
sented in the six papers 1401 to 1406. 

The length of the hydraulic jump has been the subject of much discussion 
in the past and, as shown by Fig. 7, considerable differences exist between 
the various series of observations. One reason for this is differences in 
definition of the position of the end of the jump. 

The hydraulic jump is essentially a transition from a high-velocity (super- 
critical) stream to a lesser-velocity (subcritical) stream with a higher sur- 
face level. This involves an upward expansion of flow in which transfer of 
momentum plays a major part, and it is the rate of net upward momentum 
transfer which determines the length of jump. The process of momentum 
transfer is not complete at the end of the surface roller but continues down- 
stream, the rate of upward transfer gradually diminishing as the velocity 
distribution becomes adjusted until the loss of momentum due to bed drag re- 
verses the process into one of net downward momentum transfer. Under such 
conditions of gradual change the difficulty of defining the end of the jump is 
manifest. 


Bakhmeteff and Matzke” conceive the end of the jump as the section at 


. Consultant, London, England. 
b. Bakhmeteff, B. A. and Matzke, A. E., “The Hydraulic Jump in Terms of 
Dynamic Similarity.” Transactions ASCE., Vol. 101, p. 640, 1936. 
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which the surface level reaches its maximum height. This has the drawback 
that it is affected by the surface roughness of the bed and depends on the bed 
slope—a gently sloping bed might increase the apparent length of jump in- 
definitely. The authors (page 14) have taken the jump as ending at a point 
downstream either where the high-velocity jet began to leave the floor or on 
the surface immediately downstream of the roller, though which point was 
used in each case is not stated and the former point must have depended on 
the roughness of the bed. 

The writer feels that rather than assume a length based on arbitrary con- 
ditions, on which experimentors are not unanimous, it is preferable to record 
the length of roller, as done by Kindsvater,© a length which is reasonably 
clearly defined, easily observed and a characteristic parameter of the jump. 

The length of jump may be said to depend on the rate of upward momen- 
tum transfer and on the degree or extent of expansion. These may be com- 
plex functions of Froude Number and depth, that is of Fy and Dy and D2 or Dj 
(height of jump = Dg - Dj). Lj is often expressed as a ratio of Dj and, when 
plotted against Fy, this has the advantage of giving a flatter curve in the range 
of Fy between 2 and 6 than does the ratio Lj/Dg which as may be seen from 
Fig. 7 varies sharply in this range. In fact the Authors’ data when plotted in 
the form Lj/D; : Fy indicates a linear relationship over the whole range, 
Lj/Dj reducing from approximately 7.4 with Fj = 2 to 6.0 with Fy = 16. 

The Authors’ classification of jumps (Figs. 7 and 9) is of much interest 
and is to some extent supported by experience with canal falls in India. These 
consist of a free overfall followed by a stilling basin to contain the hydraulic 
jump, generally with side expansions. Bank scour downstream has been more 
troublesome when the fall was relatively low (low Fj) than when it was high 
(high Fj). In the former case there is less turbulence to stabilize the flow 
and it is possible that jets might oscillate not only in a vertical plane, as 
illustrated by the Authors in Fig. 9B, but also in a horizontal plane, resulting 
in the formation of large eddies. Conditions could no doubt be improved by 
adopting a greater width, resulting in reduced discharge concentration and 
higher Fj, but width was generally determined by other considerations, for 
example degree of contraction required to give a desired stage : discharge 
relation upstream. 

Insufficient evidence has been produced for the conclusion at the foot of 
page 22 that baffle piers are of little value when F, lies between 2.5 and 4.5. 
When Fj is low, baffle piers may be less effective in providing back pressure 
but they can be valuable in this range in stabilizing the flow and reducing the 
length of basin required. 

It would be of Value if tests similar to those described in this paper were 
carried out on jumps in channels with side expansions to determine the basic 


relationships in regard to position and length of jump in laterally expanding 
flow. 


c. Kindsvater, C. E., “The Hydraulic Jump in Sloping Channels,” Transactions 
ASCE., Vol. 109. p. 1112, 1944. 
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HYDRAULIC DESIGN OF STILLING BASINS: 
HIGH DAMS, EARTH DAMS AND LARGE CANAL STRUCTURES (BASIN II) 


Discussion by A. Rylands Thomas 


A. RYLANDS THOMAS, M. ASCE.—The titles of Papers 1402 and 1403 
and the conclusions therein indicate that stilling basins with chute blocks and 
end sill (Type II) are recommended for high dams, earth dams and large 
canal structures and shorter basins with chute blocks and end sill but with 
floor blocks in addition (Type III) are recommended for other canal struc- 
tures, small outlet works and small spillways. In the case of high spillways 
the use of floor blocks might be undesirable because of cavitation but it is not 
clear why in all other cases advantage could not be taken of the benefits of 
floor blocks as evidenced in the Type III design, except in the case of such 
small works that the saving in cost would be immaterial. 

Where cost is the main criterion for the adoption of a Type II design, and 
the choice lies between chute blocks and floor blocks, it may be asked 
whether a single row of floor blocks might not be more effective than chute 
blocks. Comparative effect would be indicated by scour tests in a model with 
mobile bed downstream. 

The observations of jump length are of value but the real criterion of 
minimum safe length of basin is protection against scour downstream. Ex- 
periments have shown that depth of scour is reduced as length of basin is in- 
creased to as much as 12 times the downstream depth.© Such length would in 
most cases be far greater than is required, but the result indicates that 
basin length is a matter of degree of protection required, and this depends to 
some extent on the case considered. Model tests with mobile beds give a 
good indication of scour pattern and comparative depths, but there have been 
cases of deeper scour in the full-scale works than indicated by the models. 

The design shown in Fig. 10 with dentated sill in the middle of the floor 
appears to be an unnatural compromise between baffle and end sill, not likely 
to be very effective as either. To reduce scour at the end of the floor it is 
more important to provide an end sill and the Authors are certainly right fn 
recommending that the sill should be at the end (Fig. 14). The function of a 
baffle is quite distinct and if a baffle effect is required a Type III design 
should be used. 

In Fig. 15 the curves are carried to great heights of fall. Where the up- 
stream head at the jump is high a deep and long basin is required. The pos- 
sibility of dissipation of energy in the chute, for example by “splitters,” 


a. Proc. Paper 1402, October, 1957, by J. N. Bradley and A. J. Peterka. 

b. Consultant, London, England. 

c. Blench, T.: Punjab Irrigation Research Institute, Government of Pakistan. 
Reports for 1946 (p. 52) and 1947 (p. 21). 
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might profitably be explored in all cases when Fy exceeds 10 or where it is 
found necessary to excavate the basin appreciably below ground level. 

If Fig. 11 is used for determining the basin floor level, the process is one 
of trials when Fj, Dj and D2 are not known until the floor level is deter- 
mined. The level may be directly determined by the use of the diagram 
shown in Fig. I of the Writer’s discussion on Paper 1403. 
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HYDRAULIC DESIGN OF STILLING BASINS: SHORT STILLING BASIN FOR 
CANAL STRUCTURES, SMALL OUTLET WORKS, 
AND SMALL SPILLWAYS (BASIN III)@ 


Discussion by A. Rylands Thomas 


A. RYLANDS THOMAS,» M. ASCE.—The hydraulic design of a stilling 
basin may be undertaken in two stages: the first producing a two-dimensional 
design in which the floor level and length and the arrangement of baffles are 
determined for the worst combination of discharge and tailwater level and 
then tested for other combinations which may occur, and the second expanding 
this design into a three-dimensional basin with wing walls and bank protection 
downstream, taking account of conditions which may produce asymmetry of 
flow. 

The authors’ experimental work and recommendations relate to the first 
stage, in which the designer has to select from a variety of shapes, sizes and 
positions of baffles and end sills. They are to be complimented on their 
methodical attack on the problem, covering a wide range of Froude Number, 
and their results will be of value to designers not only of the works mentioned 
in the title but also of larger river works where surplus energy is to be dis- 
sipated without danger to the structure. In the following discussion the 
Writer has in view this wider field of application. 

Given a set of conditions the designer’s first step is to determine the basin 
floor level without baffle blocks. Probably in most cases in practice the 
quantities known are the discharge, upstream total energy level, tailwater 
level and total energy level, but not V1, Dy, Dg or F1 which depend on the 
floor level to be determined. If Fig. 11 is to be used a process of trials is 
involved. Curves have, however, been devised for direct determination of 
floor level, for example those of the Central Board of irrigation and Power, 
India.(1) The Writer uses curves of the type shown in Fig. I; Ey the differ- 
ence between upstream and downstream total energy levels is known; the 
critical depth De is calculated from the discharge per unit width q, and the 
floor level is found by reading E2/De. and deducting E2 from the downstream 
total energy level. In the same diagram curves are provided for Ej/D¢ and 
Dj/De, the latter representing height of jump which is of use in estimating the 
length of basin required. Dy may be calculated from Fy and De as shown and 
D2 from Eg if required. 

These curves, like those in the Paper, relate to jumps on horizontal floors. 
In many cases a floor with a gentle slope in the direction of flow is advan- 
tageous, in particular where the flow is controlled by gates; should the jump 
begin to move away from the toe of the chute there is more time to take 


a. Proc. Paper 1403, October, 1957, by J. N. Bradley and A. J. Peterka. 
b. Consultant, London, England. 
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remedial action before the jump moves out of the basin and greater protec- 
tion against scour at the lower tailwater levels if the floor is sloping. Ex- 
periments by the Writer(2) indicate that the same curves apply with sufficient 
accuracy to gently sloping floors (the slope tested was 1 in 43) when the floor 
level is taken at its mid-point. 

The provision of baffle blocks or piers in the stilling basin brings three 
benefits: - 


i) Because the blocks create back pressure, the floor may be higher 
with equal margin of safety in stabilizing the jump. 


ii) Because jet expansion is more rapid, the basin may be shorter with- 
out risk of deeper scour downstream. 


iii) Because of the additional turbulence generated in the jump, the down- 
stream flow is more stable, there being less tendency for asymmetri- 
cal flow or large vortices to develop. 


(i) and (ii) bring economies in construction costs. (iii) is of importance in 
the case of canal works in the prevention of bank erosion downstream and in 
the case of large river works in reducing the risk of deep local scour which 
might endanger the structure. 

The value of baffle blocks in enabling a higher floor level to be adopted 
may be seen from Fig. II in which the Authors’ data (Table 4) are shown to- 
gether with those of the Writer.(2) These data are not strictly comparable 
because of differences of relative areas of the blocks and distances from the 
toe of chute. In the case of the Writer’s results the floor had a slope of 1 in 
43 and no chute blocks were provided; the data represent the condition when 
the jump was just held at the toe of the chute. In the Authors’ results the 
floor blocks were nearer the toe of the chute and the tailwater had been 
lowered to the stage when the beginning of the jump had receded from the toe 
and reached the blocks (see Fig. 19); prolongued flow under this condition 
might lead to cavitation damage and dangerous scour downstream. In gener- 
al, however, it can be said that, if the dimensions and location of the blocks 
are reasonably near to those shown in Fig. 17, a reduction of 15 to 20 per 
cent in depth of basin, compared with depth required with no baffles, is 
permissible. 

The Authors would reserve this margin as a factor of safety and design 
the floor level as though there were no blocks. This amount of margin may 
suit many cases but the Writer prefers to adopt a margin in each case suit- 
able to its needs. For example, where ample safety margin has already been 
allowed in the maximum discharge and downstream retrogression of levels 
advantage may well be taken of the reduced depth of basin required when 
floor blocks are provided. A 10 per cent reduction is represented by the 90 
per cent curve in Fig. Il and this may be used for design purposes in appro- 
priate cases. 

The experimental data plotted in Fig. 12 support the conclusion that the 
length of jump is considerably reduced when floor blocks are provided. The 


length of jump is perhaps more directly a function of height of jump, Dj =D2- 


Dj, than of Dg. A suitable length of basin might be 3 Dj which is easily cal- 
culated from the curve in Fig. I. The minimum length of basin required, 
however, is not necessarily the same as the length of jump as measured by 
the Authors. It is the minimum length with which dangerous scour will not 
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occur downstream under any possible flow conditions. For this aspect of de- 
sign scour tests are required. 

The Authors’ finding that the most effective floor blocks are those with 
vertical upstream faces confirms the same conclusion reached by the Writer 
as a result of model tests with blocks of various shapes and sizes.(2) Blocks 
of total height equivalent to 2m. of which 1.5m. was vertical face and the top 
0.5m. was sloped back at 33° were found to be no better than rectangular 
blocks 1.5m. high without the sloping top. There are, however, other con- 
siderations which may decide the shape of block in a given case. Where high 
velocities may be encountered precautions must be taken against cavitation 
(Fig. 4 of Berryhill’s paper(3) is surely a warning) and streamlining may be 
necessary. If fish are to be protected from damage the corners of blocks 
should be chamfered or rounded. In such cases the effective area of the 
blocks will be less than the projected area and generally reliance can only 
be placed on the area of vertical face. 

The position of the blocks on the floor is a matter on which guidance is 
needed. The nearer the blocks are to the toe of the chute the higher the 
velocity and hence the greater the back pressure created, but if they are too 
near the jet “hurdles” over the blocks which are then ineffective. This was 
found by Inglis and Joglekar(4) who used a continuous baffle and the same 
effect has been observed by the Writer with rectangular blocks and by Harle- 
man(5) with streamlined blocks. In none of these tests, however, was the 
range of F, sufficient to establish a general rule for design, such as that of 
the present Authors that the distance from the toe of chute to floor blocks 
should be 0.8 Dg. As dimensions depend on the size of the work it is con- 
venient to express all linear dimensions as ratios of the critical depth Dc 
which is a suitable linear parameter of the discharge per unit width. The 
optimum distance ratio X/Dc (where X is the distance from toe of chute to 
blocks) is clearly a function of Fy and probably aiso of hg/De. It does not 
depend on Dg though it is possible that if expressed as a ratio of Dg it is con- 
venient for design purposes and sufficiently near the optimum value for 
practical purposes. This, however, cannot be judged from the results pre- 
sented. A curve representing the Authors’ recommendation is shown in 
Fig. II together with the ratios pertaining to the basins tested, vide Table 4. 
These are divided into three classes according to the results obtained and it 
will be seen that the best results were given by blocks relatively nearest the 
toe, though no other decided pattern emerges. 

With high velocities, cavitation imposes a limit, so that in a large work 
blocks must be placed at a greater relative distance from the toe than in 
smaller works with equal F;. Some flexibility in design rules is therefore 
required and it would be of value if the Authors gave the results of tests on 
which they based recommendations, in particular tests made with blocks at 
various distances from the toe of chute, other conditions being the same. 

The optimum height of the floor blocks is related to the depth of jet at 
their position. It therefore depends on their distance from the toe of chute, 
the greater the distance the higher the blocks required for the same effect, 
but increasing the height beyond the optimum is not likely to increase their 
effectiveness and may result in undesirable surface disturbance at lower 
stages. Though the Authors state that the height of the blocks should vary 
with Fy and in Fig. 18 the ratio hg/D, is shown to increase with Fj, the true 
relative height of blocks in their tests was almost constant, as may be seen 
from Fig. Ill. The Authors’ rule for block height could be stated in the form 
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h3 = 0.5 De. This, however, applies only to the authors’ standard design 
where X = 0.8 Do. If X is varied hg should be modified, within limits. 

In the absence of data in regard to scour downstream there is insufficient 
evidence for the conclusion that a second row of floor blocks is unnecessary 
and the Writer feels that the second row should not be so lightly written off. 
Besides constituting a reserve in case any blocks in the front row sustain 
damage, the second row may well reduce bed scour and have the effect of 
creating turbulence which helps to stabilize the flow downstream. Where 
streamlined blocks are used to avoid cavitation, a second row may be re- 
quired to provide the desired back pressure. 

The provision of chute blocks is recommended by the Authors but no data 
are produced to indicate their value where floor blocks also are provided. If 
any tests were carried out the results would be of interest. 

The end sill shown in Fig. 17 is an almost essential feature to reduce 
scour at the end of the floor. Deep scour may still occur further downstream, 
where however it is less dangerous. In model experiments by the Writer (re- 
sults unpublished) a continuous sill with two rows of rectangular blocks just 
upstream (these blocks were of the same height as the sill and were distinct 
from the baffle blocks nearer the chute) were found to be more effective in 
reducing scour than end blocks alone or a Rehbock type sill alone. But the 
continuous sill alone was almost as effective as the sill-cum-blocks in re- 
ducing scour at the end of the floor and, though less effective in respect of 
scour further downstream, is probably sufficient for most cases. 

The height of sill is shown in Fig. 18 as a ratio of Dj. It should more 
properly be considered in relation to Dg (as in Fig. 14) because the height 
required is determined by the flow at the end of the jump irrespective of Fj 
or Dj which represent conditions at the a of the jump. In India, 
where the end sill is also known as a deflector, (4) the height is generally 
1/12 to 1/10 of Dg and this has proved satisfactory. The Authors’ sills are 
in some cases within this range and in others somewhat higher. 

As the purpose of a stilling basin is to afford protection against scour, 
the value of the results would be greatly enhanced if supplemented by tests 
with a mobile bed downstream to indicate comparative depths of scour. In 
this connection mention may be made of the research into the design of canal 
falls with stilling basins which has been in progress for several years in 
India. Observations have been made both in models with mobile beds and on 
full-scale structures and accounts of these have from time to time been 
published. (6) 

In this paper a bold attempt has been made to reduce the design of the 
Type III stilling basin to a routine procedure of determining dimensions 
based on two main parameters, Dg and Fj. In the design of complete struc- 

tures it is not always desirable to adhere rigidly to a standard design and 
flexibility is required to meet local conditions. With so many features to be 
determined it is understandable that the detailed results which led to the 
adoption of the recommended dimensions have not been presented, yet with- 
out these it is impossible to see whether in fact the proportions are those 
most suited to requirements or to judge the effect of deviation. Perhaps the 
Authors will be able to include some of the more important supporting data 
in their reply to the discussion. Even with the data already given, however, 
this Paper will enable an effective initial design to be produced and much 
time to be saved in the experimental work. 
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HYDRAULIC DESIGN OF STILLING BASINS: STILLING BASIN WITH 
SLOPING APRON (BASIN V)@ 


Discussion by A. Rylands Thomas and Earl J. Beck 


A. RYLANDS THOMAS,» M. ASCE.—Sloping aprons offer advantages over 
horizontal aprons in operational control and the numerous experimental re- 
sults presented in this paper are a welcome addition to the rather few series 
of data previously on record. 

The methods of design advocated by the present and previous authors are 
concerned with the determination of tailwater depth at the end of the jump and 
this requires an assumption of jump length, a dimension which is not pre- 
cisely known. A more direct method, not involving any assumption of jump 
length, is illustrated in Fig. A. Where the initial information includes the 
discharge per unit width, q, and the upstream and downstream total energy 
head levels, and therefore Ey the loss of head in the jump, it is a simple 
matter to calculate Ey /D;, where Dg is the critical depth, and read E9/Dc. 
from the appropriate curve for horizontal or sloping floor in the diagram. 

In this diagram, E2 is the downstream total head above the point on the floor 
vertically below the beginning of the jump, and the required floor level at this 
point is therefore easily determined. 

The curve for the horizontal floor is calculated by the well established 
formulas based on the pressure—momentum theorem. The curves for slop- 
ing floors are drawn by eye from the Authors’ data (omitting those of Flume 
A which had expanding sides) and the Bakhmeteff and Matzke data relating to 
the slopes of 0.046 and 0.070. 

Apart from the experimental scatter of points within each series of data, 
the results of the series are not completely consistent among themselves. 
The points relating to Flume D lie above those of Flume B for the same 
Slope (0.10) and the latter lie rather higher than would be expected from the 
positions of the curves for 0.05 and 0.20. Also the points for the range of 
slopes between 0.16 and 0.263 lie so close together that one curve (0.20) has 
been drawn for all these. This may indicate an optimum slope of about 0.16 
beyond which less advantage is obtained in depth reduction, or it may be a 
result of a particular characteristic of the channel or experimental technique. 
There does seem to be a limit beyond which the flow pattern changes and 
this may occur between the slopes of 0.20 and 0.35, depending on Fj. There 
is considerably more scatter in this range, as evidenced by the data of 
Hickox for a slope of 0.336, possibly indicating instability of flow pattern. 


a. Proc, Paper 1405, October, 1957, by J. N. Bradley and A. J. Peterka. 

b. Consultant, London, England. 

c. Hickox, G. H., Discussion on “The Hydraulic Jump in Sloping Channels,” 
by C. E. Kindsvater, Transactions ASCE., Vol. 109, p. 1141, 1944. 
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Whether or not future experimental work will confirm that the curve 
shown for the slope of 0.20 is too conservative and is more applicable to a 
slope of 0.16, it is suggested that in the meantime the results for slopes up 
to 0.20 are sufficiently consistent for the curves to be used for practical de - 
sign. In view, however, of the scatter a safety margin of 10 per cent of E2 
is desirable in determining the floor level, i.e. the floor should be this 
amount below the level indicated by the diagram. 

Air entrainment can have significant influence on the jump profile and 
tailwater level. It introduces a scale effect which may be of importance 
where the upstream velocity is high. It would be of interest if the Authors 
would state whether appreciable air entrainment was observed in any of their 
tests, and its effect. 

In discussing the comparative benefits of sloping and horizontal aprons 
(page 28) the Authors suggest that instead of adopting a sloping apron a 
slightly lower horizontal apron be used. This would increase the cost, but 
also the chief advantage of a sloping apron would not be realised, that of a 
visual warning during gate operation. With even a gently sloping apron, when 
it is seen that the jump has begun to recede from the toe of the chute there is 
still time to control the jump before it moves off the apron. With a horizon- 
tal apron this is not so, and damage may occur due to scour and the tail end 
of the apron before the jump can be controlled. 

Case B is more complex than Case D in that an additional variable is in- 
troduced. If points representing the data of Table 9 were plotted in Fig. A 
they would, however, lie between the curve for a horizontal floor and the ap- 
propriate curve for the sloping floor, in positions approximately consistent 
with the relative proportion of jump length on the sloping floor. This appears 
to provide a suitable basis for design, points being interpolated according to 
the ratio of the jump length on the slope to 6E9. 


The authors have not suggested the use of baffle blocks on sloping aprons. 
Such blocks are, however, as useful in this case as on horizontal aprons in 
reducing the depth and length of stilling basin required. Provided the blocks 
are safe against cavitation damage, the most economical design of stilling 
basin in most cases with moderate upstream head is probably one with a 
gentle slope to the floor and provided with floor blocks and end sill. To de- 
termine the length of basin required to afford adequate protection against 
scour, tests are needed in models with mobile beds, and it is hoped that 
future work on these lines will include scour tests with sloping aprons. 


EARL J. BECK, ! A.M. ASCE.—Application of their experimental data to 
existing apron spillway designs helped the authors substantiate the practical 
value of their recommended design method. 

To extend the range of their comparisons, to spillways of lower head and 
lower Froude numbers, the three existing spillways shown in Fig. 1 were 
re-examined by this same method. These spillways, each of similar sloping 
apron design and founded on sand, were designed with the aid of hydraulic 
model studies.2,3 Although they have low Froude numbers, they have been 


1. Harza Eng. Co., Chicago, IL. 


2. Model Study of the Spillway for the Santee River Dam—US Waterways, Ex- 
periment Station, Vicksburg, Mississippi, 1944. 

3. Model Tests of Spillway for Petenwell and Castle Rock Dams—Arno T. 
Lenz—Engineering Experiment Station, University of Wisconsin, 1949. 
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operating successfully for many years. The pertinent data for these spill- 
ways is given in Table 1. The results plotted on Fig. 2, further substantiate 
the authors findings. 

The experimental information furnished by the authors would have been 
very helpful in each of these cases for establishing the most advantageous 
apron arrangement and thereby eliminating tests of many less satisfactory 
arrangements. However, model test would still have been desirable if not 
mandatory for furnishing additional information which at best the ex- 
perienced designer can only approximate. For example, model tests of the 
spillways mentioned above were used to determine the affect and critical 
limit of tailwater recession due to downstream degradation, to determine 
length of training walls downstream of apron required to prevent erosion of 
adjacent fill sections, to establish erosion patterns for determining safe 
depth of footings for training walls, to determine maximum uplift pressures 
on apron and walls for varying discharge and tailwater, and to investigate 
special operating conditions. 

These spillways were found to be a very economical design for the sand 
foundations. However, due to the erodability of the foundation material it is 
imperative that such spillways be operated with symmetrical flow. For this 
reason a separate regulating bay was provided in each case to permit small 
releases without operating the entire spillway. This has the added advantage 
of allowing tailwater to rise before opening the main gates. The regulating 
bay was provided with a longer apron and longer training walls. 

Many structures have special problems which can be solved best with the 
aid of a model. In these cases the paper will be a valuable aid to the de- 


signer in establishing preliminary designs for verification by model tests. 


For cases where model tests will not be made this paper will be invaluable 
to the designer. 


i 
7. 
~ 


$279 LE $6°O 6€ 82 
SSS Ib $6°0 80T 


(42) (92) (sz) (€z}(22) (tz) (oz) (61) 
2a faq 2q *33 °998/3F °35 
ML widep Ta Tq y3prm utseq urtseq 
b M 
UIA 


April, 1958 


$9 $*92 000‘'F9T 68 $2 9°26 S2°O UTSUOSSTM AION 
€T S6 S€ LOT 000‘'9tI 6€ 696 SET 
S°6 0000S ZI €9 GL - PUTTOTeD *g Iado0D aaqueg 


(st) (pt) (et) (zt) (tt) (or) (6) (8) (2) (9) (Ss) (z) (1) 
uoaide uoide yjdep *xeur uoade uoide uoirde urep 
Sutdojs ML ML DO uo jo jo puo jo jo jo 
jo jo peofy pus uMmop pue dn adotg 
q T Wea 


*T 


SNOUdV ONIGOTS HLIM SNISVE ONITTILS ONILSIXG 


1616 -44 HY 2 Pe 
a — Now 
“ee 
ooo 
= 
as 
a 
Ob c — 
oon 
<2 
ee. 
ow 
4 
N= 
owu- 
© 
w 
ANN 
| 
nal 


DISCUSSION 


SANTEE COOPER 


PETENWELL 


CASTLE ROCK 


FIG. | 


TW 76.5 


EL.4725- 


J 
ASCE 1616-45 
EL.75 
== 
| | [EL.22.5 
ELI35 
| EL.72 | | 
| | 
EL.92.5 
65' 
EL.42.5 
| 


April, 1958 


| 
T 


wo 


1616 -46 HY 2 
T Too 
OO ER 
| 
| 
COC 
EEE 
O I 2 1 3 4 5 
De 
FIG. 2 = 


ASCE 1616-47 


THE HYDRAULIC DESIGN OF STILLING BASINS* 


Discussion by W. H. R. Nimmo and Walter Rand 


W. 4H. R. NIMMO,! M. ASCE.—This comprehensive group of papers may 
well serve as a manual of design wherever conditions are suitable for one of 
the type of stilling basins proposed by the authors. In the case of two dams 
in Queensland, designed under the direction of the writer, however, such 
conditions did not exist. 

Somerset Dam on the Stanley River is a straight gravity structure for 
water supply and flood mitigation. A central spillway, having a crest 130 ft. 
above the horizontal apron, is divided by piers into eight openings, each 26 ft. 
long controlled by sector gates 23 ft. high. The downstream face of the spill- 
way is joined to the apron by a curve. There are no chute blocks. Overall 
length of spillway including piers, is 253 ft. 4 in. Below each bay of the spill- 
way is a 12 ft. x 8 ft. low level sluice, which has a curved profile and dis- 
charges on to the apron close to the foot of the spillway. At either end of 
the spillway are valve chambers, each containing two 7 ft. 6 in. cone type dis- 
persing valves. The width between the 40 ft. high parallel sidewalls of the 
stilling basin is 340 ft. The sluices and valves are required to rapidly lower 
the reservoir after floods which temporarily submerge a railway bridge. To 
avoid excavation of sound rock—the river bed was originally only 150 ft. 
wide —the sluices were placed below the controlled spillway although when 
construction commenced in 1935 there appeared to be no precedent —at least 
not in this country—for such an arrangement. The maximum discharging ca- 
pacity, with a head of 23 ft. on the spillway, is spillway 100,000; sluices 
50,000; and valves 10,000 second feet respectively. At low flows, tailwater 
level is insufficient to induce a satisfactory jump. At high flows, tailwater 
may exceed 50 ft. and a stilling basin then becomes almost unnecessary. The 
problem was to dissipate the energy of the 70 ft-sec. jets from the 12 ft. x 
8 ft. sluices acting singly, in groups, or all simultaneously. 

Model experiments showed that it was essential to install a continuous 
barrier sufficiently massive to not be displaced by the jets. As constructed, 
the stilling basin has a length of 196 ft. There is a baffle wall 10 ft. high by 
18 ft. wide, 106 ft. from the foot of the spillway slope and a 4 ft. x 4 ft. con- 
tinuous sill at the downstream end. The model experiments showed that under 
the worst conditions the velocity of water leaving the basin will not exceed 14 
ft-sec. which is regarded as satisfactory for the rock in the river bed. 

Assuming that the combined maximum discharge of the spillway and 
sluices, 150,000 sec. ft., is spread uniformly over the 253 ft. 4 in. length of 
the spillway, the discharge per foot of width would be 593 ft. From the charts 


a. Proc. Papers 1401-1406, October, 1957, by J. N. Bradley and A. J. Peterka. 


1. Civil and Hydraulic Engr., Brisbane, N. 8, Queensland, Australia 
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given in the papers, it is estimated that the length of a type II basin would be 
219 ft. If the action of the Somerset Dam basin is less elegant than that of a 
Basin Il, it is effective and costs are about the same. 

Tinaroo Falls Dam nearing completion on the Barron River is a gravity 
structure with spillway crest 131 ft. above the level surface of the apron. 
Distance between the side walls of the spillway is 250 ft. With a head of 12 
ft. on the crest, the maximum outflow will be 41,000 sec. ft. Model experi- 
ments showed that at no rate of outflow was tailwater high enought to de- 
velop a jump even with dentated sills. This could have been corrected by 
lowering the floor of the basin but this would involve removing sound rock 
below the elevation of the foundation of the dam itself and was considered to 
be inadvisable. The stilling basin as constructed includes a continuous baf- 
fle wall 8 ft. high by 12 ft. wide, 68 ft. from the intersection of spillway 
slope with the apron and a continuous sill, 38 ft. from the baffle wall. The 
model showed that the basin, though turbulent in action, functioned satisfac - 


torily. It is no longer than a type II basin proportioned by the rules given by 
the authors. 


WALTER RAND,! A.M. ASCE.--The effort of Messrs. J. N. Bradley and 
A. J. Peterka to supplement the 23-year experience of the U. S. Bureau of 
Reclamation by numerous new experiments, and to publish the results with 
the idea to *. . . . emphasize practical design .... so that dependable stilling 
basins may be designed without the need for exceptional judgment or exten- 
Sive experience.”, as declared in the foreword of the paper, should be appre- 
ciated by any engineer who deals with the design of stilling basins. These 
papers certainly belong among the best that have been written on this subject. 
However, because of the approach, quoted above, what is offered is not so 
much a generalization of the design of stilling basins, as intended by the 
authors (Paper 1401, page 2), but a standardization of the design. In stress- 
ing structural standards, the hydraulic features of the design have been less 
emphasized, even as much hydraulic data have been given in the form of 
Summaries and tables. Such presentation makes it difficult to see why cer- 
tain dimensions of a proposed basin are the most appropriate, and how they 
depend on the flow pattern in the given basin. Some of the proposed basins of 
finished design may be confidently adopted without much judgment or ex- 
perience, and without fully understanding their hydraulic background. If, 
however, due to particular reasons, a somewhat different design is contem- 
plated, it would be difficult to predict the resulting changes in the flow pat- 
tern, and to determine which structural changes may be acceptable. 

In the writer’s opinion, a general design method should consist of a set of 
hydraulic and structural design principles on the basis of which an individual 
design solution can be found for a particular stilling basin. Ordinarily, de- 
sign does not mean that complete structural solutions should be copied. How- 
ever, the hydraulics of stilling basins is not yet so far developed, and the 
typical design solutions, offered by the authors, are of great immediate value. 

As, in any case, it is an advantage to understand better the hydraulic back- 
ground of the proposed stilling basins, it may now be propitious to discuss 
certain aspects of it, One of the basins, for example, is provided with a 
dentated end sill (Basin II). The length of the forced hydraulic jump, pro- 
duced by this sill, is considerably shorter than the natural hydraulic jump. 


1. Asst. Prof. of Civ. Eng., The City College, New York, N. Y. 
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Is it exactly of the same length as the basin? Could this forced jump, and 
consequently the basin, be even shorter? How long is the basin relative to 
the shortest possible forced hydraulic jump? Why is the chosen basin 

length the best? There is an even shorter basin (Basin III). How are the 
dentated sill of Basin II and the baffle piers of Basin III in their action hy- 
draulically related? In connection with a general study of stilling basins, the 
writer has investigated the properties of a forced hydraulic jump, that was 
produced by a continuous, vertical end sill.“ The results seem to agree 
qualitatively with the effect of a dentated sill or baffle piers, and some of the 
questions above may be tentatively answered. 

Fig. A summarizes observations on lengths of basins. The ratio of the 
distance L, of the sill from the entrance of the stilling basin to the entrance 
depth Dy is given as a function of the Froude number F. The shortest pos- 
sible basin length, given by min.Ls/Dj, corresponds to the shortest physical - 
ly possible forced hydraulic jump. The maximum relative distance of the 
sill from the entrance is given by max.L,/D,. In this case, the natural hy- 
draulic jump of the relative length L/D, occurs in the basin. Between 
min.L, and max.Lg any particular sill of a position Lg and height hg should 
be accompanied by a particular tailwater depth to avoid a sweepout of the 
jump. Any distance Lg/D, can be expressed by a coefficient K with a value 
between zero and one. The forced jump that corresponds to K = 0 to K = 0.2 
is an effective energy dissipator. However, the flow conditions downstream 
of the sill would not be acceptable for an unprotected channel bottom. At 
K~- 0.2 to K = 0.4 an effective forced jump is accompanied by satisfactory 
flow conditions downstream of the sill. If the basin length is extended over 
what corresponds to K = 0.5, there would be only little improvement in flow 
conditions, and the corresponding stilling basin would be of considerable 
length. At K = 0.8 the basin length is equal to the natural hydraulic jump. 
However, as there is an end sill, the hydraulic jump that actually occurs is 
not a natural one. Only at K = 1 is the basin of sufficient length to produce a 
natural hydraulic jump upstream of the sill. 

Fig. A has been completed by the drawing of lines, representing the dis- 
tance from the entrance of the basin to the end sill or baffle piers, corres- 
ponding to the discussed stilling basins of type II, II, and IV, and to the 
shortest existing stilling basins of type I, and to the basin length proposed 
by Forster and Skrinde.3 The plot impresses one with its diversity. All 
these are acceptable design solutions for basin length, provided the height of 
the tailwater and of the sill are adjusted to the given length. These solutions 
are not contradictory as an analysis of flow conditions will show. Basin II, 
that corresponds to K = 0.4 to K = 0.5 makes use of an effective forced hy- 
draulic jump accompanied by a relatively tranquil flow pattern downstream 
of the sill. This is an excellent design solution. The shortest existing type I 
stilling basins correspond approximately to K = 0.1 to K = 0.2. Acceptable 
design solutions are possible also in this range. The corresponding forced 
hydraulic jump is an effective energy dissipator, and if the height of the sill, 
and the given tailwater depth, produce acceptable flow conditions downstream 
of the sill, depending on the erodability of the channel, the basin could still be 
satisfactory. 


2. “An Approach to a Generalized Design of Stilling Basins” by Walter Rand, 
The N. Y. Academy of Sciences Transact. Vol. 20, No. 2, 1957. 

3. “Control of the Hydraulic Jump by Sills” by J. W. Forster, and R. A. 

Skrinde, ASCE. Transactions, Vol. 115, p. 973, 1950. 
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The baffle piers of Basin III, corresponding to K = 0.0 to K = 0.05, pro- 
duce a jump that is practically the shortest physically possible. These piers 
are close to a critical position (min. Lg), and also their height is close toa 
critical value as shown later. The notion of the authors that there is no 
a ..... intention to give the impression that the position or height of 
the baffle blocks are critical” (Paper 1403, p. 16) results from their different 
experimental approach. Of course, the basin has to be extended beyond the 
baffle blocks to protect the channel against harmful erosion. Because the 
forced jump is produced by the baffles, the end sill serves only to improve 
flow conditions downstream of the piers, and it can be smaller. 

Fig. A also shows that Basins II and II] become relatively short (smaller 
K values) at smaller Froude numbers, and no solution is therefore possible 
for F< 4. This leads to Basin IV. Would it not be possible to have a con- 
tinuous variation from the length of Basin II to the length of Basin IV? 
Forster and Skrinde apply such a trend. However, their basin is relatively 
long at larger Froude numbers. 

This analysis shows that the range and properties of the forced hydraulic 
jump can be used to compare the relative length of different stilling basins. 
It also shows that there is no basin length that is the only correct one. Many 
choices are possible in a certain acceptable range, provided the height of the 
sill and tailwater, and the erodability of the channel are correctly related. 
Some observations will be offered about these relationships. 

The concerted action of the sill and tailwater keep the jump from being 
swept out from the stilling basin. If the problem of erosion downstream of 
the sill is provisorily left out of consideration, then a higher sill could be 
used with a lower tailwater to stabilize the jump in the basin. It is possible 
to define the height of the sill as a function of tailwater, and to determine the 
smallest sill that is sufficient to stabilize the jump without any collaboration 
from the tailwater (in the case of very low tailwater). With this critical sill, 
the required depth of the tailwater depends only on the degree of erodability 
of the channel. 

The authors have another approach. They recommend a tailwater of con- 
jugate depth Dg, that is sufficient to stabilize the jump without any action of 
the sill. This recommendation is partially based on a statistical average of 
the tailwater depth at existing structures, equal to 0.99 Dg. As sometimes 
an unnecessarily high tailwater is naturally present, this average may not 
indicate what is actually needed. 

The height of the sill of Basin II is equal to 0.2 Dg (Paper 1403, p. 7). 
is relatively low, especially at larger Froude numbers, if compared to the 
critical height of the sill, that would stabilize the jump without any action of 
the tailwater. Therefore the lowest tailwater depth that prevents the sweep- 
out is close to D2, and relatively closer at larger F values (Paper 1402, 

Fig. 11). The effect of a high sill is well illustrated by the baffle piers of 
Basin II]. Their height is close to the critical, and a lower tailwater stabi- 
lizes the jump. 

By the limitation of the height of sill to 0.2 Dg, which is an average value 
based on the existing stilling basins, and in no shown relationship to the basin 
length or tailwater, a rather rigid standard has been established. Could not 
more flexibility be provided? 

The writer hopes that these observations further illustrate the hydraulic 
aspects on which the design of stilling basins depends. 
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FRED W. BLAISDELL,* M. ASCE.—Messrs. Bradley and Peterka have 
performed a most useful service to the designer in presenting a number of 
generalized stilling basin designs in this series of papers. For the smaller 
structure, where a special model study is usually considered too costly, these 
papers make available a number of carefully tested designs. For the larger 
structure, where an individual model may be justified, these designs will pro- 
vide valuable basic information if they do not actually eliminate the need for 
the special study. 

The generalization of the results obtained from model studies on specific 
structures, supplemented by additional laboratory research to fill in the gaps, 
is something that has been long talked about but on which little has been ac- 
complished. From these examples, it is hoped that others will see the result - 
ing significant benefits and will make similar analyses of additional struc- 
tures. 

The authors show an equation for the curve representing E,/D1 in Fig. 8. 
They do not mention that this is a theoretical equation derived through the use 
of the energy equations 


the equation of continuity, and Eq. 3. The theoretical equation representing 
the curve of Ey /Ej, derived from the same basic equations, can take the form 


The curve of Ey, /E, in Fig. 8 does not represent Eq. A precisely, so the 
writer has shown it in Fig. 8A. The very good agreement between the experi- 
mental data and the theoretical curves of Fig. 8a is not unexpected in view of 
the good agreement between the theoretical and experimental conjugate depths 
shown in Fig. 5. 

The proportion of the energy entering the hydraulic jump which is lost in 
the jump is seen in Fig. 8 to increase with the Froude number. The advantage 
of high Froude numbers is apparent if a high degree of energy dissipation is 
desired. This is true not only for the hydraulic jump stilling basin but also 
for those basins where the “jump” is assisted by means of blocks, piers or 
sills. 

It is sometimes possible to exercise a degree of control on the Froude 
number. A transition between a chute or a culvert and a stilling basin may be 


a. Project Supervisor, U. S. Department of Agriculture, Agricultural Re- 
search Service, Soil and Water Conservation Research Division, St. 
Anthony Falls Hydraulic Laboratory, Minneapolis, Minnesota. 
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used to spread the flow into a thinner sheet. This may cause some change in 
Vi but the major change will be a decrease in Dj and a resulting increase in 
F; (see Eq. 1). Based on certain assumptions, the writer has shown(1) that 
reductions in the concrete and excavation volumes for chute or culvert outlets 
as high as 62 percent and 41 percent respectively can be achieved. Cost re- 
ductions might be even greater since the net reduction in concrete volume is 
obtained by a large reduction in the more expensive wall concrete and an in- 
crease in the less expensive floor concrete while the net reduction in excava- 
tion is achieved by widening the excavation and decreasing its depth, the 

deeper excavation being more expensive in the usual case. The use of a 
transition to increase the Froude number will frequently be economical for 
culverts, outlet works and chutes, but will probably not be practical for rela- 
tively wide and short structures such as dams. Further advantages resulting 
from an increase in the Froude number are that a larger percentage of the 
energy is dissipated, as shown in Fig. 8, and it may be possible to use the 
shorter Basin III instead of Basin IV or Basin VI. 

When discussing Fig. 9 the authors mention the pulsating action observed 
for Froude numbers between 2.5 and 4.5. The writer made a similar observa- 
tion during the studies leading to the development of the SAF stilling basin(2) 
for a range of the Froude number from 1.73 to 17.3. During these studies it 
was noted that the best performance of the stilling basin was achieved at 
Froude numbers greater than about 4.5. Under the heading, “Practical Con- 
siderations” in Paper 1401, the authors suggest that the 2.5 to 4.5 Froude 
number range sometimes can be avoided “by altering the dimensions of the 
structure.” The writer agrees and suggests a transition as one means of 
making this alteration. 

In Paper 1402 the authors state that wingwalls are not amenable to generali- 
zation. The writer is not in agreement with this broad statement. Three 
generalized stilling basins have been studied under the direction of the 
writer(2,3,4) and in each basin wingwalls having a sloping top and flaring in 
plan have been found to be superior to other types of wingwall in protecting the 
dam fill from scour. They, in combination with a good stilling basin, have 
been found so effective that the writer ordinarily suggests that no riprap be 
used in the downstream channel to prevent bed and bank scour, except possibly 
for that bank scour and beaching caused by surface waves. Experiments(5) 
conducted by T. M. Brett for the Soil Conservation Authority of Victoria, 
Australia verify the effectiveness of the sloping top, flaring wingwall. 

Basin III described in Paper 1403 and the SAF stilling basin are similar in 
appearance in that chute blocks, baffle piers and an end sill are used. The 
basin dimensions, however, are considerably different in some particulars. 
The similarity may be due in part to the fact that the principal features of the 
SAF stilling basin were obtained from a rectangular stilling basin design re- 
ported by Jacob E. Warnock, M. ASCE,(6) a former associate of the authors 
at the Bureau of Reclamation Hydraulic Laboratory. The differences perhaps 
can be ascribed to differences in philosophy of those who developed the pro- 
portioning rules. The authors apparently proportioned Basin III to insure that 
the energy dissipation took place entirely within the stilling basin. From the 
photographs of Flumes B, C, D and E in Figs. 1, 2 and 3, it appears that the 
authors’ tests were performed in fixed-bed channels. The SAF stilling basin 
tests were performed in moveable bed channels where the bed scour was used 
as a measure of the stilling basin performance. No attempt was made to ac- 
complish the entire energy dissipation within the stilling basin if the scour or, 
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more properly, the lack of scour was such as to not endanger the structure. 
That some energy dissipation takes place outside the SAF stilling basin is 
shown in Fig. A where it can be observed that a standing wave may occur 
downstream of the basin exit. Since this standing wave is on the surface, the 
end sill protects the stream bed, and the basin performance is satisfactory, 

it is not felt necessary to lengthen the basin and thereby increase its cost. 

This difference in philosophy and experimental methods has resulted in 
the SAF stilling basin being much shorter than Basin III]. The comparison is 
made in Fig. 12A. Relatively longer SAF stilling basins are required at low 
Froude numbers. This is in contrast to the relatively shorter lengths re- 
quired at low Froude numbers for Basins I, D and III. The writer feels that 
the relatively shorter basin lengths required for the SAF basin as the Froude 
number increases reflects the greater percentage of energy dissipation that 
occurs naturally in the hydraulic jump with increasing Froude numbers, as is 
shown in Fig. 8A. 

The authors locate the Basin II] baffle piers 0.8D2 downstream from the 
chute blocks. The basin length is given in Fig. 12 as from 2.2D9 to 2.75Do. 
The baffle piers are thus located from 1/2.75 to 1/3.34 times the basin length 
downstream from the chute blocks. The corresponding distance recommended 
for the SAF stilling basin is 1/3. The agreement is close. 

A difference is that the baffle piers in Basin III do not have to be staggered 
with the chute blocks. Staggering is recommended for best performance of 
the shorter SAF basin. With this one exception, the writer emphatically en- 
dorses the authors comments concerning baffle piers. 

The method used to determine the elevation of the stilling basin floor illus- 
trated in the “Application of Results (Example 3)” deserves emphasis. It has 
been the writer’s experience that many designers consider that the basin 
floor or end sill elevation should be at stream bed level. This is not neces- 
sarily the correct elevation. The correct elevation is determined by sub- 
tracting the required basin tailwater depth from the elevation of the natural 
tailwater downstream from the stilling basin. The authors have correctly 
shown that the floor elevation should be checked at maximum and lesser dis- 
charges. 

The recommended width and spacing of the baffle piers is shown in Exam- 
ple 3 to be 25.5 in. The basin width is computed as 50 ft. from the total dis- 
charge and the discharge per unit width. This results in about 23-9/17 
spaces. An even number of spaces is required if a half space is to be located 
next to each sidewall. There will be 24 spaces if the baffle pier width and 
spacing is 25 in. A baffle pier 25 in. wide can then be centered on each odd 
numbered 25-in. distance from the basin wall. 

In the case of the SAF stilling basin, the baffle piers are specified to occu- 
py from 40 percent to 55 percent of the basin width. One reason for giving 
this latitude was to permit the designer to select sizes that would eliminate 
fractions of inches. Presumably some such latitude is possible in Basin III, 
but the writer agrees that, for best performance, the piers should occupy 
about 50 percent of the basin width. 

In the section headed “Effect of Slope of Chute” in Paper 1405, mention is 
made that “The baffle piers on the floor [ of the Type III basin] tend to alter 
the asymmetrical jet, resulting in an overall improvement in operation.” The 
writer observed similar results for the SAF stilling basin where the depths at 
the basin entrance were non-uniform. In one case, the circular pipe outlet 
discharged onto a diverging chute. At the basin entrance the jet was high in 
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the center and there were standing waves at each sidewall. In spite of this 
very poor depth distribution, the scour pattern at the basin exit was complete- 
ly normal. This is ascribed to the action of the baffle piers in correcting for 
the non-uniform depths at the basin entrance. Attempts should be made to 
secure uniform depths and velocities in the approach to the stilling basin, but, 
if they are not uniform, the baffle piers of Basin III or other similar basins 
will help correct the adverse condition. 

No Single stilling basin can be expected to best fit all the many combina- 
tions of discharges, velocities, tailwater conditions, bed materials, sites, 
risks to life and property, and other tangible and intangible conditions that 
must be considered by the designer. The selection of the most suitable 
stilling basin for detailed investigation from those types presented by the 
authors is facilitated by the comparison of the important features of each 
basin summarized in Fig. 47. There are, of course, other types of stilling 
basins that were not studied by the authors, but this in no way detracts from 
the value of their contribution. 
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HYDRAULIC DESIGN OF STILLING BASINS: HYDRAULIC JUMP ON A 
HORIZONTAL APRON AND HYDRAULIC JUMP ON A SLOPING APRON@ 


Discussion by Solano Vega-Vischi 


SOLANO VEGA-VISCHI, ! A.M. ASCE.-—In analyzing the data obtained by 
the U.S.B.R., instead of using the Froude Number as the parameter for relat- 
ing the data, a different approach was attempted, making use of the critical 
depth as a factor in defining the dimensionless relationship of the different 


features of the phenomenon. This relationship is expressed as the ratio ante 


where dp is the depth of the flow entering the hydraulic jump and d, is the 
critical depth for that discharge. 


A correlation between the Froude Number and the parameter - can be 
Cc 
established. In the following equations this parameter will be noted as Xp. 
The Froude Number is defined as: 


gd 


where F is a dimensionless parameter, V and d are velocity and depth of 
flow respectively entering the jump, and g is the acceleration of gravity. 


The correlation of both parameters can further be demonstrated as 
follows: 
Vo 


F = 


a. Proc. Papers 1401 and 1405, October, 1957, by J. N. Bradley and A. J. 
Peterka. 


1. Cons. Eng., Santiago, Chile. 
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In this form it is evident that the parameter X, still retains the same ad- 
vantages as the Froude Number with respect the dynamic similarity of the 


models, and further allows the use of simple analitical expressions in the 
foundamental equations. 


The theoretical aspects for the derivation of the theory using the critical 
depth criteria has been taken from the treatment given by F. J. Dominguez. 2 
This derivation is explained in the following paragraphs. 

Referring to figure 1 which represents a rectangular channel with a slope 
a@, the theorem governing the variation of the linear momentum is applied 
over a projected length above the horizontal axis Y-Y, over which the rate of 
change of momentum and the forces producing that change are equated. 


Notation of symbols used in the following equations (reference to figure 1) 


Vo= velocity in the initial section of the hydraulic jump 

V, = velocity in the final section of the hydraulic jump 

do= depth of the supercritical flow at the initial section 
of the jump 

d, = depth of the subcritical flow at the final section of 
the jump. 

a = angle of the bottom grade with a horizontal plane 

L = length of the jump measured on the bottom grade 

L, = length of the jump measured ona horizontal plane 

g = acceleration of gravity 

q = Tlow of water per unit width of channel, in the same 

units of measurement used ing 
unit weight of water 


The increment of linear momentum projected on Y-Y is: 


qa(Vi-Vo) cosa 


The difference of hydrostatic pressure in the initial and final sections 
over Y-Y is: 


costa (dé - d?) 


In addition to the external horizontal hydrostatic pressure forces acting 
on the water enclosed by the initial and finalsections, the horizontal 


2. F. J. Dominguez.—Curso de Hidraulica.—Santiago de Chile, 1945. 
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component of the reaction caused by the water’s weight on the bottom of the 
channel, should also be included. This component force, projected onto Y-Y 


wf hdl sina 
0 


This integral can be simplified by substituting an average value h, for the 


mean integrated values of h. Substituting L for L, the value of the 


cosa 


integral becomes: 


wh. L, sina 
wheL, tona 


The equation for the hydraulic jump can then be written: 
| 2 2 2 
q(V,-Vo)cosa@ = >wcos a(do-d;) + wheL, tana 


_cos’a 
cosa (d?-d°) + hel, tana 


but LR 2 de in which d, is the critical depth 
2 
Dividing equation (2) by. de : 
2 


(3) (=. )cosa costa 


Noting: 


de 
d, be 


equation (3) can be written: 


(4) (x. x,) cosa > X,;) + HL, tana 


EO 
| 
he L 
= 
| de de 
Gd... h 
do_ 
| 
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If the average value for h, is represented by the arithmetic mean 
| 
> (dot d,) , equation (4) becomes: 


tot 2 2 
(x x cos*a x’) + 5 (Xo +X tana 


(5) + 5(X, cosa)’ > X,L.tana = 


+ $(Xycosa)’ + $Xol,tana 


Since the influence of cosa@ for slopes up to 0.200 is small, the use of 
equation (5) can be restricted to slopes above 0.200. 


The following simpler equation can be applied for slopes between 0.000 
and 0.200: 


X + 5X, L, tana = x,* +X? + 


In applying equation (5) to slopes above 0.200, and considering that the 
tailwater level was maintained horizontal for the test conditions, (reference 
to figure 2), equation (5) becomes: 


cosa , | 
(7) X,Letana = Xo + 3(X,cosa)* + 5 tana 


Length of jump. 


It is stated in Proc. Paper 1401, ‘the length of jump was the most difficult 
measurement to determine. Where chutes or overfalls were used the front of 

the jump was held at the intersection of the chute and the horizontal floor. 

The length of jump was measured from this point to a point downstream i 
where either the high velocity jet began to leave the floor or to a point on the 
surface immediately downstream from the roller, whichever was the longer. 

In the case where the flow discharged from a gate onto a horizontal floor, 

the front of the jump was maintained just downstream from the completed 
contraction of the entering jet. In both cases the point at which the high . 
velocity jet begins to rise from the floor is not fixed, but tends to shift up- 
stream and downstream. This is also true of the roller on the surface. It 
was at first difficult to repeat length observations within 5 percent by either 
criterion, but with practice satisfactory measurements became possible. It 
was the intention to judge the length of the jump from a practical standpoint; 
in other words, the end of the jump, as chosen, would represent the end of the 


| 
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concrete floor and side walls of a conventional stilling basin.” 
Table I includes 94 U.S.B.R. test results for slopes between 0.050 and 
0.280. (Table 8 Proc. Paper 1405). The experimental values for d,, dj, Lj, 
and q were introduced into each side of equations (6) and (7) to check the 
correlation with the equations, and to determine the acceptance of the hypo- 


thesis used for solving W hd [- The computed difference between the two 


sides of the equation for any given set of experimental values was computed 
as a percentage referred to the average value of the two sides of the equa- 
tion (table I col. 14). The average variation for all the tests was 4.45%, in- 
L dicating it to be within the same range of variance usually expected in the 
experimental observations. It appears, then, that the correlations are suf- 
ficiently close to substantiate the mathematical analogy. 

In an attempt to correlate the length of the jump with the degree of slope, 
plottings were made using various parameters and different coordinates sys- 
tems. The plottings where the parameter X, was used in relation to L¢ and 
the slope tengent showed some apparent consistency for values of X, be- 
tween 0.200 and 0.500, but became erratic for values of Xp outside that range. 
Plottings on natural and semi-log scales and the use of others parameters 
still failed to show a consistent relationship. 

When the test values were plotted on log-log paper, using the parameter 
Lc as ordinates and X, as abscissas, the correlation proved to keep some 
definite relationship (figure 4). The test slopes varied from 0.05 to 0.280; 
and the plotting appeared to follow a definite pattern, showing a consistent 
relationship between L; and Xo, irrespective of the slope. The arc of an 
ellipse very suitable approximated the locus of the plotted points and offered 
a simple analytical medium of definition in the form of a logarithmic ex- 
pression. 

The percentage of error of the curve in relation to the actual values mea- 
sured did not exceed about 5.07%, still within the range of variance expected 
in experimental observations. The expression for the assumed curve is: 


Log Le= 1.36VI-(Log 


(Log Le and Log 10Xo are common logarithms) 

This relationship is to be considered valid between Xp, = 0.100 and 0.700, 
which ordinarily includes the practical accepted range for definable hydraulic 
jumps. 

From an academic point of view the function linking Le and Xp might be 

a criticised, since it does not fulfill the condition of having a value of zero for 
the extreme case Xo = 1 (instead of that, Log Le = 0). However, considering 
that this relationship has been assumed to be acceptable only in the range of 
the experiments, namely, for values of Xp between 0.100 and 0.700, and that 

when Xo values approach to 0.700 (F = 1.7) the jump will lose its stability and 
its length will become indefinable, the above mentioned criticism can be over- 
looked, and the function accepted within the aforementioned range. 

Taking the data from the model tests and substituting the length of jump L' 
computed from the curve whose expression is equation (8) (figure 6), new 
values for each side of equation (5) and (6) were computed. In determining 
the value for Xj, the tail water depth dj also was modified to reflect the new 
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length, this new value becoming: 


(L-L)tona (igure 2) 


Similarly to the comparison of the original values, the computed differ - 
ences between the two sides of the equation were compared to the average 
values of the two sides of the equation to arrive at a percentage variation. 
This percentage variation was slightly less than the original values, showing 
an average error of 4.29% (Table I, col. 16). 

As can be noted in the expression for Le given in equation (8), the length 
of the jump is independent of the slope, L, being a function only of the depth 
of flow at the beginning of the jump related to the critical depth of flow. 

In order to substantiate the hypothesis that L, and Xo are interrelated in- 
dependently of the slope, the experiments performed on horizontal grades 
were plotted in the same manner as those for the sloping channels. (figure 
5). Superimposed on these plottings is the curve defined by equation (8). 

The consistency of the curve and plotted points are noted. The test data for 
horizontal grades are tabulated in Table II (from Table I Proc. Paper 1401). 
The difference between the recorded and computed lengths by equation (8) for 
117 tests shows an average difference of 5.86% (Table II, col. 15). 

On the basis of using the value of L,. in equation (8), equations (6) and (7) 
were solved for values of dp and dj and plotted for various values of the slope 
(tanec). The graph is given in figure 7. Having given the depth of either the 
supercritical or subcritical flow, the conjugate depth of the hydraulic jump 
can be determined. Then from figure 6 the value of Lj can be ascertained. 
By trial and error, the position of the jump can be selected which will satisfy 
the initial condition of momentum at the start of the jump, and the tailwater 
condition below the jump. 


Case of jump formed both on the slope and the horizontal portion of the apron. 
(Case B, Proc. Paper 1405) 


As to the case in which the jump starts on the sloped part of the basin 
(figure 3), the treatment has been for the intermediate case between hori- 
zontal and sloping apron. 


From Table 9 in Proc. Paper 1405, values have been taken to form Table 
III of this discussion. In this table the parameter Xo= = has been com- 


puted from the experimental data and appears in column 6. The partial length 
of the jump formed on the slope is noted as 1 as in the original paper. The 
length L of the jump have been obtained from Figure 6. The ratio 1/L is 
shown for each case in column 9, Table II. 

The conjugate depth for do, namely dj, is supposed to consist of two 
terms: one dj, as if the whole jump were to be formed on a sloping apron, 


and a second term dj as if the whole jump were to form itself on the hori- 
zontal apron. 


d,=d}x (I- + ) (9) 
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or divided by d, 


3 


The factor cos‘a@ has been introduced for sloping aprons equal or above 


0.200. 

In columns 10 and 11 are the values of x; and X} taken from figure 7 for 
the different slopes, and in column 15 the resulting conjugate depths. In 
column 16 the actual conjugate depths (tailwater) for each experience. The 
difference in % is the last column of the table. It is shown that the average 
difference in the 103 experiments, for the computed tailwater according with 
this theory and the actual values measured in the experiments is only 4.1%. 
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THE HYDRAULIC DESIGN OF STILLING BASINS: SMALL BASINS FOR 
PIPE OR OPEN CHANNEL OUTLETS—NO TAIL WATER REQUIRED 
(BASIN 


Discussion by John R. Argue i} 


JOHN R. ARGUE.!—The Bradley-Peterka impact-type stilling basin was 
designed for use in Irrigation works where flow from the structure ‘s car- 
ried away by a trapezoidal channel with invert at sill level and sides sloping 
at 1:1. Where the inflow into the structure is carried in a pipe, head in ex- 
cess of 2 diameters (above invert) would be extremely rare. The pipe would 
usually be running partly full (i.e. open channel flow in the pipe). The 
Bradley -Peterka design shows a rigid relationship between maximum dis- 
charge and basin dimensions and no variation in these dimensions in given 
with different erodible materials downstream from the structure. 

A similar problem to that considered by Bradley and Peterka, has con- 
fronted Local Government engineers in New South Wales, Australia. It con- 
cerns outlet structures for road pipe culverts. 

In many parts of the State, arid lands are occasionally subject to heavy 
rainfall, which causes gullying and severe erosion. It is common for a pipe 
culvert during flood conditions, to be subject to a head of 3 diameters above 
the upstream invert. Where discharges of this magnitude are passed onto 
erodible soils, severe damage is caused to both the roadway and the terrain 
downstream from the structure. 

The existing standard structure consists of a headwall and two wing walls 
set at 45° to the centre line. When this structure operates at high discharge 
under open channel flow conditions in the pipe, a fast moving jet passes di- 
rectly from the pipe onto the erodible bed. The structure itself provides no 
energy dissipation—the full force of the jet is taken by the moveable material 
downstream from the structure. 

The severe scour which occurs around this structure is entirely due to the 
strong eddies which form on either side of the fast moving jet. These two 
eddies rotate in opposite directions—the left anticlockwise, the right clock- 
wise —both being “powered” by the centre line jet passing from the pipe. 

A successful design must, therefore, attack this problem at the root cause 
by ensuring that flow onto the moveable bed is not concentrated into a fast 
moving jet, but evenly distributed at a tolerable velocity across as wide a 
front as possible. This is undoubtedly the criterion for the majority of hy- 
draulic structures which discharge onto moveable beds. 

In pipe culvert outlet structures it is most important that the extent of 
scour downstream from the structure be reduced as much as possible. It is 


a. Proc. Paper 1406, October, 1957, by J. N. Bradley and A. J. Peterka. 
1. Engr., Harbors and Rivers Branch, N. S. W. Dept. of Public Works, 
Sydney, Australia. 
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desirable to have the scour hole confined within the boundaries of the road- 
way reserve. Providing that the scour hole does not encroach upon the toe 
of the roadway embankment, the width of scour is immaterial. A squarish 
hole rather than an elongated hole would, therefore, best suit the condition 
required. 

For pipe culvert structures, it is preferable to have the pipe running full 
(i.e. pipe flow). The advantages which pipe flow holds over open channel 
flow in the pipe are threefold: - 


1. Greater discharge for a given head; 
2. Smaller average velocity; 
3. More even velocity distribution across the pipe. 


To investigate the problem, a research programme using a 4" pipe culvert 
model with downstream pipe invert at natural surface level is being carried 
out at the Hydraulic Laboratory of the Public Works Department, N. S. W. 
The object of this investigation is to determine a satisfactory outlet struc- 
‘ture for all pipe culverts with invert at natural surface level, in the proto- 
type range 18" to 72". 

When details of the Bradley-Peterka impact-type stilling basin were pub- 
lished an interesting comparison was made between this American structure 
and that which has evolved, to date, at the P. W. D. Laboratory. The com- 
parison was made on the basis of a 54" diameter pipe, 28 feet long, sloping 
at 1% grade. The Bradley-Peterka structure model corresponding to this 
size was built in accordance with the details given in A.S.C.E. Paper 1406. 
The alternative end sill design was used in order to reduce scour toa 
minimum. 

Comparative scour tests were carried out in a moveable bed of Sydney 
sand—a fine grained, yellow, non-cohesive beach sand. 

The following discussion gives details of the comparison: - 


Testing Procedure: 


To assess the effect of tailwater on the functioning of the structures, two 
sets of tests were carried out as follows:- 


a) Tailwater at 0.5 feet (prototype) above sill level; 
b) Tailwater at 1.0 feet (prototype) above sill level. 


This level was maintained constant, during testing, by a tailgate and a 
point gauge located in the still water clear of the structure and near the toe 
of the embankment. 

Each test was run for a model time duration of 30 minutes, the build-up of 
discharge taking the first five minutes of this period. The 30 minute test 
duration represents approximately 85% stability of scour. 

Both structures were tested with discharges ranging from 115 C.F.S. to 
330 C.F.S. The lower discharge corresponds to the head at which pipe flow 
would commence, i.e. 115 C.F.S. is the maximum open channel flow discharge. 
The latter discharge (330 C.F.S.) corresponds, approximately, to a head of 3 
pipe diameters. 

The maximum design discharge nominated for the Bradley-Peterka 54" 
diameter outlet, is 190 C.F.S. Thus, it was necessary for the test to cover a 


range which is nearly 75% greater than the design maximum for the Ameri- 
can structure. 
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RESULTS 


COMPARISON OF STRUCTURE FUNCTION 


Rating. The most obvious difference between the two outlet structures was 
the variation in discharge. For a given headwater level, the Bradley- 
Peterka structure given approximately 10% less discharge —at high flows— 
than does the P.W.D. design (see fig. 1) this is probably due to the back 
pressure from the water throttled between the hanging baffle and the head- 
wall in the Bradley-Peterka outlet, reducing the head difference through the 
pipe. 


Open Channel flow (in pipe) 
Bradley-Peterka structure: 0< Q< 240 CFS. 

The hanging baffle had no effect until the flow reached 220 CFS. Up to 
this discharge the flow became progressively more concentrated towards the 
centre of the end sill. At this discharge (220 CFS), the hanging baffle began 
to take effect. At 240 CFS, the baffle caused self-priming—it being impos- 


sible to run the model at open channel flow for higher discharges. 


P.W.D. Structure: 0< Q< 330 CFS. 


The P.W.D. structure will run at open channel flow over the whole range 
of heads up to 3 diameters (270 CFS). Upto 120 CFS the flow across the end 


sill is evenly distributed. Beyond this discharge, the flow, at high velocities, 
becomes progressively more aerated at the centre. 


Pipe flow 
Bradley -Peterka structure: 115< Q < 220 CFS. 


The hanging baffle again appeared to have very little function within this 
range. Flow across the end sill was fairly evenly distributed. 


220< Q< 330 CFS. 


At 220 CFS the hanging baffle began to take effect and separation at the 450° 
end walls commenced. Owing to this separation, high flows across the end 
sill tended to become concentrated into a band equal in width to the parallel 
walled section of the outlet structure (see fig. 7). 


P.W.D. Structure: 115< Q < 330 CFS. 


This structure diverges the flow and a steady hydraulic jump and “surface 
roller” occur at all discharges. Flow across the end sill is fairly evenly 


distributed throughout the full range of discharges. Aeration occurs for all 
flows above 270 CFS. 
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COMPARISON OF SCOUR PATTERNS 


1. Maximum Scour Depth (see figs. 2A, 4A) 


« It is clear from these graphs that the separation effect (at the 45° end 
walls), which commences in the 220 - 230 CFS range has fully developed 
by 260 CFS, and that for higher discharges the flow across the end sill is 
concentration into a width approximately equal to the width between the 
parallel sides of the structure. 7 
The P.W.D. design is less effective at small discharges but more ef- 
fective at high discharges. The end sill of this structure is 16.0 feet wide, 
compared with the 17.5 feet width of the Bradley-Peterka structure. If 
the P.W.D. structure were lengthened sufficient to give an end sill width 7 
of 17.5 feet, then it is probable that at low flows (i.e. not greater than 220 
CFS) the two curves would coincide. 
The Australian structure also has the advantage that the maximum 
depth of scour occurs at a moderate discharge (230 - 260 CFS) and that 
excessive overloading will not produce progressively deeper scour. 


. Downstream Extremity of Scour (see figs. 2B, 4B) 


This aspect is of particular importance to roadway culvert outlets, as 
it is highly desirable to have the entire area effected by scour confined 
within the boundaries of the roadway reserve. 

The curves for the Bradley-Peterka structure again show the influence 
of separation, at the 45° end walls, under high discharge conditions. The 
transition range from zero separation to full separation being again 220 - 
260 CFS. 

The P.W.D. design shows a steady increase over the whole range, 
levelling off at a little over 40 feet from the end of the pipe. 


The two sets of curves corresponding to the two tailwater conditions 
are almost identical. 


Maximum Scour at Structure (see figs. 3A, 5A) 


The safety of the structure and of the embankment depends largely on 
this property. 

An important model alteration was made here in order to give a true 
comparison. 

The American structure has an end sill which extends down to a depth 
of 3 feet. With the erodible material used in these tests, it was found that 
after some flows of short duration, the floor of the structure under the sill r 
was undercut. This produced a ‘model effect” scour pattern which made 
true comparison difficult. To prevent this undercutting, the walls of the 
structure were produced to a depth of 6 feet (prototype) below sill level— 
the internal and plan dimensions remaining the same. In this way, scour . 
depths to RL—6.0 could be tolerated. 

The Australian structure was treated similarly. 

The Bradley-Peterka design showed distinct advantages in this aspect, 
though tailwater was observed to effect the relationship between the curves 
of the two structures. 


Comparative photographs of typical scour around the two structures 
are shown in figs. 9A, 9B. 
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4. Scour near toe of Embankment. (see figs. 3B, 5B) 


DISCUSSION 


This aspect is important from the point of view of safety of the em- 
bankment. 

The serious effects of the extension of scour in the upstream direction 
along the side walls of the structure and towards the roadway, depend on 
the following: - 


a) Maximum scour at structure. 
b) The location of deepest scour. 
c) Angle of repose of the erodible material. 


The curves for the two designs lie fairly well together—the deeper 
structure scour in the P.W.D. design being compensated by its greater 
distance away from the end of the pipe. 

The points corresponding to discharges in excess of 290 cusecs are 
doubtful as the boundary walls of the moveable bed (set 4' -6"' apart), 
caused an eddy at the toe of the embankment on each side of the structure. 
This was clearly a model effect. 


Discussion on Comparison Tests. 
The comparison tests were carried out to determine: - 


Whether the Irrigation outlet of Bradley and Peterka could be successfully 
applied to Australian road culverts. 


b) Whether the structure which had evolved from the research programme at 


the P.W.D. Laboratory was superior or inferior to the American struc- 
ture. 


From the viewpoint of first cost, the P.W.D. structure involves ap- 
proximately three quarters of the quantity of concrete in the corresponding 
American outlet. Also, the formwork in the Australian design is greatly 
simplified. 

As far as performance is concerned, it is very difficult to declare that 
either design is superior. 

Comparison of hydraulic models in moveable beds is usually made on 
the basis of “maximum depth of scour.” This is unreliable here as, in 
the Bradley-Peterka structure, the scour patterns within the design range 
(0 to 190 CFS) are quite different from those beyond this range. When 
tested beyond 0-190 CFS., the scour hole tends to become elongated and 
deep along the centre line. The scour holes of the two structures are 
compared for three flows in diag. 6. 

A deep scour hole downstream from a structure is not, in itself, bad 
unless it promotes erosion of material under or around the structure; the 
P.W.D. design, although causing generally a shallower scour hole, always 
has slightly more scour at the structure. 

The downstream extremity of scour, is again related to the shape of the 
scour pattern. The elongated hole which occurs when the Bradley-Peterka 
outlet is subject to high flows, is dangerous from the point of view of 
gully erosion. Small flows, when passed through a previously scoured 
elongated hole, tend to cut a channel along the centre line. This effect is 
much reduced when small flows pass through a “square ended” scour hole, 
as produced by the P.W.D. structure. 
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Scour at the structure and proximity of scour to the toe of the embank- 
ment are both important factors to be considered. The Bradley-Peterka 
structure is superior here i.e. the tendency to scour around the structure 
is slightly less in the American than in the Australian. 

However, in both cases, rock protection would have to be carried out 
thoroughly around the structure. 

Scour along the side walls was due entirely to “slipping” of the non- 
cohesive bed material in order to lie at its appropriate angle of repose. 
The degree to which this failure occurs is due to the type of material in 
which the structure is built. If scour at the two ends of the end sill can 
be prevented or controlled by suitable rock placement, “angle of repose” 
scour will be eliminated. t 

As mentioned in the preamble, the Bradley-Peterka design shows a 
rigid relationship between maximum discharge and basin dimensions—no 
alteration being made with different erodible materials downstream from 
the structure. The P.W.D. design permits flexibility in this respect. If 
the scour velocity of the downstream erodible material is known, then the 
width of the end sill can be adjusted to ensure that this velocity is not ex- 
ceeded. This adjustment is achieved by suitably fixing the length of the 


structure—the side walls remaining at the same angle with the centre line 
for all designs. 


Further Investigation. 


The present investigation, which will also include study on multiple pipe 
culverts, flood gates, causeways and inlet structures, is far from com- 
plete. However, progress to date has been satisfactory. An attempt will 
be made to increase the angle which the side walls make with the centre 
line in order to further reduce the cost of the structure, and also to study 
the effects of different end sill shapes. 
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9A, 9B: Scour Around Structures often Test with High Discharge 
(is Q7260CFS) A: Bradley-Peterka Structure. B: PWD Structure. 


10A, 10B: Operation of PWD Structure with Discharge equal to 240CFS. 
A: Open Channel Flow in Pipe. B: Pipe Flowing Full. 


11A, 11B: Operation of Both Structures with Discharge equal to 330CFS. 
(Pipe Flowing Full) A: Bradley-Peterka Structure. B: PWD Structure. 
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DISCHARGE CHARACTERISTICS OF RECTANGULAR THIN-PLATE WEIRS#@ 


Discussion by Turgut Sarpkaya 


TURGUT SARPKAYA, | J.M. ASCE.—The study of the characteristics of 
flow over an innocent looking weir has become a basic and challenging prob- 
lem for research workers in fluid mechanics, quite apart from the interest it 
holds for its purely practical engineering applications. Therefore, any engi- 
neer who is concerned with the measurement of discharge by the use of thin- 
plate rectangular weir along with a true understanding of the limitations of 
various approximate formulas should welcome this most noteworthy addition 
of experimental data. The precision of the authors’ measurements and the 
wide variation of the parameters involved, excluding the fluid properties, and 
the qualitative physical explanations of the observed characteristic curves of 
discharge coefficient for the raison d’étre of the proposed formula are the 
outstanding features of the paper. 

The authors associate the influence of viscosity and surface tension with a 
fictitious increase in head h, and in width b, in order to obtain a discharge 
coefficient independent of the size of the weir and of the magnitude of the 
head and embody the combined effects of the several phenomena attributed to 
viscosity and surface tension in the quantities k, and k, without fully ack- 
nowledging the limitations which these quantities imposed upon the proposed 
formula. This writer does not wish to imply that a formula developed by the 
adjustments of head and width or by any other method are not valuable as 
long as they satisfy the experimental results within the range they are de- 
veloped for. He would like, however, to see that the development of the em- 
pirical formulas be closely related with the natural variations of the parame- 
ters involved and be interpreted in that manner wherever possible. And at 
all times it should be kept in mind that an empirical equation derived by fit- 
ting a curve to a certain range of data can scarcely be expected to be general, 
or physically sound. This writer would like to remain faithful to the original 
definition of the discharge coefficient and would submit herewith a criticism 
of the authors’ presentation together with an analysis which he believes will 
clarify the problem further. 

The standard form of flow formula is based on the Bernoulli proposition of 
equality of the sum of pressure and kinetic energies in the two sections con- 
cerned. Changes of energy due to fluid friction are allowed for only in their 
effect on the discharge coefficient, which is an experimentally determined 
value. In addition to the fluid friction, the discharge coefficient is influenced 
by, and allows for, other departures from ideal conditions, e.g., the nonuni- 
formity of velocity distribution, contraction of jet, and the mean curvature of 


a. Proc. Paper 1453, December, 1957, by Carl E. Kindsvater and Rolland W. 
Carter. 


1. Asst. Prof., Dept. of Eng. Mech., Univ. of Nebraska, Lincoln, Nebr. 
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the stream lines. The velocity distribution is determined by initial accelera- 
tions, fluid and boundary friction, and the geometrical characteristics of the 
weir. The increasing influence of viscosity shown by the lower Reynolds’ 
numbers influences the coefficient of discharge in three ways: (a) Internal 
fluid and boundary friction is increased, leading to greater resistance to 
shear in the acceleration through the notch, hence lowering of the discharge 
coefficient, (b) The velocity distribution of flow approaching the weir-notch 
is changed in a direction which increases the central velocity in relation to 
the mean. This results in a greater total kinetic energy in the approach flow 
and a lower requirement of pressure energy to accelerate through the notch, 
thus raising the discharge coefficient, (c) Increased wall friction on the up- 
stream face of the weir plate, causing a retardation of the layers close to the 
surface which form, by their momentum, the lower and the side boundaries 
of the nappe. This results in an enlargement of the latter quantities in com- 
parison with the ones in potential flow, hence an increased coefficient. The 
above reasoning is in close conformity with the one given by the authors. 
However, the authors have failed to set forward the limitations of their rea- 
soning. Is it true that the discharge coefficient will increase with viscosity 
and surface tension always? If it is so, then highly viscous fluids such as 
glucose, honey, and molasses would flow faster than water for a given head. 
This cannot be true. What then is the answer to the apparent reversal of the 
foregoing reasoning? The only explanation for this is that for values of 
Reynolds’ number smaller than a certain value the discharge coefficient de- 
creases rapidly. This is already a well-known phenomenon for flow in ori- 
fices, nozzles and venturi tubes, (Fig. 1). The humped shape of orifice co- 
efficient curves is quite comparable with rectangular notchweir discharge 
coefficient curves. Obviously, the predominance of the influences (b) and (c) 
over a limited flow phase for Reynolds’ numbers greater than Rer, creates 
the hump on the coefficient curve. For values of R smaller than Rey, the in- 
fluence (a) is predominant. In the case of a low ratio of area with practically 
zero approach velocity the hump is due to influence (c) alone. Since the flow 
over the weir is free surface flow, the surface tension has additional effect 
on the influences sighted above and makes the analysis of the problem even 
more complicated. However, since the viscosity seems to play a more im- 
portant role, for lower values of Reynolds’ number—or in other words as the 
viscosity increases—the coefficient of discharge will decrease instead of in- 
crease which would be expected. Therefore, the authors’ explanations over 
the influences of viscosity and surface tension should definitely be confined 
to a flow range where the Reynolds’ number is greater than a certain critical 
value. Moreover, for certain values of Reynolds number some irregularities 
can be expected in discharge coefficient due to the transitional instability of 
the flow in the approach channel at the transition between turbulent and vis- 
cous flows. The proof for-the existence of a critical Weber number and its 
experimental verification is rather difficult. Firstly because, for the liquids 
commonly used, the one with a higher viscosity has a lower surface tension 
at a given temperature, and secondly because for a given liquid the increase 
of surface tension with a decrease in temperature, is much slower than the 
increase of viscosity. Therefore, the viscosity becomes always predominant 
and makes it very difficult to obtain a wide range of Weber numbers when the 
viscosity and the flow characteristics are kept constant. 

The equation proposed by the authors as a comprehensive equation of dis- 


charge for all rectangular thin-plate weirs could be written in the following 
form: 
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2 3/2 k, k 2 
3 b h 


in which C, is assumed to be a function of b/B and h/P alone. If the last 


parentheses is opened and higher orders of small quantities neglected, one 
obtains, 


do. 3 


Experiments of various investigators seem to establish the fact that Co 
could be written in the following functional form: 


£(b/8) + e(b/s) (3) 
P 


in which f and g are used as a symbol of “the function of.” This writer pro- 
poses that f and g be expressed in the following form: 


f(b/B) = (4-a) 


e(b/s) (4-b) 


in whicha,8, ¥,4,m, and n are assumed to be constants independent both 
of the characteristics of the weir and of the properties of fluid. Hence, the 
coefficient of discharge becomes, 


2 
The mean velocity at the section of weir notch is given by, 


2 
vVe—- 
3 


The derivation of a general equation that should be applicable to any 
liquid flowing over a rectangular notch weir is aided by a dimensional analy- 
sis of the problem. Although, the Eq. (5) is empirical one, it should remove 
some of the objections to its use for liquids other than those actually tested 
if it is expressed in terms of Reynolds and Weber numbers and verified later 
by experiments. The coefficient of discharge Cg which is defined here as in- 
dependent of the acceleration of gravity is, 
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The influence of the ratio b/h will be neglected in view of the experimental 
results mentioned by the authors. The definitions of the Reynolds and Weber 
numbers for this particular problem is somewhat crucial. Since both the 
Reynolds and Weber number must be related to both the head and the width of 
the weir, and since there is no way of predicting which length is physically 
more significant for each number. In each number, velocity can be replaced 
by a dimensional equivalent /gh . However, it should be kept in mind that if 
the particular function of Eq. (6) were known, this substitution would not be 
entirely correct. But since the particular function is not known, the replace- 
ment is possible. In this paper Reynolds and Weber numbers will be defined 
as follows: 


h 2 yd h 
Vv S/p 


If now b/B is replaced by 7j,h/P by 9, 1/R by 73, and 1/W by 74, one has, 


(8) 


(11, ) (9) 


Dimensional analysis states that the general solution of the problem is in the 
form, 


in which P, r, s, t, q,...... are constants to be determined. Letting a be 
the value of Cq as 71, 79, 3, and 74 approach zero, one has 


r r 8 


1 | 2 2 
M(b/B) (h/P) (h/P) (11) 


92 


wi R W 


R 


If Ey. (11) is compared with Eq. (5), one obtains, 


in which 


3 
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Now the constants in Eq. (4-a) and (4-b) are obtained by plotting the values of 
{(b/B) and g(b/B) versus b/B, obtained by the authors’ Fig. (10). A logarith- 
mic plot of these values gives a perfect straight line. However, here only 
the plots of these values on an ordinary graph paper are shown, Fig. (2). 

The values obtained are given below: 


Bs 9.015 Me 


Fe ~0.0023 Ns. 


Substituting, 


1x10 ft“/sec, 


AA + Je 
at/ft, 


0.003 ft 


7 A? 3 
fe slug/ft 


Equation (12) takes the following final form: 


2043 
7 +0.035 (b/s) ] + [0.0023 + 0.773 (w/e) 


Gell 


This equation holds very good for all values of b/B and h/P except for b/B = 
0.9. For this latter value experimental results gives about 1% greater Co 
values. In the above equation the value of o = 0.587 deserves more attention. 
The theoretical value of o for slots, orifices, sluice gates and rectangular 
weirs is 0.611. 

With the definitions accepted for Reynolds and Weber numbers, the varia- 
tion of kp, and kp with the fluid properties can be obtained from, 


hy k k )x (15-a,b) 


in which the subscripts x and w indicate respectively “for any liquid” and for 
water.” It would be indeed very interesting to determine experimentally the 
variation of k} and k, with viscosity and surface tension. If the Eqs. (15-a,b) 
hold for various liquids tested, then the Eq. (14) could be used for all liquids. 
And only then the figures (11) and (12) presented by the authors can have a 
true and broader meaning. Moreover, the usage of highly viscous fluids will 
remove the uncertainties existing in kp and particularly in kp. Since accord- 
ing to Eq. (15-a) the value of kpy for a liquid 8 times more viscous than wa- 
ter would be 0.012 ft., which is a value sufficiently great to determine more 
accurately by means of experiments. It should be recognized that, after the 


0.0773 
4 
2.43 
| 

> 

| 

2043 h 
0.52 
(14) 
R /; 
| 


1616 -98 HY 2 


April, 1958 
demonstration of the use of R and W as suitable parameters, the treatment 
is purely empirical. It is conceivable that several different functions might 
exist which would satisfy the experimental data equally well. Extrapolation 
of these functions outside the experimental range is unwarranted and should 
be made with caution or not at all. 

In conclusion, this writer hopes that further experiments will be per- 
formed with different liquids to make the proposed equation applicable to 


liquids other than water as well as water, and held over a wide range of b/B 
h/P, R, and W. 
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